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Symmetry-protected topological (SPT) phases are gapped quantum phases with a symmetry, 
which can be smoothly connected to the trivial product states only if we break the symmetry. 
For a given symmetry, we can have many different SPT phases. But how to describe/construct 
those different SPT phases that can not be distinguished by their symmetry? It has been shown 
that different bosonic SPT phases in any dimensions and for any symmetry groups can be de- 
scribed/constructed using group cohomology theory of the symmetry group. In this paper, we 
introduce a group super-cohomology theory which is a generalization of the standard group co- 
homology theory. Using the group super-cohomology theory, we can describe/construct different 
interacting fermionic SPT phases, in any dimensions and for symmetry groups where the fermions 
form ID representations of the symmetry group. Just like the boson case, our systematic construc- 
tion is based on constructing discrete fermionic topological non-linear a-models from the group 
super-cohomology theory. Our discrete fermionic topological non-linear cr-model, when defined on a 
space-time with boundary, can be viewed as a "non-local" boundary effective Lagrangian, which is a 
fermionic and discrete generalization of the bosonic continuous Wess-Zumino-Witten term. Thus we 
believe that the boundary excitations of a non-trivial SPT phase are gapless if the symmetry is not 
broken. As an application of this general result, we construct three non-trivial SPT phases in 3D, for 
interacting fermionic superconductors with coplanar spin order (which have T 2 = 1 time-reversal 
and fermion- number-parity Z% symmetries described by a full symmetry group Z% x Z(). We also 
construct three interacting fermionic SPT phases in 2D with a full symmetry group Z2 x Z£ . Those 
2D fermionic SPT phases all have central-charge c — 1 gapless edge excitations, if the symmetry is 
not broken. 

PACS numbers: 71.27.+a, 02.40.Re 
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states can be transformed into each other via LU trans- 
formations. So all SRE states belong to the same phase 
(see Fig. la), ie all SRE states can continuously deform 
into each other without closing energy gap and without 
phase transition. 

LRE states are states that cannot be transformed into 
direct product states via LU transformations. It turns 
out that, in general, LRE states may not be connected 
to each other through LU transformations. The LRE 
states that are not connected via LU transformations be- 
long to different classes and represent different quantum 
phases. Those different quantum phases are nothing but 
the topologically ordered phases. 

Fractional quantum Hall states 13 ' 14 , chiral spin 
liquids, 15,16 Z2 spin liquids, 17-19 non-Abelian fractional 
quantum Hall states, 20-23 etc are examples of topologi- 
cally ordered phases. The mathematical foundation of 
topological orders is closely related to tensor category 
theory 9,10,24 ' 25 and simple current algebra. 20,26 Using 
this point of view, we have developed a systematic and 
quantitative theory for non-chiral topological orders in 
2D interacting boson and fcrmion systems. 9,10,25 Also for 
chiral 2D topological orders with only Abelian statistics, 
we find that we can use integer if-matrices to describe 
them. 27-30 

For gapped quantum systems with symmetry, the 
structure of phase diagram is even richer (see Fig. lb). 
Even SRE states now can belong to different phases. The 
Landau symmetry breaking states belong to this class of 
phases. However, there are more interesting examples in 
this class. Even SRE states that have the same symmetry 
can belong to different phases. The ID Haldane phase 
for spin-1 chain 31-34 and topological insulators 35-40 are 
non-trivial examples of phases with short range entangle- 
ments that do not break any symmetry. Those phases are 



I. INTRODUCTION 

Landau symmetry breaking theory 1-3 allows us to un- 
derstand a large class of phases and phase transitions. 
In 1989, it was realized that many quantum phases 
can contain a new kind of orders - topological order - 
which cannot be described by Landau symmetry break- 
ing theory. 4-6 Recently, it was found that topological 
orders are related to long range entanglements. 7,8 In 
fact, we can regard topological orders as patterns of long 
range entanglements 9 defined through local unitary (LU) 
transformations. 10-12 

The notion of topological orders and long range en- 
tanglements leads to the following more general and also 
more detailed picture of phases and phase transitions (see 
Fig. I). 9 For gapped quantum systems without any sym- 
metry, their quantum phases can be divided into two 
classes: short range entangled (SRE) states and long 
range entangled (LRE) states. 

SRE states are states that can be transformed into 
direct product states via LU transformations. All SRE 
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FIG. 1: (Color online) (a) The possible gapped phases for a 
class of Hamiltonians H(gi,g 2 ) without any symmetry restric- 
tion, (b) The possible gapped phases for the class of Hamil- 
tonians H aylnln (gi, g 2 ) with a symmetry. The shaded regions 
in (a) and (b) represent the phases with long range entan- 
glement. Each phase is labeled by its entanglement proper- 
ties and symmetry breaking properties. SRE stands for short 
range entanglement, LRE for long range entanglement, SB for 
symmetry breaking, SY for no symmetry breaking. SB-SRE 
phases are the Landau symmetry breaking phases, which are 
understood by introducing group theory. The SY-SRE phases 
are the SPT phases, which can be understood by introducing 
group cohomology theory. 
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TABLE I: (Color online) Bosonic SPT phases in d sp -spatial 
dimensions protected by some simple symmetries (represented 
by symmetry groups G). Here Zi means that our construction 
only gives rise to the trivial phase. Z n means that the con- 
structed non-trivial SPT phases plus the trivial phase are la- 
beled by the elements in Z n . Z 2 represents time-reversal sym- 
metry, U(l) represents (7(1) symmetry, Z n represents cyclic 
symmetry, etc . Also (m, n) is the greatest common divisor of 
m and n. The red rows are for bosonic topological insulators 
and the blue rows bosonic topological superconductors. 



beyond Landau symmetry breaking theory since they do 
not break any symmetry. On the other hand, their very 
existence requires symmetry in Hamiltonian. We will 
call this kind of phases Symmetry Protected Topologi- 
cal (SPT) phases. 33 ' 34 Some examples of 2D SPT phases 
protected by translation and some other symmetries were 
discussed in Ref. 41-43. 

It turns out that there is no gapped LRE state in ID. 11 
So all ID gapped states are either symmetry breaking 
states or SPT states. This realization led to a complete 
classification of all ID gapped quantum phases. 42 ' 43 ' 48 

In Ref. 49,50, the classification of ID SPT phase is gen- 
eralized to any dimensions: For gapped bosonic systems 
in d sp - spatial dimension with an on-site symmetry G, we 
can construct distinct SPT phases that do not break the 
symmetry G from distinct elements in H d+1 [G, (7^(1)] ~ 
the (1 + d)-cohomology class of the symmetry group G 
with G-module Ut(1) as coefficient. Note that the above 
result does not require the translation symmetry. The 
results for some simple on-site symmetry groups are sum- 
marized in Table I. In particular, the interacting bosonic 
topological insulators/superconductors are included in 
the table. 

In Ref. 25, the LU transformations for fermionic sys- 
tems were introduced, which allow us to define and 
study gapped quantum phases and topological orders 
in fermionic systems. In particular, we developed a 
fermionic tensor category theory to classify non-chiral in- 
trinsic fermionic topological orders defined through the 
fermionic LU transformations without any symmetry. 
In this paper, we are going to use the similar line of 
thinking to study fermionic gapped phases with sym- 
metries. To begin, we will study the simplest kind of 
them - short range entangled fermionic phases with sym- 
metries. Those phases are called fermionic SPT phases 
which can be defined through the fermionic LU trans- 
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TABLE II: (Color online) Fermionic SPT phases in d sp - 
spatial dimensions constructed using group super-cohomoloy 
for some simple symmetries (represented by the full symme- 
try groups Gf). The red symmetry groups can be realized 
by electron systems. Here Zi means that our construction 
only gives rise to the trivial phase. Z n means that the con- 
structed non-trivial SPT phases plus the trivial phase are la- 
beled by the elements in Z n . Note that fermionic SPT phases 
always include the bosonic SPT phases with the correspond- 
ing bosonic symmetry group Gb = Gf/Z^ as special cases. 
The blue entries are complete constructions which become 
classifications. Z 2 represents time-reversal symmetry, U{1) 
represents (7(1) symmetry, etc . As a comparison, the results 
for non-interacting fermionic gapped/SPT phases, 44 as well as 
the interacting ones in ID, 45-48 are also listed. Note that the 
symmetric interacting and non-interacting fermionic gapped 
phases can be the SPT phases or intrinsically topologically 
ordered phases. This is why the lists for gapped phases and 
for SPT phases are different. Here a number n means that 
there are n — 1 non-trivial phases and one trivial phase. 2Z 
means that the phases are labeled by even integers. Note that 
all the phases listed above respect the symmetry Gf. 



formations with symmetry. Fermionic SPT phases are 
T = gapped symmetry unbroken phases of fermionic 
Hamiltonians with a certain symmetry. But if we break 
the symmetry of the Hamiltonians, all those symmetry 
unbroken phases can be smoothly connected to the triv- 
ial product states. Note that in spatial dimension, the 
trivial product states can have even or odd numbers of 
fermions. So the trivial product states in OD can be- 
long to two different phases. In higher dimensions, all 
the trivial product states belong to one phase (if there is 
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no symmetry). This is why when there is no symmetry, 
there are two fermionic SPT phases in OD and only one 
trivial fermionic SPT phase in higher dimensions. 

After introducing the concept of fermionic SPT phase, 
we find that we can generalize our construction of bosonic 
SPT orders to fermion systems, by generalizing the group 
cohomology theory 50 ' 51 to group super-cohomology the- 
ory. This allows us to systematically construct fermionic 
SPT phases for interacting fcrmions in any dimensions 
and for a certain type of symmetries where the fermions 
form an ID representation of the symmetry group. The 
constructed fermionic SPT phases for some simple sym- 
metry groups are given in Table II. The rows correspond 
to different symmetries for the fermion systems. 

We note that, sometimes in literature, when we de- 
scribe the symmetry of a fermion system, we do not 
include the fermion-number-parity transformation Pf = 
(— ) in the symmetry group. In this paper and in Ta- 
ble II, we use the full symmetry group Gf which in- 
cludes the fermion-number-parity transformation to de- 
scribe the symmetry of fermion systems. 44 So the full 
symmetry group of a fermion system with no symmetry 
is Gf = Z 2 generated by Pf. The bosonic symmetry 
group Gb is given by Gb = Gf/Z^, which correspond to 
the physical symmetry of the fermion system that can 
be broken. In fact, the full symmetry group is a projec- 
tive symmetry group discussed in Ref. 52, which is a 
extension of the physical symmetry group Gb- 

The columns correspond to different spatial dimen- 
sions. In OD and ID, our results reproduce the exact 
results obtained from previous studies. 45-48 The results 
for 2D and 3D are new. 

Each entry indicates the number of non-trivial phases 
plus one trivial phase. For example, Z2 means that there 
is one non-trivial SPT phase labeled by 1 and one triv- 
ial phase labeled by 0. Also, say, 8 means that there 
are seven non-trivial SPT phases and one trivial phase. 
For each non-trivial fermionic SPT phase, we can con- 
struct the ideal ground state wave function and the ideal 
Hamiltonian that realizes the SPT phase, using group 
super-cohomology theory. 

We believe that each non-trivial SPT phase has gap- 
less boundary excitations protected by the same sym- 
metry that protects the bulk SPT phase. This belief 
is supported by the discussion in the next section, and 
by the fact that the boundary excitations are described 
by a "non-local" Lagrangian which is the fermionic and 
discrete generalization of the bosonic continuous Wess- 
Zumino-Witten (WZW) term. 53 ' 54 

In the rest of this paper, we will first compare the re- 
sults in Table II with some known results for ID and 
for free fermion systems. This will give us some physi- 
cal understanding of Table II. We then define fermionic 
path integral on discrete space-time via Grassmann ten- 
sor network. Next we discuss the conditions under which 
the fermionic path integral becomes a fixed-point theory 
under the coarse-graining transformation of the space- 
time complex. Such a fixed-point theory is a fermionic 



topological non-linear cr-model. The fixed-point path in- 
tegral describes a fermionic SPT phase. We then con- 
struct the ground state wave function from the fixed- 
point path integral. Some simple explicit examples of 
fermionic SPT phases and their boundary excitations 
are discussed. In appendix, we develop a group super- 
cohomology theory, and calculate the Table II from the 
group super-cohomology theory. We also construct the 
exactly soluble Hamiltonian of the fermionic SPT phases 
from the fixed-point path integral. 

II. PHYSICAL PICTURES OF SOME GENERIC 
RESULTS 

Before describing how to obtain the generic results in 
Table II, in this section, we will compare some of the 
generic results with known results for ID interacting sys- 
tems and 2D/3D non- interacting systems. The compari- 
son will give us a physical understanding for some of the 
interacting fermionic SPT phases. 

First, from the Table II, we find that fermion systems 
with T 2 = 1 time-reversal symmetry (or the full sym- 
metry group Gf = Z2 x Z 2 ) can have three non-trivial 
fermionic SPT phases in d sp — 1 spatial dimension and 
three non-trivial fermionic SPT phases in d sp = 3 spatial 
dimensions. 

For OD free or interacting fermion systems with T 2 = 
1 time-reversal symmetry Zj , the possible symmetric 
gapped phases are the ID representations of Z% x Z^ . 
Since the time-reversal transformation T is anti unitary, 
Z2 has only one trivial ID representation. Z| has two ID 
representations. Thus OD fermion systems with T 2 = 1 
time-reversal symmetry are classified by Z2 correspond- 
ing to even fermion states and odd fermion states. 

For ID free fermion systems with T 2 = 1 time-reversal 
symmetry Z£ , the possible gapped phases are classified 
by Z. 55 ' 56 For the phase labeled by n e Z, it has n Majo- 
rana zero modes at one end of the chain. 57 Those bound- 
ary states form a representation of n Majorana fermion 
operators 771 , rj n which all transform the same way un- 
der the time-reversal transformation T. (If, say r)i and 772 
transform differently under T: rji — > rji and rj 2 — > —112, 
ir\\r\2 will be invariant under T and such a term will gap 
out the rji and 772 modes.) Some of those ID gapped 
states have intrinsic topological orders. Only those la- 
beled by even integers (which have even numbers of Ma- 
jorana boundary zero modes) become the trivial phase 
after we break the time-reversal symmetry Zj. Thus the 
free fermion SPT phases are labeled by even integers. 

For interacting ID systems with T 2 = 1 time-reversal 
symmetry, there are eight possible gapped phases that 
do not break the time-reversal symmetry. 45-48 Four of 
them have intrinsic topological orders (which break the 
fermion-number-parity symmetry in the mapped bosonic 
model) 48 and the other four are fermionic SPT phases 
given in Table II (three of them are non-trivial fermionic 
SPT phases and the fourth one is the trivial phase). 
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For 3D free fermion systems with T 2 = 1 time-reversal 
symmetry, there is no non-trivial gapped phase. 55,56 In 
contrast, for 3D interacting fermion systems with T 2 = 1 
time-reversal symmetry, there are (at least) three non- 
trivial fermionic SPT phases. One of the non-trivial 
3D fermionic SPT phases is actually a bosonic SPT 
phase formed by tightly bounded fermion pairs. The 
other two fermionic SPT phases cannot be obtained from 
bosonic models and are totally new. So the non-trivial 3D 
fermionic SPT phases protected by T 2 = 1 time-reversal 
symmetry cannot be constructed from free fermion sys- 
tems. 

We also find that fermion systems with symmetry 
group Gb — Z 2 (or the full symmetry group G/ = 
Z 2 x Z% ) can have one non-trivial fermionic SPT phases in 
d sp = 1 spatial dimension and three non-trivial fermionic 
SPT phases in d sp = 2 spatial dimensions. (One of 
the non-trivial 3D fermionic SPT phases is actually a 
bosonic SPT phase while the other two are intrinsically 
fermionic.) 

To compare the above result with the known free 
fermion result, let us consider free fermion systems with 
symmetry group Gb — Z2 (or the full symmetry group 
Gf = Z% X Z 2 ), which contain two kinds of fermions: one 
carries the ^-charge and the other carries ^-charge 1. 
For such free fermion systems, their gapped phases are 
classified by 44 

d sp : 1 2 3 4 5 6 7 

gapped phases : Z| Z 2 . Z 2 Z 2 ; 

The four d sp = phases correspond to the ground state 
with even or odd Z 2 -charge-0 fermions and even or odd 
i^-charge-l fermions. The four d sp = 1 phases corre- 
spond to the phases where the Z 2 -charge-0 fermions are 
in the trivial or non-trivial phases of Majorana chain 57 
and the ^-charge- 1 fermions are in the trivial or non- 
trivial phases of Majorana chain. A d sp = 2 phase la- 
beled by two integers (m,n) E Z 2 corresponds to the 
phase where the Z 2 -charge-0 fermions form m (jp x + i-Py) 
states with m right moving Majorana chiral modes and 
the Z 2 -charge-l fermions form n (jp x + ipy) states with n 
right moving Majorana chiral modes, (if m and/or n are 
negative, the fermions then form the corresponding num- 
ber of {p x — ip y ) states with the corresponding number 
of left moving Majorana chiral modes.) 

Some of the above non-interacting gapped phases have 
intrinsic fermionic topological orders. (Those phases are 
stable and have intrinsic fermionic topological orders 
even after we turn on interactions.) So only a subset 
of them are non-interacting fermionic SPT phases: 

d sp : 1 2 3 4 5 6 7 , . 

SPT phases : Z| Z 2 Z Z ^ ' 

The four d sp — phases correspond to the ground 
states with even or odd numbers of fermions and or 
1 ^-charges. The two d sp — 1 phases correspond to 
the phases where the Z 2 -charge-0 fermions and the Z 2 - 
charge-1 fermions are both in the trivial or non-trivial 



phases of Majorana chain. The d sp = 2 phase labeled by 
one integer neZ corresponds to the phase where the Z 2 - 
charge-0 fermions have n right moving Majorana chiral 
modes and the Z 2 -charge-l fermions have n left moving 
Majorana chiral modes. We sec that in 0D and ID, the 
non-interacting fermionic SPT phases are the same as the 
interacting fermionic SPT phases. 

However, in two dimensions, the non-interacting 
fermionic SPT phases are quite different from the in- 
teracting fermionic SPT phases. Despite very different 
phase diagram, it appears that the above three non- 
trivial interacting fermionic SPT phases in 2D can be 
realized by free fermion systems. In fact, it appears that 
there are seven non-trivial interacting 2D fermionic SPT 
phases protected by Z 2 symmetry 58 . All the seven non- 
trivial SPT phases can be represented by free fermion 
SPT phases. This suggests that our current construction 
is incomplete, since we only obtain four fermionic SPT 
phases (including the trivial one) with the Gf = Z 2 x Z 2 
symmetry. One possible reason of the incompleteness 
may be due to our limiting requirement that fermions 
only form ID representations of Gf. 

One way to understand why there can only be seven 
non-trivial interacting 2D fermionic SPT phases pro- 
tected by Z 2 symmetry is using the idea of duality be- 
tween intrinsic topological orders and SPT orders discov- 
ered recently 59 . The key observation in such a duality 
map is that for any SPT orders associated with a (dis- 
crete) global symmetry G, we can always promote the 
global symmetry to a local (gauge) symmetry. For differ- 
ent SPT phases, the corresponding promoted (discrete) 
gauge models describe different intrinsic topological or- 
ders. For fermionic SPT phases protected by Z 2 sym- 
metry, we can let Z 2 -charge-l fermion couples to a Z 2 
gauge field. In this way, different intrinsic topological or- 
dered phases can be characterized by different statistics 
of the Z 2 -flux. According to Kitaev's classification 60 for 
different types of vortices in superconductors 69 , we know 
that there are seven non-trivial cases. Thus, we see that 
although in free 2D fermionic systems, SPT phases pro- 
tected by Z 2 symmetry is classified by an integer (Chern 
number), the interactions dramatically changes the clas- 
sifications. 

Following Kitaev's idea 60 , we can have a very simple 
way to understand the seven non-trivial types of vertices 
by counting the number of Majorana modes in the ver- 
tex core. In the corresponding free fermion model, the 
number of Majoran modes n corresponds to the Chern 
number of the Z 2 -charge-l free fermion. The seven non- 
trivial SPT phases are described by n = 1, 2, • • ■ ,7. 

Now, let us try to understand why four of the non- 
trivial fermionic SPT phases protected by the Z 2 sym- 
metry require fermions to form high dimensional rep- 
resentations of Gf = Z 2 x Z 2 . It is easy to see that 
when n = even, the free fermion models are described 
by n/2 Z 2 -charge-l complex fermions per site and these 
fermions form ID representation of G f = Z 2 x Z 2 . How- 
ever, the situations are more complicated when n = odd. 
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FIG. 2: (Color online) Two branched simplices with opposite 
orientations, (a) A branched simplex with positive orientation 
and (b) a branched simplex with negative orientation. 



For example, when n = 1, the free fermion Hamiltonian 
describes one ^2-charge-0 Majorana fermion and one Z 2 - 
charge-1 Majorana fermion per site, labeled as 7^ and 
Under the Z 2 action, 7^ does not change while 
7»4 changes to — ^4. Thus the symmetry operation 
Uz 2 can be constructed as Uz 2 — iC^ -1 )^/ 2 Y\f =1 744, 
where N is the number of sites. It is easy to check that 

u z 2 li,t u z 2 = 7*,t> u z 2 Ji,i u z 2 = -7*4 and u z 2 = !• 
However, the symmetry operator Uz 2 can be regarded as 
the Z 2 symmetry transformation only when TV = even, 
in that case Uz 2 contains an even number of fermion op- 
erators. But in our construction, the Z 2 symmetry can 
be realized as bosonic unitary transformation regardless 
N = even or N = odd. Thus, we can only construct four 
fermionic SPT phases that are labeled by even n. 

Another example requiring fermions to carry high di- 
mensional representations is the well known time reversal 
symmetry with T 2 — Pf. We believe that the princi- 
ple/framework developed in this paper should be appli- 
cable for all these interesting cases and the results will 
be presented in our future work. 



III. GRASSMANN TENSOR NETWORK 

After having some physical understanding of the 
fermionic SPT phases protected by Z 2 x Z 2 or Z 2 x Z 2 
symmetries, in the following, we will discuss a generic 
construction that allow us to obtain fermionic SPT 
phases protected by more general symmetries in any di- 
mensions. 

One way to obtain a systematic construction of bosonic 
SPT phases is through a systematic construction of topo- 
logical bosonic path integral 50 using a tensor network 
representation of the path integral. 33 ' 61 Here we will use 
a similar approach to obtain a systematic construction of 
fermionic SPT phases through topological fermionic path 
integral using a Grassmann tensor network 25 ' 62 represen- 
tation of the fermionic path integral. The Grassmann 
tensor network and the associated fermion path integral 
are defined on a discretized space-time. So first, let us 
describe the structure of a discretized space-time. 



A. Discretized space-time 

A e?-dimensional discretized space-time is a d- 
dimensional complex E^ formed by many (/-dimensional 
simplices Eu _ i(J i. We will use ... to label the ver- 
tices of complex E^. S[ io id j is the simplex with vertices 
io, i d - 

However, in order to define Grassmann tensor net- 
work (and the associated fermion path integral) on the 
discretized space-time E^, the space-time complex must 
have a so called branching structure. 50,63 A branching 
structure is a choice of orientation of each edge in the 
(i-dimensional complex so that there is no oriented loop 
on any triangle (see Fig. 2). 

The branching structure induces a local order of the 
vertices on each d-simplex. The first vertex of a rf-simplex 
is the vertex with d outgoing edges, and the second vertex 
is the vertex with d — 1 outgoing edges and 1 incoming 
edge, etc . So the simplex in Fig. 2a has the following 
vertex ordering: 0, 1, 2, 3. 

The branching structure also gives the simplex (and 
its sub simplices) an orientation. Fig. 2 illustrates two 
3-simplices with opposite orientations. The red arrows 
indicate the orientations of the 2-simplices which are the 
subsimplices of the 3-simplices. The black arrows on the 
edges indicate the orientations of the 1-simplices. 



B. Physical variables on discretized space-time 

To define a path integral on a discretized space-time 
[ie on a complex E^ with a branching structure), we as- 
sociate each vertex i (a 0-simplex) in E^ with a variable 
gi. So gi is a local physical dynamical variable of our 
system. The allow values of gi form a space which is de- 
noted as Gb- At the moment, we treat Gb as an arbitrary 
space. Later, we will assume that Gb is the space of a 
group. 

We also associate each (d— l)-simplex (i^..ii..id) of a 

(i-simplex [io---*d] with a Grassmann number. Here the 

(d— l)-simplex (io..ii..id) is a subsimplex of the d-simplex 

[if). ..id] that does not contain the vertex ij. Relative 

to the orientation of the d-simplex [«o---V]> the (d — 1)- 

simplex (io..ii..id) can have a "+" or "— " orientation (see 

Fig. 2). We associate (d— l)-simplex (io--ii--i<i-i) with a 

Grassmann number 9, io 4 > if the (d— l)-simplex has a 

"+" orientation and with a Grassmann number 9 r ~ . 

(i ..ii..i d ) 

if the (d — l)-simplex has a "— " orientation. However, 
the Grassmann number 6u„.k) (or may or may not 

present on the (d— l)-simplex (i...k). So on each (d — 1)- 
simplex (i...k) we also have a local physical dynamical 
variable rii...k = 0, 1 or fii...k = 0,1, indicating whether 
9(i...k) or 9(i...k) is present or not. So each (d— l)-simplex 
(i...k) is really associated with a dynamical Grassmann 
variable O^ - .K or 6>?.*"A. 

We see that each (d— l)-simplex (io---i<i-i) on the sur- 
face of the cZ-complex is associated with two types of 
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FIG. 3: The 2-simplex (012) is associated with #(012) which 
belong to the top 3-simplex [0123]. The 2-simplex (012) is 
also associated with 6(012) which belong to the bottom 3- 
simplex [0124]. A rank 4-tensor V3 on the top tetrahedron 
is given by Vg(go, 9i, 92, 53) and on the lower tetrahedron by 

V3" (go, 91, 52,54) ■ 

Grassmann variables 8,*°'"* d ~\ and 6,*°'"* d ~\, each be- 
longs to the ci-simplex on each side of the (d — l)-simplex 
(io...id-x) (see Fig. 3). We also see that each d-simplex 
has its own set of Grassmann variables that do not over- 
lap with the Grassmann variables of other simplices. 

C. A constraint on physical variables 

In general, the physical dynamical variables gi, rii...fc, 
and ni...k are independent. We can formulate fcrmionic 
topological theory and obtain fermionic topological 
phases by treating gi, rii...u, and ni...u as independent 
variables. This is essentially the approach used in Ref. 25 
in (2+l)D, and we have obtained a system of non-linear 
algebraic equations, whose solutions describe various 2D 
fermionic topological phases. However, the non-linear al- 
gebraic equations are very hard to solve, despite they 
only describe 2D fcrmionic topological orders. As a re- 
sult, we have obtained very few fermionic solutions. 

In this paper, we are going to study a simpler and 
less general problem, by putting a constraint between 
the physical dynamical variables gi and ni...u (fii... k ): 

rii...k = n(g l ,...,g k ), fii.„ k = n(g i} g k ). (3) 

That is when m...k 7^ n (9t, ■■■,9k) or fk...k n (9t, -,9k), 
there will be a huge energy cost. So at low energies, 
the system always stays in the subspace that satisfy the 
above constraint. In this case rii...k, and fii... k are deter- 
mined from <7j's. 

It turns out that the constraint indeed simplifies the 
mathematics a lot, which allows us to systematically de- 
scribe fermionic topological phases in any dimensions. 
The resulting non-linear algebraic equations also have 
nice structures which can be solved more easily. So we 
can also find many solutions in any dimensions and hence 
obtain many new examples of topological phases in any 
dimensions. However, the constraint limits us to describe 
only SPT phases. To describe phases with intrinsic topo- 
logical orders, we have to use more general unconstraint 



formulism, which is discussed in Ref. 25. In this paper, 
we will only consider the simple constraint cases. 

The above constraint can be generalized a little bit. 
In general, we may have several kinds of fermionic, or 
even bosonic, variables in the (d — l)-simplices. So we 
can associate (09. .) ni o-U-i r (69. ) ni o-*d-i 
with each (d — l)-simplex. Here nf o t =0,1 indi- 
cating the presence or the absence of the a th fermionic 
variable in the (d — l)-simplex (iQ...id-i)- We can 
also associate bosonic variables (<bf. ) i o- -*d-i or 

(4>? ia id t i o--- i <«-i with each (d— l)-simplex, where 4> a 
and (f> a are complex phases. Here Z° ^ = 0,1,2,... 
indicating the number of the a th bosonic variable in the 
(d — l)-simplex (iQ...id-i)- Also, a, a — 1,2,... label dif- 
ferent kinds of fermionic or bosonic variables, which may 
be in different ID representations of the symmetry. 

In this more general situation, we impose the following 
constraint between the physical dynamical variables gi 
and nf k (n?. k ): 

„«(«.-,») = n( 9i ,...,g k ), <;l n(g g k ), 

other n9, k ,fi9 k = 0, 

&{99h—,99k) =a(j9i,-,9k)- ( 4 ) 

So for each fixed set gi , . . . , g k , only one kind of fermionic 
variable with a = a(gi, g k ) can appear. We also im- 
pose the following constraint between the physical dy- 
namical variables gi and i° k (lf k ): 

ll. k =l a (9i,-,9k), ll. k = l a ( gi ,...,g k ), 
a(99i,-;99k) =a(g i ,... ) g k ). (5) 

So for each fixed set gi, ...,g k , the number of the bosonic 
fields for each kind cannot fluctuate. They are given by 
fixed values l a (gi, ...,g k ) and l a (gi, ...,g k ). This weaker 
constraint is enough to simplify the mathematics, since 
after the constraint, only gi on vertices can fluctuate in- 
dependently at low energies. The degrees of freedom on 
the (d— l)-simplices are tied to g^s and cannot fluctuate 
independently. 

We note that different degrees of freedom are attached 
to the (d — l)-simplices for different choices of gi, ...,g k - 
Those degrees of freedom may carry different ID rep- 
resentations of the symmetry. Later, we will use this 
picture to construct path integrals that respect the sym- 
metry. 

D. Grassmann tensor on a single d-dimensional 
simplex 

In the low energy constraint subspace, we can associate 
each d-simplex of "+" orientation with a Grassmann rank 
(1 + <i)-tensor Vj~, and each d-simplex of "— " orientation 
with a different Grassmann rank (1 + (i)-tensor V^. 

Let us first assume that a d-simplex [0...d] has a "+" 
orientation, and discuss the structure of the Grassmann 



tensor Vj~. The Grassmann tensor Vj~ is a map from 
space G\ +d -> My: 



Vj~(9o,9i,-,9d) e M/, g t eG b . 



(6) 



The order of the variables go, 9i, 9d in Vd is the same as 
the order of the branching structure: vertex-0 < vertex- 1 
< ... < vertex-d. 

Here an elements in Mf is a complex number 
times an even number of Grassmann numbers. Thus 
Vj" (90)51, -,9d) has a form 



Vj~ (#o,5i, -,5d) = "a (go, 9u92,9a,94, -,9d)x 

a«d-i(9i,92,93,94,-..,9d)/i™d-i (9o,9i,93,94,-..,9d) 
"(1234.. .d) "(0134.. .d) — X 

n n d-l{go, g2,93, 94, 9d) n n d-l{ga, 9l,92, g4---,gd) 
"(0234. ..d) "(0124. ..d) 



(7) 



where fj" (90,91, ■■■,g d ) is a pure phase factor and Ouj„.k) 
or Ouj...u) is the Grassmann number associated with (d — 
l)-dimensional subsimplcx on E ,, (ij---k). The constraint 



q n 1234...d 



(1234. ..d) 



as 



is implemented explicitly by writing, say. 

/V»d-l(ffl>92, 93,94,, ■ ■■,9d) 
"(1234. ..d) 

Note that the vertices on the branched simplex E , have 
a natural order: vertex-0 < vertex- 1 < ... < vertex- 
d. This leads to the ordering of the Grassmann num- 
bers: the first Grassmann number is associated with the 
subsimplex that does not contain the vertex-0, the sec- 
ond Grassmann number is associated with the subsim- 
plex that does not contain the vertex-2, etc . Then it is 
followed by the Grassmann number associated with the 
subsimplex that does not contain the vertex- 1, the Grass- 
mann number associated with the subsimplex that does 
not contain the vertex-3, etc . 

The above example is for simplices with a "+" orien- 
tation. For branched simplex E^ with a "— " orientation, 
we have 



V d (go,9i,-,9d) = v d (go,9i,g2,93,94,-,9d)x 
n n d-i(go,gi,g2,g4---,gd) a n d-i(go,g2,g3,94,---,gd) 

■■•"(0124. ..d) "(0234. ..d) X 

n n d-i(go,9i,g3,g4,---,gd) n n d-i(gi,g2,g3,g4,---,gd) 

"•"(0134.. .d) "(1234. ..d) 



(8) 



where the order of the Grassmann numbers is reversed 
and 9's are switched with (9's. 

The number of the Grassmann numbers on the (d— 1)- 
dimcnsional simplex (ij...k) is given by 



n>d-i(9i,9j,—,gk) = 0, 1. 



(9) 



which depends on d variables gi, gj, 9k, and must sat- 
isfy 



d 

E 

i=0 



nd-i(go, — ,9i,—,9d) = even, 



(10) 



so that the Grassmann rank (1 + d)-tensor always has an 
even number of Grassmann numbers. Here the sequence 
go, ...,g"i, ...,gd is the sequence <?o,---,<?d with gi removed. 



All the solution of eqn. (10) can be obtained using an 
integer function m c i~2{gi, 92, ■ g d -i) =0,1 with d — 1 
variables (while fid-i(9i, 92, ■■■,gd) is an integer function 
of d variables): 

nd-i{gi,g2, -,9d) 

d 

= ^2 m d-2{gi, gi, g d ) mod 2. (n) 



For example, a rank 1-tensor has a form 

Vo(ffo) = f (<?o) (12) 

which contains no Grassmann number. A rank 2-tensor 
has a form 

Vt(9o, gi) = 4(9o,gi)e^ 9l) ~e n { °l 9o) (13) 

where no (9) has only two consistent choices 

no(g) = 0, Vg G G b , and n (g) = 1, Vg G G b 



(14) 



As another example, a rank 4-tensor V% for the tetrahe- 
dron in Fig. 2a is given by 



^(90,51,92,93) = v£(j9o, 9i, 92, 93) x 

/i«2 (gi ,92,93) n«2 {90,91,93) a n 2 (go ,92 ,33 ) 7p2 (go ,gi ,92 ) 

"(123) "(013) "(023) "(012) 



(15) 



Note that the four triangles of the tetrahedron have dif- 
ferent orientations: the triangles (123) and (013) point 
outwards and the triangles (023) and (012) point inwards. 
We have used 6 and 6 for triangles with different orien- 
tations. 

For the tetrahedron in Fig. 2b, which has an opposite 
orientation to that of the tetrahedron in Fig. 2a, we have 
a rank 4-tensor VJ" which is given by 



V 3 (90,91,92,93) = v 3 (90, 9i, 92, 93) x 

fl"2 (so ,91 ,92)/)l2(ffo ,9: 

"(012) "(023) "(013) "(123) 



"2 (go ,Sl ,92 ) a n 2 (go ,92 ,93 ) a n 2 (go ,91 ,93) g"2 (91 ,92 ,93 ) ^g-j 



Note the different order of the Grassmann numbers 

o n 2 (go ,91 ,92 ) a n 2 (go ,92 ,93 ) rjn2 (90 ,91 ,93 ) /)«2 (91 ,92, gs) 
"(012) "(023) "(013) "(123) 



E. Evaluation of the Grassmann tensors on a 
complex 

Let Yid be a d-complex with a branching structure 
which is formed by several simplices [ai>...c]. As dis- 
cussed in the last section, a simplex is associated with 
one of the two Grassmann tensors Vjr. Let us use 



s(a,fc,...,c) 



to represent the evaluation of 



Jin(Sa) ri[o6...c] 

rank (1 + d) Grassmann tensors V d on the d-complex 
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I II y; (oA -' o) 

v a > [ab...c\ 



(17) 



/n 



cl(9 



nd-l(gi,gj,---,gk) < ^0n d - 1 (g i ,g j ,...,g k ) 



(ij...k) 



(ij—k) 



(ij...k) 
{xy...z} 



where Y[[ab c] multiplies over all the d-dimensional sim- 

plices [ab...c] in E^ and V d b '"' (g a , gb, ■■-,9c) is a Grass- 
mann rank (1 + (f)-tensor associated with the simplex 
[ab...c]. s(a,b, ...,c) = ± depending on the orienta- 
tion of the d-dimensional simplex [a6...c]. Also Yiuj k) 
multiplies over all the interior (d — l)-dimensional sim- 
plices, (ij...k), of the complex E^ [ie those (d — 1)- 
dimensional simplices that are not on the surface of 
Erf]. Yl{ xy z j multiplies over all the interior (d — 2)- 
dimensional simplices [ie those (d — 2)-dimensional sim- 
plices that are not on the surface of EJ, {xy...z}, of 
the complex E^. We note that, when E^ is a sin- 



gle simplex [ab...c], / in(Ed) II 



V 



s(a,b,...,c) 



V. 



s(a,6,. 



is given by 



d 

We 



(9 a, 9b, -,9c 
see that 



the formal notation, 



L(z d )T\[ab...c] V d > Sin(E d ) ^presents a Grass- 

mann integral over the Grassmann numbers on all the 
interior (d— l)-simplicies in the d-complex E^. So in(E,j) 



actually represents the collection of all those 



interior {d — l)-simplicies in the d-complex E^. The 
Grassmann numbers on the (d — l)-simplicies on the 
surface of the d-complex E^ are not integrated over. In 
fact / in ( Sd ) has the following explicit form 



nd-l(gi,gj,---,gk) ^Qn d - 1 (g i ,gj,...,g k ) ^ 



, , , II <%••..*) 
^ J {ij...k) 

TJ (_yna-2{g*,gy,---,gz) 

{xy...z} 



(ij...k) 



(18) 



We note that d9 always appears in front of d9. We also 
note that the integration measure contains a non-trivial 
sign factor (— ) m <*-2(s*,9!/,— ,s*) on the interior (d — 2) sim- 
plices {xy...z}. 

The sign factor ]J {xy (-)™*-2(9x,g y ,-,g*) \ s included 
to help us to define topological fermionic path integral 
later. Choosing orientation dependent tensors Vf also 
help us to define topological fermionic path integral. 
Adding the sign factor H. {xy ^y {-) m "-^a^s v , -,9*) an d 

choosing orientation dependent tensors V d appear to be 
very unnatural. However, they are two extremely impor- 
tant features of our approach. We cannot obtain topolog- 
ical fermionic path integral without these two features. 
Realizing these two features is one of a few breakthroughs 
that allows us to obtain topological fermionic path inte- 
gral on discrete space-time. The two features are related 



to our choice of the ordering convention of 9 and 9 in 
the Grassmann tensors V d , and the ordering convention 
of d9 and d9 in the integration measure. We have not 
find an ordering convention to remove the sign factor 
n, , c„ N |"id-2(gx,ff !J ,..-,9 2 ) 

ll{xy...z}\ / 

Now let us consider two d-complices E^ and E^, which 
do not overlap but may share part of their surfaces. Let 
Ed be the union of the two d-complices. From our defini- 
tion of the Grassmann integral, we find that the Grass- 
mann integral on E^ can be expressed as 



I II ltf a ' W) 

Jin (^d) [ab...c] 



(19) 



(a,b,. 



n ° A - ,c) / n v * 

d d \ d' [ab...c] v dJ [ab...c\ 



where E^ n T? d is the intersection of the two complices, 
which contains only (d — l)-simplices on the shared sur- 
face of the two complices E^ and E^. More precisely 



II d %,..k) 

(ij...k) 

TJ [_^rn,d-2(gx,gy,..;g*) 
{xy...z} 



(ij...k) 



(20) 



where Yiuj k) * s a P r °duct over all the (d — l)-simplices 
in E^ n E^ and Yi{xy z} ^ s a product over all the interior 
(d-2)-simplices in E d nE^ (note that E d nE^ is a (d-1)- 
complcx). Eqn. (19) describes the process to glue the two 
complices E d and Y? d together. 



F. An example for evaluation of Grassmann tensor 
network 

Let us consider an example to help us to understand 
the complicated expression (17). When E is formed by 
two tetrahedrons as in Fig. 3, the top tetrahedron with 



a "+" orientation is associated with a rank 4-tensor V. 



3 ' 



(90, 9i, 92, 93) = v£ (go, 91, 92,93) ■ 



q«2 (gi ,32 ,93 ) a"2 (go ,91 ,93 ) n n 2 (go ,92 ,93) /)«2 (go ,91 ,92 ) 



'(123) 



'(013) 



'(023) 



'(012) 



(21) 



The lower tetrahedron with a "— " orientation is associ- 
ated with a rank 4-tensor V7 : 



V3 (9o,9i, 92, 9i) = v z (,g , 51,52, 94) x 

n n 2 {go ,91 ,92 ) fl ™2 (90 ,92 ,94) H n 2 (go ,9i ,94 ) n n 2 (91 ,92 ,94 ) 



'(012) 



'(024) 



'(014) 



'(124) 



(22) 



Note that V3" (50 , 51 , 92 , 53 ) and V 3 (50 , 51 , 52 L 54 ) have in- 
dependent Grassmann numbers #(012) and #(012) on the 
shared face (012). #(012) belongs to V^(go, 5i, 52, 54) and 
0(oi2) belongs to V£(g ,gx,g 2 ,g3)- Since in(E) = (012), 
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the evaluation of V* on the complex S is given by 
/ V 3 + (g , 9i, 52, 9z)Vz (go, 9i,9i, 9a) 

J(012) 

/jQH2(9o,gi,g2) j/7«2 (go, 31,32) v 
Qt, (012) a °(012) X 

Vs" (ffo, 3i, 92, 93W3 {g , 91, 92, 54) (23) 
which is a special case of eqn. (17). 

IV. GRASSMANN TENSOR NETWORK AND 
FERMIONIC PATH INTEGRAL 

Given a closed space-time manifold, M$t, with a time 
direction, its triangularization £ is a complex with no 
boundary. The time direction gives rise to a local order of 
the vertices. So the complex E has a branching structure. 

Since £ has no boundary, the evaluation of Grassmann 
tensors V d on it gives rise to a complex number: 




Such a complex number can be viewed as the action- 
amplitude in the imaginary-time path integral of a 
fermionic system. The different choices of the com- 
plex functions ^(50, 5d) an d the integer functions 
m d-2(go, ■■■,9d-2) correspond to different choices of La- 
grangian of the fermionic system. The partition function 
of the imaginary-time path integral is given by 

Z = H[ II V d • (25) 

{g,} Jin ^ [ab...c] 

where each vertex of the complex £ is associated with a 
variable gi. We note that 



the fermionic partition function is determined 
from two (1 + G?)-variable complex functions: 

v d(9i,9j,-), v didi,9j,-), and one (d - In- 
variable integer function md-2(gi, gj, ■■■)■ 



V. TOPOLOGICALLY INVARIANT 
GRASSMANN TENSOR NETWORK AND 
FERMIONIC TOPOLOGICAL NON-LINEAR 
a-MODEL 

Now we would like to study the low energy effective 
theory of a gapped fermion system. The fixed point 
low energy effective action-amplitude of a gapped system 
must describe a topological quantum field theory. So the 
fixed point action-amplitude must be invariant under a 
renormalization group (RG) flow which generates coarse- 
graining transformation of the space-time complex £. In 
fact, such a fixed-point requirement under RG flow com- 
pletely fixes the form of the fixed-point action-amplitude. 



t t 




#3 Si §3 8l 



(b) 



FIG. 4: (a) Consider two triangles which are a part of a 
space-time complex E with a branching structure. The first 
type of the RG flow step changes the two triangles to other 
two triangles. The branching structure leads to the following 
vertex ordering (0, 1, 2, 3). (b) A different space-time complex 
E which corresponds to a different ordered vertices (0, 3, 1, 2). 
The corresponding complex does have a valid branching struc- 
ture because some triangles overlap. 

The fixed-point action-amplitude, in turn, classifies the 
possible fermionic SPT phases. 

The detailed RG flow steps depend on the dimensions 
of the space-time. We will discuss different space-time 
dimensions and the corresponding fixed-point action- 
amplitude separately. Here we will only give a brief dis- 
cussion. A more detailed discussion will be given in Ap- 
pendix A. 

A. (1+1)D case 

We first consider fermion systems in (1 + 1) space-time 
dimension. Its imaginary-time path integral is given by 

z =e/ n v 2 (Q - b ' c \ (26) 

{ gi } J ' m ^ [abc] 

where £2 is the (1+1)D space-time complex with a 
branching structure (see Fig. 4), and are rank-3 
Grassmann tensors. 

The first type of the RG flow step that changes the 
space-time complex is described by Fig. 4a. If the Grass- 
mann tensors describe a fixed point action-amplitude, 
their evaluation on the two complices in Fig. 4a should 
be the same. This leads to the following condition 

f 2 (.90, 5i ,53)^2(31, 32, #3) 
= ^2(90, 9i, 92) v% (9o,92, m) (27) 

where ^2 is a complex phase related to and nio(g) 
through 

^2(90,91,92) = ^2(30,51,52), 

^(5o,5i,52) = R mo(9l) M(5o,5i,52). (28) 
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Si (a) 8\ 
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FIG. 5: (a) Consider three triangles which are a part of a 
space-time complex £ with a branching structure. The sec- 
ond type of the RG flow step changes the three triangles to 
one triangle. The branching structure leads to the following 
vertex ordering (0,1,3,2). (b) The vertex order (0,1,2,3) 
does not correspond to a valid branching structure. 




FIG. 6: (a) A 3D complex formed by two 3-simplices 
(do, 51,52,53) and (51,52, 53,54)- The two 3-simplices share 
a 2-simplex (51,52,53). (b) A 3D complex formed by three 
3-simplices (50,52,53,54), (50,51,53,54), and (50,51,52,54). 

Note that the above condition is obtained for a par- 
ticular vertex ordering (0,1,2,3) (a particular branch- 
ing structure) as described in Fig. 4a. Different valid 
branching structures, in general, lead to other conditions 
on ^2(50, 5i, 52)- It turns out that all those conditions are 
equivalent to eqn. (27). For example, for vertex ordering 
(1, 3, 0, 2), we obtain 

^2(51,50, 52)^2(51,53,30) 
= ^2(51, 53, 52W53,5o,52) (29) 

which is the same as eqn. (27) after replacing g\ — > g , 
53 -> 5i, 5o -> 52, and g 2 -> 53- 

Some vertex orders do not correspond to valid branch- 
ing structure due to the overlap of the simplices, and are 
not considered. One of the invalid ordering is described 
by Fig. 4b which has a vertex order (0, 3, 1, 2). The in- 
valid orderings can give rise to conditions on ^2(50, 5i, 52) 
that are not equivalent to eqn. (27). 

The second type of the RG flow step that changes the 
space-time complex is described by Fig. 5. Different 
valid branching structures (see Fig. 5a for one example) , 
lead to different conditions on ^2(50,51,52) (which is de- 
fined in eqn. (28)). It turns out that all those conditions 
are equivalent to eqn. (27). (Some vertex orderings do 




FIG. 7: (a) A 3D complex formed by one 3-simplex 
(50,51,52,54). (b) A 3D complex formed by four 3- 
simplices (51,52,53,54), (50,52,53,54), (50,51,53,54), and 
(5o,5i,52,53)- 




FIG. 8: (Color online) (a) A 4D complex formed by two 
4-simplices (50,51,52,53,54) and (gi, g 2 , 53, 54, 5s)- The two 
4-simplices share a 3-simplex (51,52,53,54) (color red), (b) 
A 4D complex formed by four 4-simplices (50,52,53,54,55), 
(5o,5i,53,54,55), (5o,5i,52,54,55), and (5o,5i,52,53,55)- 
The two simplices (50, 51, 52, 54, 5s), and (g , gi, g 2 , 53, 5s) 
share a 3-simplex (50,51,52,55) (color blue). 

not correspond to valid branching structures. One of the 
invalid ordering is described by Fig. 5b.) 

Since ^2(50, 5i, 52) is related to ^(50,51,52) and 
mo (5) through eqn. (28), eqn. (27) is actually a condition 
on V2 (go , 5i , 52 ) and mo (5) that determines the fermion 
path integral. The fermion path integral described by 
v 2 (go j 51,52) and mo (5) that satisfies eqn. (27) is a topo- 
logical fermion path integral which is a fixed-point theory 

We would like to point out that eqn. (27) is the stan- 
dard 2-cocycle condition of group cohomology. 50 ' 51 



B. (2+l)D case 

We next consider fermion systems in (2 + 1) space-time 
dimension described by imaginary-time path integral 

Z = Y,[ II V 3 s(aAc ' d) , (30) 
{g,} ' /in(E3) [abed] 

where £3 is the (2+l)D space-time complex with a 
branching structure and are rank-4 Grassmann ten- 
sors. 

The first type of the RG flow step that changes the 
space-time complex is described by Fig. 6. In order 
for the Grassmann tensors V3 to describe a fixed point 
action-amplitude, their evaluations on the two complices 
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JU &4 branching structures give rise to the same condition (31). 

X/-J \v /vi 1 n. The second type of the RG flow step that changes the 

/ flS^T ^-^Vg space-time complex is described by Fig. 7. If the Grass- 

^ _t_/V - - — ~Z^^ 3 — Z^- ^' ^-A*^ 5 ~ " " -J^^ 3 mann tensors describe a fixed point action-amplitude, 

So^^,/^^^'^ 8o^^},[^^^^ their evaluations on the two complices in Fig. 7 should 

g 2 82 also be the same. This leads to the same condition (31). 

(a) (b) 



FIG. 9: (Color online) (a) A 4D complex formed 

by one 4-simplex (50,51,52,53,54). (b) A 4D complex 
formed by five 4-simplices (51 , 52 , 53 , 5s , 9a ) , (go , 52 , 53 , 5s , 54 ) , 
(30,51,53,35,54), (50,51,52,55,54), and (50,51,52,53,55)- 
The two simplices (50,51,52,55,54), and (50,51,52,53,55) 
share a 3-simplex (50,51,52,55) (color blue). 

in Fig. 6 should be the same. This leads to the following 
condition 

^3(31, 32, .93, 94)^3(90, 5i,33,34)^3(3o,5i, 92, 93) = (31) 

where 1^3 is a complex phase and is given by 

4(30,31,32,33) = (-r i(9o ' ff2 V 3 (3o,3i,32,33), 
^(30,31,32,33) = R mi(9l ' ff3) M(3o,3i,32,33) (32) 

Different vertex orders give rise to different branch- 
ing structures. It turns out that all the different valid 



C. (3+l)D case 

Last, we consider fermion systems in (3 + 1) space-time 
dimension, described by imaginary-time path integral 

Z=J2 [ II Vl (aMe) , (33) 

{Si} (S4) [abede] 

where £4 is the (3+l)D space-time complex with a 
branching structure and V± are rank-5 Grassmann ten- 
sors. 

The first type of the RG flow step that changes the 
space-time complex is described by Fig. 8. The second 
type of the RG flow step that changes the space-time 
complex is described by Fig. 9. If the Grassmann tensor 
V4 describes a fixed point action-amplitude, its evalua- 
tion on the two complices in Fig. 8 or Fig. 9 should be 
the same. This leads to the following condition 



^4(31,32,33,34,35)^4(30,31,33,34,35)^4(30,31,32,33,35) 

r_^n3(90,9i,S2,93)n3(ffo,93,94,55)+n3(si,ff2,93,ff4)n3(so,Si,S4,ff5)+n3 (92,33,94, Ss)«3 (9o,9l,ff2,fl5) . 
^4(30,32,33,34,35)^4(30,31,32,34,35)^4(30,31,32,33,34) 



(34) 



where is a complex phase and is given by 
4(30,31,32,33,34)= (35) 

(-)^(so^.S3)+m 2 (9 I ,93,94)+m 2 (9i,92,93) z , 4 ( 50)ffl)52i53)54 ) ) 

^(30,31,32,33,34) = (-r 2(90 ' 92 ' 94) M(.90,3i,32,33,34) 

The sign factors in above relations between v± and ^4, 
and the sign factor in eqn. (34) are very strange. Obtain- 
ing those strange sign factors is another break through 
that allows us to obtain topological fermionic path in- 
tegral beyond (1+1) dimensions. It appears that those 
sign factors are related to a deep mathematical structure 
called Steenrod squares. 64 ' 65 Different vertex orders give 
rise to different branching structures. It turns out that 
all the valid branching structures give rise to the same 
condition (34). 



D. The fixed-point action on a closed complex and 
the unitary condition 

We know that the fixed-point theory of a bosonic 
SPT phase is a discrete bosonic topological non-linear 
cr-model. 50 The bosonic topological non-linear cr-models 
are characterized by action-amplitudes that is equal to 
1 on any closed complex. The action-amplitude being 
equal to 1 on any closed complex is the key reason why 
the models describe bosonic SPT phases without intrinsic 
topological orders. 

In Section V, we have developed a fixed-point theory 
of fermionic SPT phases, and defined discrete fermionic 
topological non-linear cr-model. A discrete fermionic 
topological non-linear cr-model is described by a pair of 
functions [v d (g , ...,ga), TOd-2(3o> —,9d-2)], tnat satisfy 
the conditions eqn. (27), eqn. (31), or eqn. (34). Accord- 
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ing to the branching moves discussed above, the evalua- 
tion on any closed manifold is reduced to an equivalent 
evaluation 

/ VjVJ (36) 

where the symbol J means integrating out all the 
Grassmann variables on the corresponding d — 1 sim- 
plices (see eqn. (18)). We note that the sign factor 
U{ xy z} (-) md -^" g v-' 9 ^ should be included in the 
fermionic path integral. In Appendix B we will show 
that J V d V d IS always equal to 1. 

We also note that the bosonic topological non-linear 
cr-models are characterized by act ion- amplitudes that 
are pure U(l) phases \v d {go, ■•■,9d)\ — 1 (where the 
action-amplitude on "+" oriented simplices is given by 
v d (go, ...,g d ), and the action-amplitude on "— " oriented 
simplices is given by v d l (go, —,9d) = ^(sfo> -,9d))- The 
action-amplitude being a pure U(l) phase ensures that 
the model defined by the path integral to be unitary the- 
ory. 

But what is the analogue of the pure U(l) phase con- 
dition on the Grassmann amplitude (or on the pair of 
functions [u d (g , g d ) 1 m d ^ 2 (go, ...,g d -2)W Clearly, the 
Grassmann amplitude V d are not even complex numbers. 
It is hard to say when V d behave like U(l) phases. 

To address this issue, let us introduce complex conju- 
gate of a quantity that contains Grassmann numbers: 

(W 2 03-)* = u*...h6 2 h- (37) 

We see that under complex conjugate a) the complex 
coefficients are complex conjugated, b) the order of the 
Grassmann number is reversed, and c) 9 and 9 is ex- 
changed. With this definition, one is tempting to require 
Vj(Vj)* = 1, in order for V+ to be a 17(1) phase. But, 
V d {V d )* still contains Grassmann numbers and we can- 
not require it to be 1. So we try to require J Vj~(Vj~)* = 1 
where the Grassmann integral is defined in eqn. (18). 
Now JVj(Vj)* is a complex number. But / Vj~(Vj~)* 
is not non-negative. So we cannot treat fVf(Vt)* as a 
norm-square ofVj". After some considerations, we find 
that we need to define a different complex conjugate: 

[Vt(9o,9i)} f = [Vt(9o,9i)]*, 

[Vt {90, 91,92)? = (-) mo(9l) [V 2 ± (.go,ffi,52)]*, (38) 

[Vf{ 90 ,..,9z)] ] = (-) mii9o ' 92)+mi{9l ' 93) [Vi(go,..,g3)}*, 

[Viigo,-^^ = [Vt(9o,-,94)]*x 

^_\OTi(90iffi,sa)+'Tii (91 ,g3,94)+TOi (si ,g 2 ,S3)+mi (90,32,94) 

Using the new complex conjugate and the evaluation of 
/ described above, we find that 

/v d + [V+]t= |v7[V7]t=^. (39) 



Thus we would like to require the Grassmann amplitude 
V d to satisfy 

/v d ± [V d ± ]t = l. (40) 

which is to require \v d (go, 9d)\ — 1- Eqn. (40) is the 
analogue of the pure U(l) phase condition on the Grass- 
mann amplitude V d . 

For the Grassmann amplitude V d that satisfy the pure 
U(l) phase condition (40), we have 

[V^=VJ- (41) 

E. Symmetry and stability of the fixed-point action 
amplitude 

We have seen that a fermionic system can be described 
by a fermionic path integral on the discretized space- 
time. In d dimensional space-time, the fermionic path in- 
tegral is determined by two (d-l-l)-variable complex func- 
tions v d (go, g d ) and one (d — Invariable integer func- 
tion TOd-2(ff0j <?d-2)- In the last few sections, we have 
shown that if v d (go,—,gd) and 1714-2(90, ... ,fld-a) sat- 
isfy some conditions, eqn. (27), eqn. (31), and eqn. (34), 
the fermionic path integral determined from v d (go, g d ) 
and 771,^-2(505 •••) 9d~2) is actually a fixed-point theory un- 
der the coarse-graining transformation of the space-time 
complex. 

A fixed-point theory can be used to describe a phase 
if the fixed-point theory is stable. To see whether 
a fixed-point theory is stable or not, we perturb the 
fixed-point action amplitude v d (go, g d ) by a small 
amount v d (go,..., gd) -> v d (go, -,9d)+8v d (go, -,gd)- H 
Vd(9o, -i9d) + Sv d (go, -,gd) flows back to v d (go, —,9d) 
under the coarse-graining transformation of space-time 
complex, the fixed-point theory is stable. 

From this point of view, our constructed fixed- 
point theories described by v d (go, ...,gd) and 
m d-2(9o, gd-2) are not stable. If we add a per- 
turbation of form 81/4 (go,..., gd) = ( Ilt=o^g 1 ,g to 
increase the weight of g state in the action amplitude, 
we believe that the action amplitude will flow to the 

simple form v d (g , —,gd) = ]li=o <W where only the g 
state has a non-zero weight. 

However, under a symmetry condition, our constructed 
fixed-point theories can be stable. To impose a symmetry 
condition, we assume that our fermion system has full 
symmetry group Gf. Let Gb = Gf/z£, we will take gi in 
v ^d(9o,--,9d) and md- 2 (go, — ,9d-2) as a group element 
in G;,. From Section III C, we see that different degrees of 
freedom are attached to the (d— l)-simplices for different 
choices of gi, ...,<?ft. The degrees of freedom on a (d — 1)- 
simplex form a ID representation u 9 d _ 1 of the symmetry 
g S Gf, which may depend on gi, ...,gk- 

u d -i(9i:9j,-,9k) = u 9 d 1 _ 1 (gi,g j ,...,gk)uf_ 1 (g i ,g j ,...,gk) 

u d -i(g'gi,gj,.-,g'gk) =u 9 d _ l (g i ,g j ,...,gu), (42) 
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We also require the variables with bar and without bar 
to carry opposite quantum number (ie their ID repre- 
sentations are inverse of each other). In this case, the 
symmetry of the path integral (ie the invariance of the 
action amplitude) can be implemented by requiring: 

6(ij...k) -^ U g 1 (9i,9j, ■■■,9k)d(ij...k), 

%?...fe) ->Ug(j9i,9j,—,9k)Q(ij...k), (43) 

and 

i^(99o,-,99d)) s{9) 
= v d (go, ■■,gd)Y[(u d _ 1 ) { ~ r (g , -,9i, -,9d), 

i 

m d - 2 (gg ,..,gg d - 2 ) = m d - 2 (go,-,gd-2), g€G f . (44) 

where s(g) = —1 if g contains the anti unitary time- 
reversal transformation and s(g) = 1 otherwise. If we 
view V d (go 7 ••) g<i)) as a whole object, it is invariant under 
the symmetry action g: 

Vd(ggv,~,ggd)) s{9) = v^(g ,..,g d ) (45) 

Indeed, the above symmetric v d (g , g d ) and 
m d-2(go, ■••, gd-2) describe a fermion system with Gf 
symmetry. If v d (g ,...,g d ) and m d _ 2 (g , ■•■, gd-2) fur- 
ther satisfy eqn. (27), eqn. (31), or eqn. (34), we believe, 
they even describe a fixed-point theory that is stable 
against perturbations which do not break the symme- 
try. But further consideration suggests that such a sym- 
metry condition is too strong. The fermionic systems 
that are described by the symmetric v d (g , g d ) and 
"^-2(50, ■ ••, gd-2) appear to be essentially bosonic. At 
least, their classification is the same as the classification 
of bosonic SPT phases. 

So in this paper, we will use a weaker symmetry condi- 
tion. We only require v d (g , ...,g d ) and n d ^i(g , ...,g d -i) 
(which are certain combinations of is d (go, g d ) and 
m d - 2 (go, -,9d-2)) to be symmetric 

Vd 9 \ggo, -,99d) 
= ^d(.9o,--,gtf)n( w d-i) ( ~ r (3o, ->9i, -,9d), 

i 

n d -i(gg , ..,gg d -i) = n d -x(go,..,g d -\), g^G f . (46) 

Since the fermion fields always change sign under the 
fermion- number-parity transformation Pf e Gf, the ID 
representations satisfy 

u Pf (g ,..,g d - 1 ) = (-) n *-^°>->^. (47) 

Note that we only require n<j-i(5o, 9d-i) to be sym- 
metric. In general, m d _ 2 (gQ, g d - 2 ) is not symmet- 
ric m d ^ 2 (gg ,...,gg d ^ 2 ) ^ m d _ 2 (5o, g d - 2 ). Although 
the fermionic path integral appears not to be symmet- 
ric under the weaker symmetry condition, it turns out 
that the symmetry breaking is only a boundary ef- 
fect and the symmetry can be restored by adding an 



additional boundary action to the space-time path in- 
tegral. In fact, we will show that the ground state 
wavefunction and the Hamiltonian obtained from the 
path integral with the weaker symmetry condition are 
indeed symmetric. So a fermion path integral deter- 
mined by u d (g , ...,3d) and n d ^ 1 (g , ...,g d -i) that satisfy 
eqn. (46) as well as eqn. (27), eqn. (31), or eqn. (34) 
does describe a symmetric stable fixed-point theory. Such 
"d(go, -,9d), n d -!(g Q , ...,g d -i), and u 3 d _ x (go, -,g d -i) de- 
scribe a fermionic topological phase with symmetry. 
Let us introduce 

h (5o,3i) = 0; 
h(9a, 51,32) = 0; 
/ 3 (5o,5i,-,5 3 ) = 0; (48) 
h{9o,9l, ■•■,54) = n 2 (g Q , gx, g 2 )n 2 (g 2 , g 3 , g 4 ); 

h(go,gx,-,g^) = "-3(50,51,32,33)^3(50,33,34,55) 

+ "3(51, 52, 53, 54W50, 5i, 54, 55) 

+ "3 (52 , 53 , 54 , 55 ) "3 (50 , 51 , 52 , 55 ) • 

Using f d , we can define a mapping from (d + 1)- 
variable functions v d (go, ■■■,3d) to (d + 2)-variable func- 
tions (6v d )(g Q , ...,g d+1 ): 

(8vd)(go,-,gd+i) (49) 

d+1 

= H^+^°--^) [J u d - }i (g , ..,g d+1 ). 

4=0 

Then using (5v d ), we can rewrite the conditions, 
eqn. (27), eqn. (31), and eqn. (34), in a uniform way 

(8vd)(9o,-,gd+i) = 1. (50) 

Here [v d (g , g d ), n d -i(g , ...,g d _x), u^^go, 9d-i)] 
that satisfies eqn. (46), eqn. (10), and eqn. (50), will be 
called a d-cocycle. The space of d-cocycle is denoted as 
^ d [Gf, Ut(1)]- The fermionic path integral obtained 
from a d-cocycle (i/ d , n d -i, w d _ x ) (via v d (go, — , gd) and 
m d - 2 (go, ...,g d _ 2 )) will be called a fermionic topological 
non-linear cr-model. 

In the following, we will show that we can construct 
an exactly soluble local fermionic Hamiltonian in d spa- 
tial dimensions from each (d + l)-cocycle (v d +i, n dl u 9 d ) 
in £F d+1 [Gf, Ut(1)]. The Hamiltonian has a symme- 
try Gf. The ground state wave function of the con- 
structed Hamiltonian can also be obtained exactly from 
the (d + l)-cocycle. Such a ground state does not break 
the symmetry Gf and describes a fermionic SPT phase. 

VI. IDEAL GROUND STATE WAVE FUNCTION 

We have shown that from each element (f d +i, n dl u 9 d ) 
of 2f d+1 [Gf, C/r(l)], we can define a fermionic topologi- 
cal non-linear cr-model in (d+ 1) space-time dimensions. 
The action amplitude V d+1 of the model [obtained from 
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FIG. 10: A 9-site triangular lattice on a torus, where each 
site i has physical states \gi) where gi € Gb, and each triangle 
has fermionic states \riijk) where riijk = 01 is the fermion 
occupation number. The orientations on the edges give rise 
to a natural order of the three vertices of a triangle (i, j, k) 
where the first vertex i of a triangle has two outgoing edges 
on the triangle and the last vertex k of a triangle has two 
incoming edges on the triangle. The triangular lattice can be 
viewed as the surface of solid torus. A discretization of the 
solid torus can be obtained by adding a vertex-0 inside the 
solid torus. The branching structure of the resulting complex 
is indicated the arrows on the edges. 



pointing from vertex-0 to the vertex on the surface. We 
note that the simplices associated with the down trian- 
gles in Fig. 10 have a "+" orientation while the simplices 
associated with the up triangles have a "— " orientation. 

From each 3-cocycle (z/3, U2,u 2 ) in ^ 3 [G/, J7y(l)], we 
can construct a fixed-point action amplitude. For sim- 
plices with "+" orientation, the fixed-point action am- 
plitude is given by 

V3" (50 ,9i,9j,9k) = "3" (.9o , 9i 1 9j , 9k ) x 

a n 2 (gi ,9j ,9k ) n n 2 (go ,9i ,gk ) n n 2 (go ,9j ,gk ) ff 1 ? (9o,9i,gj) 1 r 1 \ 
"(123) "(013) "(023) "(012) \ 01 ) 

and for simplices with "— " orientation given by 

V3" {9o,9i,9j,9k) ="3 (.9o , 9i , 9j , 9k ) x 

/f»2 (go ,9* ,gj ) n n 2 (go ,9j,gk)nn 2 (90,9i ,9k ) /p2 (9i ,9j,9h) /r \ 
"(012) "(023) "(013) "(123) ' V 0Z J 

where 



( v d+ii n d, u d)\ * s a fixed-point action amplitude under the 
coarse-graining transformation of the space-time com- 
plex. The fermionic path integral is supposed to give 
us a quantum ground state. We claim that such a quan- 
tum ground state is a fermionic SPT phase described 
by (fd+ii n di u d)- I n the section, we will construct the 
ground state wave function, in (2+l)D as an example. 



A. Construction of 2D wave function 

We will assume that our 2D system lives on a trian- 
gular lattice which forms a 2D torus (see Fig. 10). On 
each lattice site i, we have physical states \gi) labeled by 
gi € Gb- On each triangle (ijk), we have two states: no 
fermion state |0) and one fermion state |1). 

The ideal ground state wave function can be obtained 
by viewing the 2D torus as the surface of a 3D solid 
torus. The fermionic path integral on the 3D solid torus 
with the action amplitude V^ +1 obtained from the cocy- 
cle (fd+ii "d, Uj) should give us the ground state wave 
function for the corresponding fermionic SPT state. To 
do the fermionic path integral, we need to divide the 3D 
solid torus into a 3D complex with a branching structure. 
Due to the topological invariance of our constructed ac- 
tion amplitude, the resulting wave function should not 
depend on how we divide the 3D solid torus into 3D com- 
plex. So we choose a very simple 3D complex which is 
formed by the triangular lattice on the surface of the 3D 
solid torus and one additional vertex-0 inside the 3D solid 
torus. The resulting 3D complex is formed by simplices 
[Oijk], where ijk is a triangle on the surface (see Fig. 
10). The branching structure of the 3D complex is given 
by the orientations on the edges. Those orientations for 
the edges on the surface are given in Fig. 10. For the 
edges inside the 3D solid torus, the orientation is always 



^3(9o,9i,9j,9k) = (-)" ll(9o * : '1/3(50, 9i, 93,9k), 

v 3 (9o, 9i, 9 V 9k) = (-)" llfe ' 9fc) / v 3 (g , 9i,9j, g k ), (53) 

Since "2(30, 31, 92) is in Ti 2 (G, Z2), we can always write 
"2(50,31,32) in terms of mi(# , 3i) : 

"2(30,31,32) = m 1 (g 1 ,g 2 ) + mi (30,32) + ™i(3o,3i) mod 2, 
mi (go, 3i) = 0,1. (54) 

However, in general m 1 (gg Q ,gg 1 ) ^ mi (30,31), 
even though "2(30,31,32) satisfies n 2 (gg , 331, 332) = 
"2(30,31,32)- 

Now, the wave function is given by 

*({3*},{%fc)},{%fc)}) 

= / n dff (SSr ,<M,iw)dg ss D,iM,w) nH mi(flo,! " ) x 

(Oij) {0i} 

, 9i, 9], 9k) Y[ V 3 (9o, 9i,gj, 9k), (55) 

A V 

where Il(oij) * s the P r °duct over all links of the triangular 
lattice, Il{oi} ^ s tne P r °duct over all sites of the triangular 
lattice, YIa ^ s the product over all up-triangles, and J| v 
is the product over all down-triangles. Clearly, the above 
wave function depends on gi through 1/3, n 2 and mj. The 
wave function also appears to depend on g 0l the variable 
that we assigned to the vertex-0 inside the solid torus. 
In fact, the go dependence cancels out, and the wave 
function is independent of go- We can simply set go = 1. 

If all the mi dependence also cancels out, 
the wave function will have a symmetry de- 
scribed by Gf. *({33j,{%jfc)},{%fc)}) = 
* s(9) ({3i},{%jfe)},{%jfe)})- But does the m x de- 
pendence cancel out? Let us only write down the mi 
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dependence of the wave function: 

*({&}> {%*)}-{%*)}) 
= J^J (_)™i(so,g;) TT(_)tni(ffi,flfc) TT(_)mi(so,3i) 

{Oi} A v 



(56) 



A 



where ijk around the up- and down-triangles are ar- 
ranged in a way as illustrated in Fig. 10. We also have 
used the relation U{ 0l }(-) mi{9o ' 9z) n v (") mi(9o * ) = L 
We see that mi does not cancel out and the wave function 
is not symmetric, since mx(gg ,ggi) ^ mx(g ,gi). 

This is a serious problem. But the symmetry breaking 
is only on the surface and it can be easily fixed: we simply 
redefine the wave function by including an extra factor 
n u p-icft(-) mi(9l ' 9fc) on the surface: 



{*(«*)}. {%*)}) 



(57) 



"2 (So ,ffi ,gj) ^ff 1 ^ (9a ,S»,3j) 



(Or?) 



(Oij) 



= n (-) mi(3i ' 9k) / n de 

up-left J (Oij) 

II R mi(ff0 ' 9l) n v * (do, gi,9i,gk) II V 3 + (<? , gi, g h g k ) 

{Oi} A V 



n ^oij) ,9t * ) d ^™ 2 -° ,g% ' 3i ] x 



(Oij) 



(Oij) 



J| v z 1 (go, 9t, 9j, 9k) Y[ v a(9o, 9i, 9j,9k) 



A 



,"2 (SO iSi iff j ) fl "2 (ffO,Sj ifljs ) a"2 (ffO,ffi ,gk) O n2 (Bi >9j iSfe) 



'(Ojfc) 



(Oi/c) 



A 



V 



"2 (si ,gj ,gk ) o n 2 (go ,s; ,gh) (so ,3? ,gh) n n 2 (so ,ff; ,3j ) 

P (Oifc) >jfc) °(0ij) 



where Ilup-ieft ^ s a product over all the links with 
the up-left orientation (note that Ilup-icft(-) mi(9, ' fffc) = 
n A (-) mi(ff,,fffc) and see Fig. 10). The redefined wave 
function is indeed symmetric, and is the wave function 
that corresponds to the fermionic SPT state described by 
the cocycle (v 3 ,n 2 ) G 3? 3 [G f , U T (1)]. 

We note that the wave function 
{%jfc)M0(ijk)}) described above contains 
Grassmann numbers. Indeed, ^({gi}, {0(ijk)}i {^(ijk)}) 
can be regarded as the wave function in the fermion 
coherent state basis. After we expand the wave function 
in power of the Grassmann numbers, we obtain 



*({$<}.{%*)},{%*)}) 



E 

n ijk =0,l 



A 



Then $({gj}, {nyfc}) is the amplitude of the ground state 
on the fermion number basis ®i\gi) ®(ijk) \ n ijk), where 
|<7j) is the state on site-i and \fiijk) is the state on tri- 
angle (ijk) (where riijk = 0, 1 is the fermion occupation 



number on the triangle (ijk)). Note that the sign of 
^({9i}i i n ijk}) will depend on how the Grassmann num- 
bers are ordered in U.A^ijk)Us7 d my We also notc 
that the wave function vanishes if mjk 7^ n-z(gi,gj,gk)- 

In Appendix C, we study the fermionic cocylces 
( v d+\i n d, u d) systematically. In particular, we will 
study the equivalence relation between them. This 
will lead to the notion of group super-cohomology 
class Jf d+1 [Gf,U T (l)]- In Appendix Gl, we will 
show the elements in J^ d+1 [G /,C/t(1)] correspond to 
distinct fermionic SPT phases. We will calculate 
Jf d+1 [Gf, Ut(1)] for some simple symmetry groups G/ 
in Appendix F. This allows us to construct several new 
fermionic SPT phases. The results are summarized in 
Table II. 

In Appendix H, we further show that for each ground 
state wave function derived from (vd+i, nd, u?), one can 
explicitly construct an exact solvable hermitian (local) 
Hamiltonian H — — Hj. Similar as the bosonic case, 
Hi acts on site i as well as its 6 neighbor. However, since 
there are 6 more fermionic degrees of freedom in the 6 tri- 
angles surround i, Hi also acts on these 6 triangles. More- 
over, all the Hi are hermitian (unconstraint)commuting 
projectors satisfying Hf = Hi and HiHj = HjHi. Such 
a nice (frustration free)property make it exact solvable, 
with ty({gi}, {Oujk)}, {0(ijk)}) as its unique ground state. 



B. No intrinsic topological orders 



In this section, we are going to show that the wave 
function constructed in the last section $({<7j}, {n^/.}) 
contain no intrinsic topological orders as a fermion sys- 
tem. In other words, starting with the following pure 
bosonic direct product state 



|*o) =<8>i|& 



fa) = \G b r 1/2 E 1* 

SiGGb 



(59) 



we can obtain the fermionic state 1$) constructed in last 
section after a fermionic LU transformation defined in 
Ref. 25. 

To show this, we start with an expand fermionic 
Hilbert space, where we have four fermionic orbitals 
within each triangles. We also have bosonic state \gi) 
on each vertex (see Fig. 11). In the Grassmann number 
form, the constructed fermionic wave function is given by 
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FIG. 11: (Color online) A triangular lattice with four 
fermionic orbitals in each triangle (ijk). The fermions on 
the two solid dots in the yellow triangle are described by op- 
erators C(yfc) (the center one) and C(oifc) (the side one), and 
the fermions on the two open dots are described by operators 
croij) and CfojfcV The fermions on the two open dots in the 
blue triangle are described by operators cujk) (the center one) 
and C(oifc) (the side one), and the fermions on the two solid 
dots are described by operators C(oij) and C(ojk)- Each vertex 
has bosonic states described by \gi) , g% G Gb- 



eqn. (60): 

{6(ijk)}, {%*)}) 



(Oij) 



(60) 



II ^3 1 {9o, 9%, 93,9k) n ^3(50, 9i, 9j,9k)x 

A V 

n/)"2 (go ,9i ,9j ) n n 2 (go ,9j ,9k ) n n 2 {go ,gi,gk)n n 2 (gi ,gj,9k) 
°(0ij) "(Ojk) U (0ik) °(ijk) 



A 



n 



q«2 (gt ,gj ,sk ) fl "2 (90,9* ,gk) a n 2 (go ,gj ,gk) a n 2 (go ,gi ,Sj ) 



(ijk) 



(Oik) 



'(Ojk) 



(Oij) 



We would like to point out that although the above ex- 
pression contains 50 , the topological invariance of the 
fermion path integral ensures that the Grassmann wave 
function on the left-hand-side does not depend on go- 
lf we treat 9's and #'s as the following complex fermion 
operators 



\abc) — c (aoc)' ® (abc) ~ c (abc) ' 



(61) 



d9(abc) = C(abc), d#(abc) = C( a bc) , 

the expression eqn. (60) can be viewed as an operator: 

(62) 



n 2 (go,gi,gj)-n2(go,gi,gj) 

C (0ij) C (0ij) 



(Oij) 

]J v 3 1 (-9o . 9i > 9 3 9k ) Y[ "3 (9o , 9% , 9 j , 9k 

A V 

n 



„t™2 (9o ,9i ,9j )J n 2 (go,9j ,9k ) -t»2 (9o ,9* ,9k ) =t"2 (gi ,9j ,9k) .. 



-(Ojfc) 



'(ijk) 



A 



n 



t»2 (9i ,93 ) t«2 (go ,g% ,9k ) ^t«2 (go ,9j ,9* ) -t«2 (go ,9< >9j ) 
C (ijfc) C (0ifc) C (0jfc) C (0ij) 



Again U is independent of 50 1 despite the appearance 
of go on the right-hand-side. Then the fermionic state 



constructed in the last section can be obtained from the 
bosonic product state \&o): 



I*} = U\$o) 



(63) 



Now we would like to point out that U itself 
is formed by several layers of fermionic LU trans- 
formations. Since 1/3 is a pure U(l) phase, thus 
11a u 3 1 (9a, 9i, 9j,9k) Il v ^3(30, 9i, 9j,9k) represents lay- 
ers of bosonic LU transformations. Also the U has a 
property that when acting on \&o), and & always act 
on states with no fermion and c and c always act on 
states with one fermion. So those operators map a set 
of orthogonal states to another set of orthogonal states. 
In this case, an even numbers of c's and c's correspond a 
fermionic LU transformation. Therefore, the terms in U, 
„„_>, „„ J«2 (gi ,9) ,9k ) Ana (90,9. ,9k ) -t«2 (go ,9j ,9k ) -t«2 (go ,gi ,9j ) 
suui as c.^.^ c (oik) (Ojk) (Oij) 

and c^°' m ' ai) ^°' 9i ' ai \ all represent fermionic LU 
transformation as defined in Ref. 25. Note that none 
of the above transformations changes 5;. So we can show 
those transformations to be unitary within each fixed set 
of {gi}- Therefore, the state \^f) has no fermionic long 
range entanglement as defined in Ref. 25 (ie no fermionic 
intrinsic topological order). I'F) is the fermionic SPT 
state described by a cocycle (v3,n 2 ,u g 2 ). 



VII. EXAMPLES 
A. 2D fermionic SPT state with Z2 symmetry 

In Appendix F, we have calculated the group super- 
cohomology classes Jf z [Z 2 x z(, U(l)] = Zf. This means 
that interacting fermion systems with a Z 2 symmetry can 
have (at least) 4 different SPT phases: a trivial one plus 
three non-trivial ones. This result is described by the 
Z 2 x Z 2 row and d sp — 2 column of Table II. 

In the following, we will list the explicit fermionic 
3-cocycles [1/3 (g , 31,32,53), n 2 (go, 31,52)], 9% G G b = 
Gf/z£, that describe the above four phases. The first 
two SPT phases are given by 

"2(50,51,52) = 0, 

1/3 (e, 0, e, o) = 1*3(0, e, o, e) = ±1, other z^ 3 = 1. (64) 

where we have assumed that the elements in = Z 2 
are described by {e, 0} with e being the identity element. 
The two 3-cocycles actually describe bosonic SPT phases 
[the Z n row and d sp — 2 column of Table I (with n = 2)]. 
The 3-cocycle f 3 (e,o, e, o) = i/ 3 (o, e,o, e) = 1 describes 
the trivial SPT phase and the 3-cocycle vz(e, o, e, o) = 
^3(0, e,o, e) = —1 describes the non-trivial bosonic SPT 
phase, protected by the on-site Z 2 symmetry 



W g = H\gJ(ggi\, g, gi eZ 2 . 



(65) 
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If we choose the bulk Hamiltonian to be the ideal 
Hamiltonian H — J^i^i with Hi given by eqn. (H9), 
the low energy edge excitations will have a zero energy 
(zero bandwidth). So we will add the following Z 2 sym- 
metric edge Hamiltonian to give those zero energy edge 
excitations a non-trivial dynamics: 



-Hedge = y^[Jofgf +1 + haf] 



(66) 



where a\ are the Pauli matrices acting on the two states 
|e) and |o) on site-z. From eqn. (G23) and eqn. (G28), 
we find that the low energy effective edge Hamiltonian is 
given by 



H eS = J2[Ja?a? +1 +haf] 



(67) 



for the trivial SPT phase with v 3 (e,o, e,o) = 
1/3(0, e, o, e) = 1; and by 

H eS = ]T[Jafaf +1 + £of (1 - <tU°Ui)] ( 68 ) 



for both the fermionic SPT phases with ^(e, o, e, o) 
1/3(0, e, a, e) = ±i; where 

















9 



(72) 



Note that the term vanishes in the low energy edge 
sector. 

Let us call a link + 1) with cr?af +1 = — 1 a do- 
main wall. We note that all the three terms in the above 
i? e ff cannot change the numbers of domain walls. So the 
model has an effective U(l) symmetry. It has two phases: 
gapped phase for large J where a z = ±1 and there are no 
domain walls; gapless phase for large /loo where the do- 
main walls form a "superfluid" . The gapped phase breaks 
the Z 2 symmetry while the gapless phase is described by 
a central-charge c = 1 conformal field theory. 

From the expression eqn. (G22) of the low energy edge 
state, we see that the total number of fermions in the low 
energy edge states |{gj e dge) are given by 



N f 



E 



m(9a,9i+Xi9i)- 



(73) 



for the non-trivial SPT phase with i*j(e, o, e, o) = 
1/3(0, e,o,e) = -1. 

The other two SPT phases are given by 

n 2 (e, o, e) = n 2 (o, e, o) = 1, other n 2 = 0, 
i/3(e, o, e, o) = 1/3(0, e, o, e) = ±i, other 1/3 = 1. (69) 

They describe non-trivial intrinsic fermionic SPT phases. 
The ideal bulk Hamiltonian is given by H = J^i Hi with 
Hi given by eqn. (H9). 

Again, for the ideal bulk Hamiltonian, the low energy 
edge excitations all have a zero energy. Let us add the 
following Z 2 symmetric edge Hamiltonian to give the zero 
energy edge excitations a non-trivial dynamics: 



H 



edge 



[-Jofof +1 + hootf 



(70) 



+ (/i+-of c[ i(i+1) c (i _ M) + /i++<cJm+i) c (;-m) + h -°^' 

where a fermion lives on each link (i, i + 1). From 
eqn. (G23), eqn. (G28) and eqn. (G30), we find that the 
low energy effective edge Hamiltonian is given by 



H 



elf : 



E 



-0-^(1 + 0-1^(1-0-!^) 
-af(l -ot x )(l + + h - 



E - J °>l+r 



h m + Rc(/t + _) 

o i v^i—X - 



Im(h+-) 



(V? - CTU°l°i + l) 



(71) 



We note that the above sum does not depend on g a and 
is only a function of gis. If we choose g a — e, we see that 
only the domain wall (gi,9i+i) — (e, o) has one fermion, 
while the domain wall (ffj, ffj+i) = (o, e) has no fermion. 
Thus the number of fermions is equal to the half of the 
number of the domain walls. 

To consider a more general edge effective Hamiltonian 
that does not conserve the domain wall numbers, we add 
the following edge terms: 

#cdgo = ^ fep+Q+gf 0~j + 2 c ( j C (i+2.i+3-\ + h - C - 
i 

The corresponding edge effective Hamiltonian is given by 
^ iIm(/i 0+0 +) , + + _ _ w 

-ffedge = 2^ 9 «°"i+2 - °"i Cr S :+ 2 ) X 



(l + af_ l( Tf +1 )(l + a.f +1 af + 3). (74) 
So the more generic edge effective Hamiltonian has a form 

#cdgc = E [ _ Ja i a i+1 + h a i ( a i-l _ °f+l) 

+ /iiK- C 7?_ 1 ^af +1 ) (75) 
+ ih 2 (a+a+ +2 - urar +2 )(l + of_ l0 f +1 )(l + ^+i< +3 ) 

The new term can only change the domain wall number 
by ±4, which can be irrelavent. So the edge excitations 
still be gapless. 

From the structure of n 2 (go , 9 1 , 52 ) in eqn. (69), we 
find that the number of fermions on the edge is given by 
half of the domain wall number. So the fcrmion-number- 
parity conservation only allows the domain wall number 
to change by a multiple of 4. 
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We note that \ X^=i(l — a f (T i+i) coun ts the number 
of domain walls, where N is the number sites on the ID 
edge. Therefore, the following unitary operator 

(N \ N 

yX^-^+i) n°? ( 76 ) 
i=l / i=l 

satisfies S 2 = 1 in the subspace where the domain wall 
number is zero mode 4. So we can regard S as a Z 2 
transformation in such a subspace. One can also check 
that the J, ho and hi terms in eqn. (75) commute with 

both exp (^ -ofof +1 )) and ritLi °f 1 whilc 

the h 2 term in eqn. (75) anti-commute with both 

ex P (if EiliC 1 ~ °"l °f+i)) and ntLi °f • So the ed g e ef- 
fective Hamiltonian _ff (75) is invariant under S. Similar 
as the bosonic Z 2 SPT case 49,59 , it can be shown that any 
ID system satisfies such a non-local Z 2 symmetry must 
be gapless if the symmetry is not broken(spontaneously 
or explicitly). 

B. 3D fermionic SPT state with T 2 = 1 
time-reversal symmetry 

In Appendix F, we have also calculated the group 
super-cohomology classes Jff 4 [Z 2 x Z(,U(1)] = Zf. So 
interacting fermion systems with a T 2 — 1 time-reversal 
symmetry can have (at least) 4 different SPT phases: a 
trivial one and three non-trivial ones. This is described 
by the Z 2 x Z 2 row and d sp = 3 column of Table II. 

Let us list the fermionic 4-cocycles 
M50, 91,92,93, 94), n 3 (g , 51,32,53)] that describe 
the above four phases. The first two SPT phases are 
given by 

^3(5o,5i,52,53)=0, (77) 
ZAi(e, o, e, o, e) = 1/4(0, e, o, e, 0) — ±1, other 1/4 = 1. 

Again, the two 4-cocycles actually describe bosonic SPT 
phases (the Z 2 row and d sp — 3 column of Table I) . 
The other two SPT phases are given by 

ri3(e, 0, e, o) = 713(0, e, 0, e) = 1, other 713 = 0, (78) 
^(e, o, e, o, e) = —1/4(0, e, o, e, o) = ±i, other 1/4 = 1. 

They describe non-trivial intrinsic fermionic SPT phases. 

Using the above cocycles, we can write down the cor- 
responding wave functions and exactly soluble Hamilto- 
nians. However, the explicit Hamiltonian is very compli- 
cated. We wonder if there exist a better basis, in which 
the Hamiltonian will have a simpler form. 

VIII. SUMMARY 

It was shown in Ref. 50 that generalized topological 
non-linear cr-models with symmetry can be constructed 



from group cohomology theory of the symmetry group. 
This leads to a systematic construction of bosonic SPT 
phases in any dimensions and for any symmetry groups. 
This result allows us to construct new topological insu- 
lators (with symmetry group U(l) x Zj) and new topo- 
logical superconductors (with symmetry group Z 2 ) for 
interacting boson systems (or qubit systems). It also 
leads to a complete classification of all gapped phases in 
ID interacting boson/fermion systems. 

In this paper, we introduce a group super-cohomology 
theory which is a generalization of the standard group 
cohomology theory. Using the group super-cohomology 
theory, we can construct new discrete fermionic topo- 
logical non-linear cr-models with symmetry. This leads 
to a systematic construction of fermionic SPT phases 
in any dimensions and for certain symmetry groups G/ 
where the fermions form an ID representation. The dis- 
crete fermionic topological non-linear cr-model, when de- 
fined on a space-time with boundary, can be viewed as 
a "non-local" boundary effective Lagrangian, which is a 
fermionic and discrete generalization of the bosonic con- 
tinuous Wess-Zumino-Witten term. Thus the boundary 
excitations of a non-trivial SPT phase is described by 
a "non-local" boundary effective Lagrangian, which, we 
believe, implies that the boundary excitations are gapless 
if the symmetry is not broken. 

As an application of this general result, we con- 
structed three non-trivial SPT phases in 3D, for inter- 
acting fermionic superconductors with coplanar spin or- 
der. Such superconductors have time-reversal Z 2 and 
fermion-number-parity Z 2 symmetries described by the 
full symmetry group Gf = Z 2 x Z^. 44 The three non- 
trivial SPT phases should have gapless excitations on the 
2D surface if the time-reversal symmetry is not broken. 
It is known that non-trivial 3D gapped topological phases 
do not exist for such superconductors if the fermions are 
non-interacting. So the constructed non-trivial fermionic 
SPT phases are new. 

We also constructed three non-trivial SPT phases in 
2D, for interacting fermionic systems with the full sym- 
metry group Gf = Z 2 x Z 2 . We show that the three 
non-trivial SPT phases indeed have gapless excitations 
on the ID edge described by central-charge c = 1 confor- 
mal field theory, if the Z 2 symmetry is not broken. 

Clearly, more work is needed to develop a group super- 
cohomology theory for the cases when the fermions do 
not form an ID representation. That will allow us to 
construct more general interacting fermionic SPT phases. 
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Appendix A: Topological invariance of the partition 
amplitudes 



1. 1 + 1D 

Similar as we prove the topological invariance of the 
partition amplitudes that describe bosonic SPT phases 50 , 
we can check all the admissible branching 2h2 and 1 O 
3 moves for the Grassmann graded 2 cocycle function: 

(Al) 

For admissible branching 2 f-> 2 moves, we can have 
three different equations(up to the orientation conjugate 




(b) 

FIG. 12: The admissible branching 2 <-» 2 moves. 
(a)Branching moves that can be induced by a global time 
ordering. (b)Branching moves that can not be induced by a 
global time ordering. 



g 3 g 3 




(b) 



FIG. 13: The admissible branching 1 -O- 3 moves. 
(a)Branching moves that can be induced by a global time 
ordering. (b)Branching moves that can not be induced by a 
global time ordering. 

+ ->■-; — >+)■ 

J de^^de^ a2) v+( gi , 92 , ff3 )v 2 -( 5o , gi, 92) 

= J dO$° - g3) d^$ - 93) V+ (go,92,9 S )V^(g ,9i,g 3 ) 

(A2) 

j Mffi**)^^ .93) V + {gu ^ g 3)V + {g0j g u g 3 ) 

= / <ot 0,92) <02 ( r ,S2) V 2 + (3o,52,53)V 2 + ( ff o,5i,52), 

(A3) 

and 

J de^ S3) de^ 3 f' g3) v^( 9l ,9 2 , g 3 )v 2 + ( 50 , 32, ff3 ) 

= / %\ ( ) 9 " 9l) <o\ ( ) 9 " 9l)v 2 + (3o,3i,33)V 2 -( ffo ,3i,32). 

(A4) 



In this section, we will prove the topological invariance 
of the partition amplitudes Eq. (25) in 1 + ID and 2 + ID 
under their corresponding fermionic group cocycle con- 
dition. The 3 + ID case is much more complicated but t 
can still be checked by computer. 
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For admissible branching 1 H 3 moves, we can have four 
different equations(up to the orientation conjugate H — > 

-;--►+): 

V 2 + (3o,ffi,ff3) 

ja«i (90,92) ja"i (90,92) jnni (91,92) j/Jni (91,92) 
Ut7 (02) UO {02) U(7 (12) U(7 (12) 

x d^ 32 ' 93 ^^ 2 ' 935 (-r° (ff2) 

x ^2 (9i, 92, 93W2 (go, 92, 93W2 (do, 9i, 92) (A5) 

V 2 + (30,ff2,ff3) 

j/)«i(90,Si) jfl"i(90,9i) (31,92) (91,92) 

d >l) d >l) de, (12) d6l (12) 

x [gi,g2,g?,)Vt (go, gi,g3 )V 2 ~ (50, 31, 52), (A6) 

and 

V 2 + (.gi,32,ff3) 

j/)"i (90,91) j/7«i(9o,9i) j/)"i(9o,92) i/j«i (90,92) 
at, (01) QP (01) QP (02) Qt, (02) 

i/ 9 ™i(9o,93) r iHni(9o,93) ( ' \m (go) 

x V 2 + (go,.g2,5 , 3)V 2 _ (5o,gi,ff3)V 2 + (g ,3i,ff2) (A7) 

V 2 + (30,ffl,ff2) 

— / AQ ni ( go > 9a } Jfl"i(90,93) j/)»i(gi,93) j5«i(9i,93) 

- / d ^(03) dy (03) dy (13) dy (13) 

x V 2 + (gi,g 2 ,g 3 )V 2 "(3o,32,ff3)V 2 + (go,gi,g3)- (A8) 

Note that the first two equations of 2 «-» 2 and 1 <-» 3 
moves can be induced by a global time ordering(see in 
Fig. 12 (a) and Fig. 13 (a)) while the rest can not(see 
in Fig. 12 (b) and Fig. 13 (b)). Here g° is defined on 
the vertex with no incoming edge, g 1 with one incoming 
edge, etc. 

Let us use the definition of V 2 and integral out the 
Grassmann variables. For the branching 2 o 2 moves, 
we have: 

"2 (gi , 52 , S3>2 (50 , 91 , 92) ="t (go ,92,93) "2 (go ,91,93) 

(A9) 

"2 (91,92, 93)^} (go, 91, 33) =v% (g , 92,93)^2 (go, gi, 32), 

(A10) 



and 



"2 (91,92,93)^2(90,92,93) 
=(-) nii3l ' 92) 4(9o,9i, 93)^(90,91, 92). (All) 

Similarly, for the branching lf>3 moves, we have: 

4(9o,9i,93) (A12) 
=(-) m ° (92 W(5i, 92, 93)^2 (go, 52, 9a)v£ (90,91,92) 
"2(90,92,93) (A13) 

=(-) T " o(9l) f2 f (5i 5 52, 93)^2(90,91,93)^2 (go, gi, 52), 

and 

^ 2 + (.9i,52,5 3 ) (A14) 
=(-) m ° {gi) v% (50, 52,53)^(50, 91,93)^2 (90, 91,92) 

"2(90,91,92) (A15) 
=(-) mo(ff2) f2 f (5i, 52,53)^(50, 52,53)^2" (5o, 51,53)- 

If we use the definition of v 2 : 

"2(90,91,92) = "2(90,91,92), 

"2(90,91,92) = (-) mo{9l) /"2(go,gi,92), (A16) 

All the above admissible branching moves will be equiv- 
alent to a single fermionic 2 cocycle equation of ^ 2 

^2(50,51,53)^2(51,52,53) = ^2(50,51,52)^2(50,52,53)- 

(A17) 

We note that this equation is the same as the 2 cocycle 
equation in bosonic systems. 



2. 2 + ID 

In 2 + ID, there are in total 10 admissible branch- 
ing 2 O 3 moves and 5 admissible branching 1 f) 4 
moves. Let us show all theses moves will lead to the 
same fermionic 3 cocycle condition. 

For example, the admissible 2 f-> 3 in Fig. 6 represents 
the following equation for V3 : 



J &6< 9 lf 2,93) d^ 93 ' 93) V 3 + (5o, 5i, 92, 9 3)Vt(gi, 52, 53, 54) 

/r\a n (go, 91 ,94) j 7jn(go, 91,94) j/) ra (90, 92,94) Jfjniao, 92, gi) jfl«(9o, 93,94) j2)™(90, 93,94) ( _\m 1 (g ,g 4 ) 
Ut7 (014) Ul7 (014) Ut7 (024) Ut7 (024) Ut7 (034) UP (034) V ) 

x (,g , 5i, 52, 54)V 3 f (50, 52, 53, 54)Vf(5o,5l,53, 54) (A18) 
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Note that in the above expression, we put the sign fac- 
tor (— ) m i(fo,94) on interior link and integral out the 
Grassmann variables on the interior faces. To simplify 
the representation, we can formally rewrite the above 
equation as (see eqn. (18)): 



V3 (.go , 5i , 92 , 93 ) V 3 + (31 , 52, 53, 54) ( A19) 

Va (9o, 9i,92,9i)Vt (90, 92, 33, 94) V3" (go, 3i, 33, 34) 



and 



In such a way, we can formally write down all admis- 



sible branching 2 f-> 3 moves in terms of V 3 (up to the 
orientation conjugate H — ► — ; > +): 



^3 (.90 
^ (3o 

V 3 "(3o 



^3 (50 

V 3 -(3o 

V 3 + (3o 
V 3 ~(so 

V 3 + (3o 



V 3 + (3o 



V 3 (30 
V 3 ~(3o 

V 3 -(3o 



V 3 + (3o 



V 3 + (3o 



V 3 + (3o 



3i,32 
3i,32 

3i,32 
3i,32 

31:32 

3i,52 

3i,33 
3i,32 

3i,32 
3i,32 

3i,32 
3i,32 

3i,33 
3i,32 

3i,32 
3i,32 



34)V 3 (30 
33)V 3 "(3o 

9iWt (3i 
33)V 3 "(3o 



33) V 3 + (3o 
9iWti9o 

34) V 3 + (3i 

33) V 3 + (3o 

93)^3(90 

34) V 3 + (3o 

34)V 3 f (30 
33)^(30 

9iWti9o 
33)V 3 _ (3o 

33)V+(3i 



34)V 3 + (3o 



32,33,34) 



(A20) 



0i ,53, 54) V 3 (31,32,33,34) 

32,33,34) (A21) 
31,33,34)1^^(30,32,33,34) 

32,33,34) (A22) 
31,33,34)1^^(31,32,33,34) 

32,33,34) (A23) 
31,32,34)^(30,32,33,34) 

31,33,34) (A24) 
32 , 33 , 34) Vg" (31 , 32 , 33 , 34) 

3i,33,34) (A25) 
32, 33, 34)1^3" (31, 32, 33, 34) 



32,33,34) 



(A26) 



3i , 32 , 34)^ (31 , 32 , 33 , 34) 

32,33,34) (A27) 
32 , 33 , 34) Vg" (30 , 31 , 33, 34) 



V3X00, 52, 53, 54) V 3 (31,32,33,34) (A28) 
= / ^(50, 3i, 32, 33)V3" (30, 3i,32,34))V 3 f (30,31,33,34) 



V 3 (30,51,52,53)^(00,51,52,54) (A29) 
V 3 " (30 , 32 , 33 , 34 ) V 3 + (30 , 31 , 33 , 34 ) V3 (31 , 32 , 33 , 34 ) 



Similarly, we can formally write down all the admis- 
sible branching 1 f-> 4 moves in terms of V 3 (up to the 
orientation conjugate H — > — ; > +): 



V 3 + (3o,3i,32,34) = / V 3 f (3o,3i,32,3 3 )V 3 + (3o,3i,33,34) 
x ^(30,32,33,34)^(31,32,33,34) (A30) 



V 3 (50,52,53,54) = / V3 (30, 3i, 32, 33) V 3 (30,31,32,34) 
x 1^ (.90,31,33,34)1^ (31,32,33,34) (A31) 

V 3 + (5o,3i,33,34) = / V 3 ~ (30, 01, 32, 33)V 3 + (30,31,32,34) 
x V 3 f (3o,52,53,54)V 3 "(5i,5 2 ,53,54) (A32) 

and 



V 3 + (3i,32,33,34) = / V 3 (3o,3i, 32, 33)V 3 f (3o, 3i,32, 34) 
x V3"(0o,0i,0 3 ,04)V^ (30,32,33,34) (A33) 



^(00,31,32,33) = / V 3 f (0 O ,3i,52,54)V 3 (30,31,33,34) 

x 1^(30,32,33,34)1^(31,32,33,34) (A34) 

Here the symbol J" means that we put the sign factors 
(— ) J "i(fi'Si) on all four interior links and integrate over 
the Grassmann variables on all six interior faces. Note 
that the first 8 2 «-> 3 moves (Fig. 14) and first 3 1 <-» 4 
can be induced by a global time ordering while the rest 
can not. Again here 3 is defined on the vertex with no 
incoming edge, 3 1 with one incoming edge, etc. 

Expressing in terms of and integrating out all 
the Grassmann variables, the admissible branching 2o3 
moves lead to the following equations: 

"3 (5o , 3i , 52 , 54) V3 (go ,92,93,94) ( A35) 

—f_\n2(go,gi,g2)n 2 (g2,g3,g4)+mi(gi,g3) 

(50, 01, 52, 53)^ (50, 51, 53, 54)^3" (51, 52, 53, 54) 
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FIG. 14: (Color online) The admissible branching 2 o 3 
moves that can be induced by a global time ordering. 



FIG. 16: (Color online) The admissible branching 1 -o 4 
moves that (a) can be induced by a global time ordering and 
(b) can not be induced by a global time ordering 



^3(50,31,32,34)^(51,32,33,54) (A36) ^3(30,31,32,33)^(30,32,33,34) (A37) 

=(— ) "2(90 ,gu S2)n2(ff2,S3,fl , 4)+mi (90,53) _/_\"2 (go, 91, 92)^2 (92, 93, S4)+mi (91,94) 

x ^3" (So, 3i,32, 33)^3" (3o, 9i, 33, 34)^(30, 32, 33, 34) x ^(30,31,32,34)^ (30,31,33,34)^(31,32,33,34) 
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v t (.9o, 91,93,94)^3 (9i,92, 93,94) (A38) 

=(— 1 ™ 2 (30,91 ,S2)n2(92,fl , 3,ff4)+mi (90,32) 

X^cT (50, 51, 52, 53)^ (50, 51, 52, 54)^ (50, 52, 53, 54) 

"3" (5o ,9i,92,93)^(9q,9i,93,94) ( A39) 

— ( — )™ 2 ' 9l >32 )"2 (32 '93 ,94)+mi (g 2 ,94) 

x^s" (50, 9l,92,94) v t (50, 52, 93,94)^3 (5l, 52, 53, 54) 

^(50,51,52,54)^(50,51,53,54) (A40) 
_/_\"2 (go, si ,92)12(32, 93, 94)+«2 (gi, 92, S3)+na(ai, 92, gi)+mi (92,93) 

x ^(5o, 5i, 52, 53)^3" (5o, 52, 93,94)1^3 (51, 52, 53, 54) 
^ (9o,9i,93, 94)^3 (9o,92, 93, 94) (A41) 

— (__)"2 (90, 91, 92)12(92, 93, 94)+«2 (9o,9l ,92)+"2 (go, 91, 93)+mi (91,92) 
x ^(50, 51, 52, 53)^3" (50, 51, 52, 54)^(51, 52, 53, 54) 

^3" (5o , 5i ,92,93)^t(9i, 92 ,93,94) ( A42) 

\"2 (90,91, 92)n2(g2,93,94)+n2(go,ga,g4)+n2(gi,g3,94)+mi(go,94) 

x ^3" (50 , 5i , 52 , 54 ) v£ (9o , 92 , 53 , 54 ) V3 (9o , 9i , 93 , 94 ) 

and 

v t (5o , 52 , 53 , 54) H (51 , 52, 53 , 54) ( A43) 

_-^_yi2(go,gi,g2)«2(g2,93,94)+™2(go,gi,g2)+«2(gi,g2,g4)+mi(go,9i) 

x v£ (go, 51, 52, 53)^3" (50, 91,92,94))^ (90, 91,93, 94) 

»3 (50, 51 , 52 , 53)^3" (90, 51 , 52 , 54) ( A44) 

_^_yi2(go,gi,g2)«2(g2,93,94)+™2(go,g2,g4)+'i2(go,93,g4)+mi(g3,g4) 

x ^(5o, 52, 53,54)^3" (5o, 5i, 53, 54)^3" (5i, 52, 53, 54) 

We note that the sign factor (— ) m i(9*,9;,) comes from the 
definition of the fermionic path integral. Later we will see 
such a factor is very important to make all the branching 
moves to be self consistent. 

Similarly, all the admissible branching 1 O 4 moves 
lead to the following equations for vf : 

v t '(9o,9i,92,94) (A45) 
_/_\n2 (90,91, g2)n2(g2,g3,g4)+n 2 (g2,g3 ,94) 

x (_)mi{ao,fl3)-h7ii(ffi,ga)+mi( 9a ,fl 3 )+m 1 (ffs,94) I/ +( fl0)fll)fl2)ff3 ) 
* v i (9o, 9i, 53, 54)^3" (90, 52, 53, 54)^3" (5i, 52, 53, 54) 



v t (9o, 92, 93, 94) (A46) 

_^_^"2(go,gi,g2)n2(g2,g3,g4)+™2(go,gi,g2) 

x (_)roi(go,9i)+mi(gi,g2)+mi(gi,g 3 )+mi(gi,s 4 ) I ,+ ^ 0i g lj g 2 ^g 3 ) 

x ^3~(9o, 9i,92, 94)^3 (9o, 9i,93, 94)^3 (9i,92, 93, 94) 

^3(90,91,93,94) (A47) 
_f_\"-2 (go ,91,92)112 (92 ,g3,94)-Hi2(gi,g2 ,93) 

x (_)m 1 (ffa,fl 2 )+mx(g 1 ,ff a )+m 1 (g 3 ,g 3 )+mi(ff2,94) I/ -( 50j5li52)53 ) 
x^ 3 f (5o, 5i, 52, 54)^3" (50, 52, 53, 54)^3" (51, 52, 53, 54) 



and 



^(51,52,53,54) 



(A48) 



_f_\n 2 (go ,gi, 92)712 (g2,g3,94)+n2(go,gi ,94) 

x (•_)mi(9o,9i)+mi(9o,92)+mi(go,g 3 )+m 1 (g ,9 4 ) z/ -(- ff0)ffi) ff2j g ^ 

XV3 (9o, 91, 52, 54)^3" (50, 5i, 53, 54)^ (5o, 52, 53, 54) 
v t (9o, 9i, 92, 93) (A49) 

_^_\»2 (90 ,91 ,92)112 (92 ,g3,94)+n 2 (90,93,94) 

x (_)m 1 (ffa,fl4)+mx(g I , ff 4)+m 1 ( fl3 ,g4)+mi(ff3,fl4) I/ +( 50)5l)52)ff4 ) 

x ^(5o, 5i, 53, 54)^ (5o, 52, 53, 54)^3" (91,92,93,94) 
Amazingly, if we define: 

4(50,51,52,53) = (-r i(9o ' 92) ^ 3 (5o,5i,52,53), (A50) 
^3"(5o,5i,52,53) = (-) mi(9l ' 93) M(5o,5i,52,53) 

we find all the above equations are equivalent to the fol- 
lowing single equation, which is the fermionic 3 cocycle 
condition of V3. 

V3 (9l, 92, 93, 54)^3 (50, 51,53,54)^3 (50, 51,52,53) = 
( _ ) «2(9o,gi,g2)«2(g2,g3,g4) l , 3(ff0jff2;53j54)z/3(50]5ijff2j34) 

(A51) 



3. 3+ ID 

The admissible branching moves for 3+ ID case is much 
more complicated. There are in total 10 3 -f-> 3 moves, 
15 2 O 4 moves and 6 lf>5 moves. 

For example, Fig. 8 represents one admissible branch- 
ing 2 -f-)- 4 move, which leads to the following equation 
for Vf: 
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d ^(1234) dy (1234) 4 Iffl ) 52 , 53 , 34, 35 J V 4 (30, 3l, 32, 33, 34) 

j/)«(9o,9l,92,g5) j/p(flo ,91,92 ,95) j/)«(so,9l,S3,S5) j/)"(ffo,Sl,S3,95) j/j™(9o,9i ,94,95) j n™(9o,gi ,94 ,9s) i fl "(go ,92 ,g3,gs) i /p(go ,g2 ,g3 ,95) 
U(7 (0125) Ut7 (0125) U(7 (Q135) Ut7 (0135) UP (0145) Ut7 (0145) Ut7 (0235) U(7 (0235) 

V A() n ( 90 ' g2 ' 94 ' 9 ^ Afl n ( 90 ' 92 ' 94 ' 9 ^r\() n ( 90t93 ' 9i ' 9 ^ r\fl n ( 9a ' 93 ' 9 ' 1 ' 95 ) ™2 (90 ,91 ,95 ) +™2 (90 ,92 ,95 ) +™2 (90 ,93 ,95 ) +™2 (90 ,94 ,95 ) 

A UI7 (0245) U(7 (0245) Ut7 (0345) Ut7 (0345) ^ J 

x V/ (30, 32, 33, 34, 35)1^(30, 31, 32, 34,35) V4 (go, 91, 33, 34, 95W4 (3o, 3i, 32, 33, 35) (A52) 

Here we integrate out the Grassmann variables on the interior tetrahedra and put the sign factor on the interior 
surfaces. We can formally denote the above equation as (see eqn. (18)): 

V4 (31 , 32 , 33 , 34 , 35 ) V4 (30 , 3i , 32 , 33 , 34 ) 

Vt (3o, 32, 33, 34, 35) (30, 31, 32, 34, 95W4 (30, 31, 33, 34, 35)V^ (30, 31, 32, 33, 35) (A53) 

Similarly, Fig. 9 represents one admissible branching 1 •<-> 5 move and can be formally written down in terms of 
as: 

V4~(3o,3i,32,33,34) 

V 4 f (.go,32,33,34,35)V 4 f " (3o,3i,32,34,35)V^(3o,3i,33,34,35)V z f (3o,3i,32,33,35)V4r(3i,32,33,34,35) (A54) 

We note that here the symbol J means integrating over all Grassmann variables on 10 interior tetrahedra and put 
the sign factor (— ) m 2(9i>9j>9fc) on \q interior surfaces. 

In addition, there are 3f>3 moves. For example, one of such moves gives rise the following formal equation for V 4 : 

V 4 f (3i,32,33,34,35)V 4 + (3o,3i,33,34,35)V 4 h (3o,3i,32,33,35) 

^(30,31,32,33,34)^(30,32,33,34,35)1^/(30,31,32,34,35) 

(A55) 

In this case, the symbol J on both sides mean integrating over all Grassmann variables on 3 interior tetrahedra and 
put the sign factor (— ) m 2(g;,9j,9fc) on \ interior surface. 

After integral over all the Grassmann variables, we can express the above three equations in terms of v±: 

v t(9\, 92, 33, 34, 35)^4" (So, 3i, 32, 33, 34) 

_/\n3(ffo ,91, 92t93)»3 (90,93, 94.95)+n3(gi ,92 ,93,94)n3(9D ,91 i94,g5)+n3 (92, 93.94,g5)n3 (90,91 l92,95)+n3 (90 ,91,94,95) 
x ^_^m2(goi9i,95)+m2(9o,92,g5)+Tn2(go,93,95)+m2(9o,94,g5) 

x ^(30,32,33,34,35)^(30,31,32,34,35)^(30,31,33,34,35)^(30,31,32,33,35), (A56) 

^(30,31,32,33,34) 

— /_ \n 3 (go,9l, g2,93)n. 3 (9o,93,94,g5)+n3(gi,92,93,94)n3(go,9i,94,9 , 5)+n3 (92, 93, 94, 95 )«3 (go, 91, g2 ,gs)+«3 (90,91 ,94,g5)+«3(gi,92,94,95)+«'3 (92,93,94,95) 
x (^y^ (go, 91 , 95 )+m2 (go, 92, 9s)+m 2 (go, 93, 9s)+m 2 (go, 94, 9s)+m2(gi, 92,95)4-^2(91,93,95)+^^ 

x "4(90, 32, 33, 34, 35)^(30, 3i, 32, 34, 35V4T (So, 3i, 33, 34, 35)^(30, 3i, 32, 33, 35)^4" (3i, 32, 33, 34, 35), (A57) 

^(31,32,33,34,35)^(30,31,33,34,35)^ (3o, 31,32,33, 35)(-) m2(5l ' 93 ' 95) 
= ^_yi3(go,9i,g2,93)n3(9o,93,94,95)+n3(9i,g2,g3,94)n3(9o,9i,g4,g5)+n3(g2,93,g4,g5)n3(go,gi,g2,g5)+m2(go,g2,94) 

x ^(30,31,32,33,34)^(30,32,33,34,35)^(30,31,32,34,35)- (A58) 
Surprisingly, if we define z/^ as: 

^4 + (3o,3i,32,33,34) =(-) m2(9o ' 9l ^ )+m2(! ' 1 ' 93 ^ )+m2(9l ' 92 ' 93) ^4(30,3i,32,33,34), 
^(3o,3i,32,33,34) =(-) m2(9o ' 92 ' 94) /^(3o,3i,32,33,34), (A59) 
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we find all the above three equations will be equivalent to the following single equation of 1/4, which is the fermionic 
4 cocycle condition: 

"4,(91, 02, 03, 04, 05)^4(00, 01, 03, 04, 55)^4(50, 01, 02, 93,95) (A60) 

_ r_^n3(90,9i,fl , 2,ff3)n3(90,93,94,S5)+n3(si,S2,S3,ff4)n3(so,Si,S4,ff5)+n3 (92,33,34, Ss)«3 (90,91 ,32 ,35) x 
^4(0O,02,03,04,05)^4(0O,01,02,04,05)^4(0O,01,02,03,04)- 

Indeed, it can be verified by computer that all other admissible 3 ■<-> 3, 2 4 and lf>5 branching moves will give 
rise to the same fermionic 4 cocycle condition! Again, the phase factor (— ) m 2(9<,9j,9fe) m the fermionic path integral 
is crucial for the self consistency of all admissible branching moves. 



Appendix B: The fixed-point action on a closed complex and the unitary condition 

In the following we will show that J V^VJ is always equal to 1. 
In 1 + ID, we have: 

V2 (00 , 0i , 02 ) V 2 ~ (,g , 0i , 02 ) 

Jfl»l (9o,9l) jflil (90, 9l) J nil (90.92 ) Jfl«l (90,92) (9l,92 ) Jfl™l (9l, 92) / \m (9o)+"io(gi)+mo(92) 

Q "(01) Q "(01) Q "(02) Q "(02) Q "(12) Q "(12) ( > 

v n n \fl™i(9l,92) fl ni(9o,9l)fl".l(9o,92) -/ \ n™i (90,92) /j™i(9o,9i)n™i (91,92) 

X v 2 (00, 01, 02 J^) t/ (01) t/ (02) ^2 (00, 01, 02 J f( 02 ) f( 01 ) f( 12 ) 

Jfl"l(9o,9l) rlfl n i (90,9l) (90,92) (90,92) (91, 92) j 311 (91.92) / \ 711(90,92) 

a "(01) Q "(01) Q "(02) Q "(02) Q "(12) Q "(12) ( J 

v fl"i (9i ,92) (9o,9i) H"i (90,92) (90,92) Hni(g ,gi)Hni(gi,g2) 
X "(12) "(01) "(02) "(02) "(01) "(12) 

= 1 (Bl) 

In the above calculation, we can move the pair d6?$°' 9l) d6^? ' 9l) to the front of O?^ ' 92) 6?^ ' 92) with- 
out generating any signs. We then can evaluate fd^^ ,9l) 6^ l ^ ,Bl) S^ 0,9 ' ) G^ a,9 ' ) = 1. We next 
move de?^f ' 9l) d6?^ ' 9l) to the front of 6?$°' 9l) d?$°' 9l) without generating any signs. We then evalu- 
ate / d^^f ,sri } d^™ 1 /^ 5 ^p^ 531 } ^A 30 ) = (-)"i(9»>9i), which has a sign factor. The Grassmann inte- 
gral fdtf$ 1 '» ) df$f 1 ^ ) 6ffi M) §ffi' a ' ) = (-)»i(9i,92) also generate a sign factor. Such two sign factors 
(90, 9i )(-)"•! (91 ,92) ca ncel the sign factor (90,92) due to the condition eqn. (10) on n d . 

Similarly, in (2 + 1)D and in (3 + 1)D, we find that the sign factors in the integration measure (18) and in the 
expressions (35) or (35), when combined with those from exchanging the Grassmann numbers, just cancel each other: 

V3" (00 , 01 , 02 , 03 ) V 3 " (g , 01 , 02 , 03 ) 

j/)"2(go,gi,92) j/7«2(go,9i,92) ,/in2(go,gi,93) ^^"2(90,91,93) j/112 (90,92^3) j/p2(go,92,93) j/ina (91,92^3) j/p2(gi,g2,g3) 
U "(012) U "(012) U "(013) U "(013) U "(023) U "(023) U "(123) U "(123) 

x r_\mi(go,gi)+mi(9o,92)+mi(go,93)-l-mi(gi,g2)-l-TOi(gi,g3)+'mi (92,93) 

V ,,+ tr, r, r, n \ fl™2 (9l ,92 ,93 ) /)»l2 (90 ,gi ,93 ) fl™2 (90 ,92 ,93 ) n™2 (90 ,9l ,92 ) 

X v 3 (00,01,02,03^(123) y ( 13) "(023) "(012) 

v r, „ r , N ifl n 2(90,9l,92)/)™2(90,92,93)/i™2(90,9l,93)fln2(gi,g2,93) 

x "3 {9o, 9i,92, 93)0 (m2) fc» (023) y (013) « (123) 

j/)"2(9o,9l,92) j/)"2 (go, 91,92 ) j/i»2 (go, 91, 93) j n n 2 (90,91 ,93) in n 2 (90, g2, 93) j/p2 (90, 92,93) j/)»2 (gi.92,93) 1/112(91,92^3) 
U "(012) U "(012) U "(013) U "(013) U "(023) U "(023) U "(123) U "(123) 

v (_ V l 2 (9o. 9i -93)+™2 (91, 92,93) z)"2(gi, 92 , 93)/jn2(go, 91 ,93)/i«2(9o ,92, g3)/T"2 (go, 91,92) (fi? (go, 91, 92) a n 2 (90,92, g3)h n 2 (go, 91, 93) /)"2 (91,92,93) 
A V J "(123) "(013) "(023) "(012) "(012) "(023) "(013) "(123) 

1; (B2) 
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Vt (#0, 91,92, 93, 9i)V 4 (30, 9i, 92, 33, 34) 

j /)"3 (so, Si ,32, S3) j n n 3 (go ,31,32,93) j/)™3 (so, 91, S2.S4) j/)«3 (30,31,92,34) j/)«3 (go, Si ,93,94) -in n 3(go, 91,93, 94) 
Ut7 (0123) Ut7 (0123) Ut7 (0124) Ut7 (0124) U(7 (0134) Ut7 (0134) 

v ^^"3 (30, 32, S3, S4) 1^13(30,32,33,34) i fl ™3 (31, S2,S3,34) ifl"3(si, 32 ,33, S4) / _ \ m 2 (90 ,9i ,92 )+in 2 (go ,9l ,93 ) 
X U(7 (0234) Ut7 (0234) Ut7 (1234) Ut7 (1234) ^ J 

x ^_^m2(go,9l,94)+ni2(go,92,g3)+m2(go,92,g4)+ni2(go,93,g4)+ni2(gi,92,93)+ni2(gi,92,g4)+ni2(gi,93,94)+ni2(g2,93,g4) 

v ^, ^, r, r, \Q n 3(gi,g2,g3,g4)nn3(go,gi,g3,g4)n n 3(go,gi,92,g3)n n 3(go,g2,g3,g4)n n 3(go,gi,g2,g4) 

x "4 \90,gi,92,g3,g4)O( 1234 ) ^(Om) P (0123) P (0234) P ((H24) 

v n „ „ , ifl™3(go, 31,32, g4)/in3(go,g2,g3,34)fln3(3o, gi, g2,33)/)«3(3o, 3i,33, g4)/)"3(gi,g2,g3,g4) 

X ^4 (90,91,92, 93, S4jfc' (01 24) ^(0234) ^(0123) W (0134) ^(1234) 

/J/)«3(30,3l,32,33) J /T™3 (30 ,31 ,32, 33) J /)«3 (SO, 31 ,32,34) J/)™3 (SO, 31 ,32,34) J/)™3 (30, 31 ,33, 34) J /)"3 (SO, 31 ,33, 34) 
Ut/ (0123) Ul7 (0123) Ul7 (0124) Ut7 (0124) U(7 (0134) Ut7 (0134) 

V ^^"3 (SO, 32 ,S3,S4) Jfl" 3 (30, 32, 33, S4) 1/113 (31, 32,33,34) Jfl"3(Sl,32,33,S4) I _ \ n 3 (g ,31 ,32 ,34)+"3 (30 ,32 ,33 ,34 ) 

X U(7 (0234) Ut7 (0234) Ut7 (1234) Uf7 (1234) ^ J 

v n n 3(gi,g2,g3,94) o n 3(9o,9i,g3,g4) Qn 3 (g ,gi,g2,93) Hn3(go,92,g3,g4) n n 3(go,gi,92,g4) 
X P (1234) P (0134) (7 (0123) t7 (0234) (7 (0124) 

v /i«3 (So, Si, 32, 34) n n 3 (So, 32, 33, 34) n n 3 (go, Si ,32, 33) n n 3 (go,9i,93 ,34) n n 3 (31 ,32, 33, 34) 
A t7 (0124) I7 (0234) (7 (0123) I7 (0134) (7 (1234) 



(B3) 



Appendix C: Group super-cohomology 

In Section V, we studied the condition on v d (go,9i, •••) 
so that the fermionic path integral constructed from 
v d{go,9i, • ■•) correspond to a fermionic topological non- 
linear cr-model which is a fixed point theory under the 
RG flow. Those conditions on v d (go, 9i, •••) actually de- 
fine a generalization of group cohomology. 50 ' 51 We will 
call such a generalization group super-cohomology. In 
this section, we will study the group super-cohomology 
in detail. 

A d-cohomology class of group super-cohomology is 
a set that depend on a full symmetry group Gf and a 
fermionic G-module M. We will denote the fermionic 
d-cohomology class as M' d (Gf,M). 

We note that a fermion system always has a Z 2 sym- 
metry which corresponds to the conservation of fermion- 
number parity. Such a symmetry group is denoted as 
Z 2 which is generated by Pf = (— ) Nf where Nf is the 
fermion number. The full symmetry group Gf always 

contain Z 2 as a subgroup. We will define G& = Gf/Z 2 . 

1. Graded structure of a group Gb 

First, let us introduce the graded structure of a group 
Gb- A rfD-graded structure of a group Gb is an integer 



function 71^-1(31, ...,gd) of d variables, whose values are 
0, 1. Here we choose gi to be an element of Gb, rather 
then an element of the full symmetry group Gf. This is 
because gt correspond to the bosonic fields in the path 
integral, and the bosonic field is invariant under Z 2 . 
The function rid-i{gi, ■■■,gd) satisfies 



nd-i(99i,-,99d) =n d - 1 (gi,...,g d ), V3 e G b , (CI) 

d 

^n-d^xigo, .., gi-x,g i+1 , ..,g d ) = even, Vg ,..,g d eG. 

t=0 



We see that a c?D-graded structure of a group Gb is a 
(d — l)-cocycle n d -\ e Z d_1 (G&, Z 2 ). If a function n d -\ 
only satisfies the first condition in eqn. (CI), then it will 
be called a (d — l)-cochain. The space of all (d — 1)- 
cochains is denoted as C d_1 (G&, Z 2 ). 

The coboundary is given by 



B d 1 (G b ,Z 2 ) = -,g d ) =^2m d _ 2 (gi,-,gi-i,gi+i,-,gd) mod 2 

i 

m ' d -2{ggi, -;ggd-i) = m d _ 2 (gi,—,g d -i), Vg,gi, -,9d-i e G b 
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We say that two graded structures differ by an above 
coboundary are equivalent. Thus different classes of 
dD-graded structures are given by "^^(G;,, Z 2 ) = 
Z d - 1 (G b ,Z 2 )/B d - 1 (G b ,Z 2 ). 

2. G-module 

For a group Gf, let M be a G/-module, which is an 
Abelian group (with multiplication operation) on which 
Gf acts compatibly with the multiplication operation 
(ie the Abelian group structure): 

g-(ab) = (g-a)(g-b), geG f , a,beM. (C2) 

For the most cases studied in this paper, M is simply 
the U(l) group and a an U(l) phase. The multiplica- 
tion operation ab is the usual multiplication of the C/(l) 
phases. The group action is trivial: g ■ a — a, g G Gf, 
a € t^(l). We will denote such a trivial G/-module as 
M = U(l). 

For a group Gf that contain time-reversal operation, 
we can define a non-trivial G/-module which is denoted 
as Ut(1)- Ut(1) is also a U(l) group whose elements are 
the U(l) phases. The multiplication operation ab, a,b € 
J7t(1), is still the usual multiplication of the U(l) phases. 
However, the group action is non-trivial now: g-a — a s ^ 9 \ 
g € Gf, a 6 Ut(1)- Here s(g) = 1 if the number of time- 
reversal operations in the group operation g is even and 
s(g) — — 1 if the number of time-reversal operations in g 
is odd. 



3. Fermionic d-cochain 

A fermionic d-cochain is described by a set 
of three functions u d (g , ... ,g d ), n d _i(fl<i, g d ), and 
~,9d)- — , 9d) is (d - l)-cocycle 

in Z d ^ 1 (Gb, Z2) which has been discussed above. 
v d (go, •■•,5d) is a function of 1 + d variables whose value is 
m a G-module M, v d : G\ +d -> M. A gain, note that g t 
is an elements of Gf, rather then the full symmetry group 
Gf. For cases studied in this paper, M is always a £7(1) 
group (ie v d (go, ■ ■■■,9d) is a complex phase). 

Since the action amplitude is invariant under the sym- 
metry transformation in Gf, v d (go,---,gd) must satisfy 
certain conditions. Note that a transformation g in 
Gf will generate a transformation in Gf,, gGb — > Gb, 
since Gb — Gf/Z%. So v d (go, ...,g d ) transform as 

v d (go,-,9d) -> v" d i9 \gg ,...,gg d ). 

The invariance under the symmetry transformation in 
Gf is discussed in Sections VE and VIA. This motivates 



us to require that the fermionic d-cochain v d (go, g d ) 
to satisfy the following symmetry condition 

v s d [9 \ggo,-,ggd) 
= Vd(go,-,9d) X\[u d ^ 1 {go,-.,gi,.-,g d )] { ~ y , (C3) 

i 

nd-i(ggo, -,99d-i) = n d -i(g , ...,g d -i), g € G f , 

where the sequence go , . . , gi , . . , g d is the sequence go, ■■■,g d 
with gi removed, and u d _ x (gi, gj, gk) is an ID repre- 
sentation of Gf 

u d a L- L (ggi,ggj,...,ggk) = [u 9 d -i(9i,9j, ■■■,9k)] s{9) , (C4) 

such that ^ 1 ( 5l , 5 ,,..., 5fc ) = (_)"*-iC»i.»i,-,»»). The 
triple functions (v d , n d -\, v? d _^) that satisfy eqn. (C3) 
and eqn. (C4) are called fermionic cochains. We will 
denote the collection of all those fermionic cochains 
(v d ,n d ^,u 9 d _ x ) as^ d [G f ,U T (l)}. 

4. Fermionic d-cocycle 

With the above setup, now we are 
ready to define fermionic d-cocycle. The 
fermionic d-cocycles are fermionic d-cochains 
Wd(go, ■■,gd),n d ^ 1 (g 1 , ..,g d ),u 3 d _ x (gi, ..,g d )} in 
< io d [Gf, Ut(1)] that satisfy some additional condi- 
tions. Fermionic 1-cocycles are fermionic 1-cochains 
M5o,5i),«-o(5i), u o(5i)] that satisfy 

vi(9i, 52)^r 1 (5o, 52)^1(50, 5i) = 1 (C5) 

Note that no (gi) = always. Fermionic 2-cocycles are 
fermionic 2-cochains [^(go, 9i, 92),n 1 (g 1 , g 2 ),u 9 1 (g 1 , g 2 )} 
that satisfy 

v2{gi,g2,g?,)v 2 - 1 (go,g 2 ,gz)v 2 (go,gi,g?,)v 2 1 {go,gi,g2) = 1 

(C6) 

Fermionic 3-cocycles are fermionic 3-cochains 
Mscfll, 92, gz),n 2 {g\, g 2 , . g 3 ), w|(5i> 92, 53)] that satisfy 

( _ ) « 2 ( ff 0, 9l , 92 )n 2 ( 92 , 33 , ff4 )^ (5ij52;53;54)l/ ~l (50j52;53;54)x 

^3(90, gi, g3, g^v^ 1 (go, gi, g2, g^sigo, gi, 92, 93) = 1 

(C7) 

Fermionic 4-cocycles are fermionic 4-cochains 

Wi (go, 91, 52, 53,54), "3(51 , 92, 53, 54), ul (51,52,33,54)] 

that satisfy 
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.^"•3(90,31,32, S3)«3(9o,93,fl , 4,S5)+n3(Si,92,ff3,94)n3(ffo,Si,94,95)+n3(fl , 2,S3,S4,S5)n3(9o,5l,92,S5)j /4 ^ 1 g 2 g 5 )x 

V 1 X (.9o, 92, 53, 54, 55)^4 (50, 5i ,33, 54, 95)^ 1 (go, 9i, 32,54,5s)^4 (50, 5i .92, 33, 55)^ (5o, 3i ,52, 33,54) = 1 (C8) 

I 



We will denote the collections of all d-cocycles 
Wd(j9o,—,9d),nd-i(gi,...,gd),u^_ 1 (jg 1 ,...,g d )] that sat- 
isfy the above conditions as iF d [G/, Ut(1)]- 

5. An "linearized" representation of group 
super-cohomology 

Since the coefficient of the group super-cohomology 
is always U(l), we can map the fermionic cochains 
[^d{9o,-,9d),n d -i{9i,-;9d),u g d ^ 1 {gi,...,g d )\ into two 
real vectors and an integer vector (v d , n<j-i, 
where the components of v d and h 9 d _ x are in [0, 1) 
and the components of n^-i are 0, 1. Then, we 
can rewrite the cochain and the cocycle conditions on 
Wd(9o,--,9d),n d ^ 1 (g 1 ,...,g d ),u 9 d ^ 1 (g 1 ,...,g d )} as "linear 
equations" on (v d , n d -i, ^-l)- 

To map v d (go, ■■■,g d ) into a real vector v d , we can in- 
troduce the log oi v d (g , g d ): 

Vd(go,9i,-,9d) = x-rlog^(5o,5i,---,5d) mod 1, (C9) 
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So we define v d as a |Gb| d+1 dimensional vector, whose 
components are given by 

(v d )( go ,... : g d ) = v d (l, go, g d ). (CIO) 
Similarly, n d -\ is a \Gt\ d dimensional integer vector 

(«d-l)(90,...,9i-i) = n d(5o, ■■■,5<i-i), (Cll) 
and h dl is a |Gfc| d dimensional vector 

( h d-i){9o,-,g d -i) =K-x(9o,-,9d-i), (C12) 
/i2-i(3o,-,3d-i) = YJ lnu d-i(9o, -,9d-i) mod 1. 

Now, the cochain conditions (ie the symmetry condi- 
tion), eqn. (C3), eqn. (CI), and eqn. (C4), can be written 
in a "linearized" form 



s(g)s d v d 


= v d + D d - X h 9 d _ x mod 1, 


S d -l n d-l 


= n d -i, 


D d -!n d -i 


= mod 2 


{9)SUK-x 


= K-v 




= h d U + s( 9l )hf_ x mod 1, 


h Pf 


= 2 n <J-i mod 1, 



3,31,32 GG/, (C13) 



where is a |Gf,| rf+1 x |Gb| d+1 dimensional matrix, 
whose elements are 

(C14) 

= ^(9 ,— ,9i),(99o,— ,99d+i) - 

Also -D^ is a |Gf,| <i+2 x |Gt| d+1 dimensional integer matrix, 
whose elements are all or ±1: 

( D d)(g ,...,g d+1 ),(g' ,...,g d ) ( C15 ) 
— (~) 3{g' ,...,g' d ),(go,---,gi,-..,gd+i)> 

where (50, 9i, g d +i) is the sequence obtained from 
the sequence (30, 3d+i) with the element gi removed. 
The matrices D d satisfy 

D d+1 D d = 0. (C16) 

A triple (v d , n d -i, h^-i) tnat satisfies eqn. (C13) corre- 
sponds to a fermionic d-cochain in < ^ 5d [G/, C/t(1)]- 

We can also rewrite the cocycle conditions eqn. (50) as 
a "linear" equation: 

D d v d + ^f d+1 = mod 1. (C17) 

Here f d is a |Gb| d+1 dimensional integer vector whose 
components are given by 

(fd)( go ,...,g d ) = fd(go,-,gd) mod 2. (C18) 

where f d {g , g d ) is obtained from n d _ 2 (5o, 5^2) 
according to the relation eqn. (48). A triple 
(v d ,n d -i,h d l ) that satisfies eqn. (Cf3) and eqn. (C17) 
corresponds to a fermionic d-cocycke in 3f d [Gf, £/r(l)]- 

6. Fermionic d-cohomology class 

As discussed in section V, each (d + l)-cocycle will 
define a fermionic topological non-linear a model in d 
spatial dimensions which is a fixed point of the RG flow. 
Those (d + l)-cocycles give rise to new fermionic SPT 
phases. However, different (d+ l)-cocycles may give rise 
to the same fermionic SPT phase. Those (d+ l)-cocycles 
are said to be equivalent. To find (or guess) the equiv- 
alence relation between cocycles, in this section, we will 
discuss some natural choices of possible equivalence rela- 
tions. A more physical discussion of this issue is given in 
Section G 1. 

First, we would like to show that (for d < 4) 

D d f d = mod 2, (C19) 
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where fd is determined from n<i-i(gi, gd) in 
Z d - 1 (G f ,Z 2 ) (see eqn. (48)). For 1 < d < 3, the 
above is trivially satisfied since the corresponding f d = 
0. To show the above for d = 4, we note that 
n d - 2 (g ,..,g d - 2 ) (which give rise to f d ) is a (d - 2)- 
cocycle in H d ~ 2 (Gb, Z 2 ). Let us assume that a 2- 
cocycle n 2 (g , gi, 52) gives rise to /4. The expression, 
74(30,51,52,53,54) = "-2(50, 5i, 52)^-2(52, 53, 54) implies 
that the 4-cochain 74(50,51,52,53,54) is the cup prod- 
uct of two 2-cocycles 712(50,51,52) and n 2 (g ,gi,g 2 ) in 
2 2 (G/,Z 2 ). So / 4 (5o, 5i,52,53, 54) is also a cocycle in 
Z 4 (G/,Z 2 ). Therefore D4/4 = mod 2. We can also 
show that -D5/5 = mod 2 by an explicit calculation 
using Mathematica. Therefore, for d < 5, 

A^ZViWd-i + \fd] = mod 1. (C20) 

for any fermionic (d— l)-cochains (v d -i, n d -i, 

There is another important relation. Let us assume 
that a (i-cocycle n d (go, gi, g d ) gives rise to / d+2 , and 
another d-cocycle 712(50,51, ■•• > 5d) gives rise to /d+2 (as- 
suming d < 3). If n d (g a ,g 1 ,...,g d ) and n d (g , g x , g d ) 
are related by a coboundary in B d (Gt, Z 2 ) 

«d(5o,5i, -,5d) +yj^d-i(gQ) -,5d) 

= nd(go,gi, -,gd) + (^d-i)(.9o, -,5d) 

= n d (g ,gi,-,gd) mod 2, (C21) 

then what is the relation between f d + 2 and /d+2? In the 
Appendix D, we will show that f d + 2 and /d+2 also differ 
by a coboundary D d+1 g d+1 in B d+2 (G b ,Z 2 ): 



d+2 



Dd+igd+i mod 2, 



(C22) 



where g d +i is determined from and m^-i (see Ap- 
pendix D). 

Using the above two relations, we can define the first 
equivalence relation between two fermionic d-cocycles. 
Starting from a fermionic (i-cocycle, (v d , n d +i, ^_i), we 
can use a fermionic (d— l)-cochain (v d -i, n d - 2 , h 9 d _ 2 ) to 
transform the above G?-cocycle to an equivalent d-cocycle 



Vd =v d - \^d + Dd-xVd-x, 



hd-i =n d - 1 + D d - 2 n d - 2 



D d - 2 h 3 d _ 2 , 



h 9 h 9 

n d-l —' l d-l 



(C23) 



where g d is determined from n d -i and n d _ 2 (see Ap- 
pendix D). We note that g d — if n d - 2 = 0. 

First, we would like to show that if (v d , n d _i, ft^-i) 
satisfies the cochain condition eqn. (C13), then 
(v d , n d+ i, ft?,) obtained above also satisfies the same 
cochain condition. Let us first consider the first equation 



in eqn. (C13): 

s (g)S d v d = s(g)S 9 v d - ^g d + s(g)S 9 D d _ 1 v d _ 1 

=v d + D d - 1 h 9 d _ 1 - ^g d + D d ^ 1 s(g)S 9 d _ 1 v d ^i mod 1 
=v d + D d ^ 1 (h 9 d _ 1 + D d ^ 2 h 9 d _ 2 ) mod 1, 
=ti d + D (i _i^_ 1 modl 1 (C24) 
where we have used 

S 9 d g d = g d , (C25) 

and the relation 

S 3 d D d -i = D d . x S 9 d _ x . (C26) 

It is easy to show that (v d , n d +i, fo^-i) a ^ so satisfies other 
equations in eqn. (C13). 

Now we like to show that if (v d , n d _i, h 9 d _A satisfies 

the cocycle condition eqn. (C17), then (^,^+1)^-1) 
obtained above also satisfies the same cocycle condition: 



D d v d + -fd+i 

= D d v d - ^D d g d + ^(fd+i + D d g d ) 
= mod 1. 



(C27) 



We can also define the second equivalence relation be- 
tween two fermionic d-cocycles. Let h d _ 1 (go, ...,g d -i) be 
representations of G b that satisfy 

h 9 d -i{go,gi, -)+s(g a )h 9 d b _ 1 (go,gi, ■•) = h 9 d ^{go,gi, ■•), 
s(g)hd-x(ggo, 55^-1) = ^-i(fo, -,gd-i), 
g,g a ,g b £G b . (C28) 

Then, starting from a fermionic d-cocycle, 
(vdirid+uh 9 ^), we can use such a /i^_ 1 (5o, ...,5<i-i) to 
transform the above G?-cocycle to an equivalent d-cocycle 
(v^hd-i,^^); 

"d(5o,5i, = Wd(5o,5i, —,9d) 

d 

+ ^-d-i(5i, -,9d) + ^2(-yhf_ 1 (go, -,9i, »,9d), 

i=l 

«d-i(5o,5i, —,9d-i) = nd-i(go,gi, -,gd-i), (C29) 
h d -i(go, -,gd-i) = ^_ 1 (5o,-,5rf-i) + h 9 d _ 1 (g ,-,gd-i)- 

The above can also be written in a more compact form: 

Vd = v d + Dd-ihd-i, 
ftd-i = n d-i, 

K-x=K-x + h g d-v (G30) 
where h d _x and h 9 d _ Y the \G b \ d dimensional vectors: 

(hd-\){ go ,..., gxi _ x ) = hf^igo^.-^d-i), 
(K-i)(go,-,g d -i) =K-i(S0,-,9d-i)- (C31) 
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Let us first show that (v d , fid-i, h d _-.) is a d-cochain: 
s(g)vd{ggo,9gi,-;99d) = s(g)v d {gg ,gg 1 ,...,gg d ) 

d 

+ s{g)h 9 d 9 J: x {ggi, ..,gg d ) + s(g) ^{-fh^iggo, -,99%, -, 

i=l 

d 

= v d (go,gi,-,gd) +^2(-TK-i(9°> ••»&,-, 9d) (C32) 

i=0 

+ K-\(9i, -.,g d ) + ^2(-yh g l\(g , -,k, -,5d) 

i=0 

= v d (go,gi,...,g d ) + ^(-TK_ 1 (g , 9d)- 

i=0 

Second, since = v d + Dd-ih d -i, it also clear that 
(v d , n d -i, ^-1) is a rf-cocycle. 

The equivalence classes of fermionic d-cocycles ob- 
tained from the above two kinds of equivalence rela- 
tions, eqn. (C24) and eqn. (C29), form the group d-super- 
cohomology J>? d [Gf , U T ( 1 )] • 

We note that, when n d -i = 0, the representations 
of Gf, h a d _ x , are also representations of So the 

second equivalence relation implies that a d-cocycle 
{v d ,rid-i = 0,h 9 dl ) is always equivalent to a simpler 
d-cocycle (vd, n d-i = O'^d-i = 0)- Such a kind of d- 
cocycles are described by group cohomology and form 
group cohomology class H d [Gb, Ut(1)}- 

Wc also note that, due to the requirement h d _ 1 = 
\n d -\ mod 1, for each fixed nd-i, all allowed choices of 
h 9 d _ 1 belong to one class under the second equivalence 
relation. So we just need to choose one allowed h 9 d _ x for 
each nd-i- 

7. Group super-cohomology classes and fermionic 
SPT phases 



non-linear er-model defined via the fermion path inte- 
gral in imaginary time (see section IV). The fermion 
path integral over a (d sp + l)-dimensional complex £ will 
give rise to a ground state on its d sp boundary. Such a 
)( gjound state represents the corresponding fermionic SPT 
state. So each element in group super-cohomology class 
J^' d [Gf, Ut(1)] corresponds to a fermionic SPT state. 

Note that, in d sp = spatial dimension, the ele- 
ments in Jf? 1 [Gf, Ut(1)] can have the trivial graded 
structure no(<7o) = 0, or the non-trivial graded struc- 
ture n (g ) = 1. The corresponding fermionic gapped 
states can have even or odd numbers of fermions. So we 
can have two different fermionic SPT phases in d sp = 
spatial dimension even without symmetry. 

8. A simplified "linear" equations for group 
super-cohomology 

The cochain condition (or the symmetry condition) 
eqn. (C13) can be solved, which implies that only the fol- 
lowing components v d {l,g 1 , ...,g d ), n d -i{l,gi,...,gd-i), 
and h 9 d _^{l,g\, ...,g d -i), are independent. Other compo- 
nents with go 7^ 1 can be expressed from those indepen- 
dent components: 

n d -i(go,gi,..gd-i) = nd-i{l,go 1 gi,-g ~ 1 9d-i), (C33) 
s{ga)h 9 d _ x {go,gi, ..,g d -i) = K-ii 1 ^ 9q X 9i, -, 9o l 9d-x), 

and 

Vd(go,gi, ■■•,3d) = s(g )(v d (l,g Q ~ 1 g 1 , ..^g^gd) 

d 

+ £(" W * (So -1 So, go 1 ^, ffo" 1 ^)) • (C34) 

i=0 



Each element in group super-cohomology class Using the independent components, we can rewrite the 
J%? d [Gf, Ut(l)] correspond to a fermionic topological co cycle condition as 



Vd(gi,92, -,9d+i) - Ud(l,P2,».,5d+i) + Vd0-,gi,-,gd+i) + fd+i0-,9i,92, -,9d+i) = mod 1, (C35) 

which can be rewritten as 

s(gi)(y d {l,g^ 1 g 2 , ■-,gT 1 9d+i) + s(g- 1 g 2 )h 9 d 1 _ 1 (l,g- 1 g 3 , ...,g~ 1 g d+1 ) - h 9 d 1 _ 1 (l,g- 1 g 3 , ...,g~ 1 g d+1 ) 

+ h 9 d -i( 1 '9i 1 92, ■■■,gi 1 gd+i) ) -Vd{l,g2,—,gd+i)+vd(i,gx,—,gd+i) + fd+iO-, 91,92, —,9d+i) 

= s(gi)v d (l,gi 1 g2, -, g^ga+i) ~ «d(l,32, -,9d+i) + v d {l,gi, -,gd+i) + 31,32, —,9d+i) 

-1 _ d+1 . -X 

+ s(92)h d l 1 (l,g 2 1 g 3 ,-,g 2 l 9d+i) - s{g x )'Y^{-) % h 9 d 1 _ l {l,g 1 1 g 2 ,...g 1 1 g l ,...,g 1 1 g d +i) = mod 1 (C36) 

i=2 
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These are the simplified "linear" equations for group 
super-cohomology. They include both the cochain (C13) 
and cocycle (C17) conditions. 

Appendix D: The mapping rid —> fd+2 induces a 
mapping H d (G b , Z 2 ) ->■ U d+2 (G b , Z 2 ) 

We note that eqn. (48) defines a mapping from a 
d-cochain rid € C d (G b ,Z 2 ) to a (d + 2)-cochain fd G 
C d+2 (Gb,Z 2 ). We can show that, if n d is a cocycle 
rid € Z d (G&,Z 2 ), then the corresponding fd+2 is also 
a cocycle f d £ Z d+2 {G b ,Z 2 ). Thus eqn. (48) defines a 
mapping from a ci-cocycle 6 Z d {G bl Z 2 ) to a (d + 2)- 
cocycle £ Z d+2 (G b , Z 2 ). In this section, we are go- 
ing to show that eqn. (48) actually define a mapping 
from d-cohomology classes to (d + 2)-cohomology classes: 
U d (G b ,Z 2 ) -> n d+2 (G b ,Z 2 ). This is because if n d and 
differ by a coboundary, then the corresponding fd+2 
and /d+ 2 also differ by a coboundary. 

To show this, let us first assume d = 2. Since the 
4-cochain / 4 (50,51, 52, 53, 9i) in C 4 (G b ,Z 2 ) is the cup 
product to two 2-cocycles n 2 (g , gi, g 2 ) and n 2 (g ,gi,g 2 ) 
in Z 2 (Gf,,Z 2 ), so if we change the two cocycles by 
a coboundary, the 4-cochain f 4 will also change by a 



coboundary. Thus there exist a vector 53 such that 

U = h + D 3 g 3 mod 2. (Dl) 

and eqn. (C22) is valid for d — 2. 

To show eqn. (C22) to be valid for d = 3, we need 
to show when n 3 (g , g 1} g 2 , g 3 ) and n 3 (g , gt, g 2 , g 3 ) are 
related by a coboundary in B 3 (G b , Z 2 ): 

"■3 (go,. 91, 92, 53) = "3(00, 51,52,53) + 013(51,52,53) 
+ m' 2 (g ,g 2 ,g 3 ) + m 2 (g , 51,53) + m 2 (g ,gi,g 2 ) mod 2, 

(D2) 

the corresponding / 5 and are also related by a 
coboundary in B 5 (G b , Z 2 ): 

f 5 = h + D 4 g 4 mod 2. (D3) 

Let us introduce 

«3(5o,5i,52,53) = ^(51,52,53) +m 2 (g Q ,g 2 ,g 3 ) 

+ m 2 (g 0) g 1 ,g 3 ) +m' 2 (g , g u g 2 ). (D4) 

Then we can express / 5 as f s =fs + f£ + f" + f" 1 , 
with: 



J5 (50,51,52,53,54,55) 

=n3(50, 51, 52, 53)n'3 (50, 53, 54, 55) + "3(5l,52,53,54)n3(50,5l,54,55) + "3(52,53,54,55)n3(50,5l,52,55), 

/5 /7 (5o,5i,52,53,54,55) 

="3(50, 51, 52, 53)^3(50, 53, 54, 55) + "3(51, 52, 53, 54)"3(5o, 5i, 54, 55) + "3(52, 53, 54, 55)"3(5o, 51, 52, 55), 
J5 77 (5o,5i,52,53,54,55) 

="3(50, 51, 52, 53)n' 3 (5o, 53, 54, 55) +"'3(51, 52, 53, 54)"3(5o, 5i, 54, 55) + n 3 (g 2 , g 3 , g^, g 5 )n 3 (g , g 1 , g 2 , g 5 ). (D5) 
In the following, we would like to show 54 can be constructed as 54 = 54 + 54 7 + 54 7/ , with: 

54 (5o , 5i , 52 , 53 , 54) = "3 (50 , 5i , 52 , 53 )"4 (50 , 53 , 54 ) + "3 (51 , 52 , 53 , 54 ) m' 2 (50 , 5i , 54 ) , (D6) 
54 J '(50, 5i, 52, 53, 54) = "3(50,52,53,54)^3(50,51,52) +"3(50, 5i, 52, 54)m 2 (52,53,54) + "3(5o,5i,53,.94)m' 2 (.gi,.g 2 ,5 3 ), 
54 /7 (5o,5i,52,53,54) = "3(50,51,52,53)^2(50,53,54) +"3(5i,52,53,54)m 2 (5o,5i,54) +"4(50, 5i, 52)^2(52, 53,54)- 

Let us prove the above statement in three steps: 

First step: After plug- in, we find that — -D454 i s given by 

(5o,5i,52,53,54,55) -5K5o,5i,52,53,54) -54(50,51,52,53,55) - 54(5o,5i,52,54,55) 
- 54(50,51,53,54,55) - 51(50,52,53,54,55) -54(51,52,53,54,55) 
= "3(50,51,52,53) [f"2(5o,53,54) +m' 2 (g ,g 3 ,g 5 ) + m' 2 (g a , g 4 , g 5 ) + m' 2 (g 3 , g 4 , g 5 )) 
+ "3(51,52,53,54) [m-s (50, 5i, 54) +m 2 (5o,5i, 55) + m 2 (5o,54,55) + "4 (51,54,55)] 
+ "3(52,53,54,55) [m 2 (50, 5i, 52) +m' 2 (5o,5i, 55) + m 2 (5o,52,55) + "4(51, 52, 5s)] 

- "3(5o,5i,52,53)"» 2 (5o,53,54) - "3(51, 52, 53, 54)"4(5o, 51, 54) - n 3 (g a , g l7 g 2 , g 3 )m' 2 (g , g 3 , g 5 ) 

- "3(51, 52, 53, 55)"4(5o> 5i,55) ~ "3(50,51,52,54)^2(50,54,55) - n 3 (g 1 ,g 2 ,g i ,g 5 )m 2 (g ,g 1 ,g 5 ) 
~ "3(50, 5i, 53, 54)r"2(5o,54,55) - "3(51, 53, 54, 55)f"2(5o,5i,55) ~ "3(50, 52, 53, 54)™ 2 (5o, 54, 55) 

- "3(52,53,54,55)^2(50,52,55) - "3(51,52,53,54)^2(51,54,55) - n 3 {g 2 ,g 3 ,g^g b )m' 2 {g 1 ,g 2 ,g b ) 
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We note that there are 5 terms containing m' 2 (go, 54, 55) and 5 terms containing m' 2 (go, <7ij 55)- Using the condition 
that "3(50,51,52,53) is a 3 cocycle: 

"3(50,51,52,53) + "3(50, 5i, 52, 53) + ^3(50,51,52,53) +^3(50,51,52,53) +"-3(50,51,52,53) = mode 2, (D7) 

we find that those terms cancel out. Also, there are two terms containing m' 2 (go, 53, 54), two terms containing 
m 2 (5o, 54, 55), e t c - Those terms also cancel out under the mod 2 calculation. But the m 2 (5o, 51,52) term and the 
77^2(53,54,55) term appear only once. Thus we simplify the above to 

H (5o,5i,52,53,54,55) - 54(50,51,52,53,54) ~ 54(50,51,52,53,55) - 54(5o,5i,52,54,55) 

- 54(5o,5i,53,54,55) ~ 54(50,52,53,54,55) - 54(51,52,53,54,55) 
="3(50,51,52,53)^2(53,54,55) + "3(52,53,54,55)^2(50,51,52) mod 2 (D8) 

Second step: Similarly, we have 

/5 7 (5o,5i,52,53,54,55) - 5^(50, 5i, 52, 53, 54) - 54 7 (5o,5i,52,53,55) - 54 7 (5o, 5i, 52, 54, 55) 

- 54 7 (5o,5i,53,54,55) - 5^(50, 52, 53, 54, 55) - 5^(51, 52, 53, 54, 5s) 

= - "3(50,51,52,53)^2(53,54,55) - "3(52,53,54,55)^2(50,51,52) mod 2 (D9) 

Third step: Finally, by using the results derived in the first step, we have 

fl H (go, 51, 52, 53, 54, 55) - gl H (go, 51, 52, 53, 54) - gl H {go, 51, 52, 53, 55) - 54 7/ (5o,5i,52,54,55) 

- 54 // (5o,5i,53,54,55) ~ 54^(50, 52, 53, 54, 55) - 54" (5i, 52, 53, 54, 55) 

="3(50, 5i, 52, 53)"T-' 2 (53, 54, 55) + "3(52,53,54,55)^2(50,51,52) 

+ "* 2 (5o, 5i, 52)f"2(52,53, 54) + m 2 (g , g l7 g 2 )m 2 (g 2 , g 3 , g 5 ) + m 2 (g , gx, g 2 )m 2 (g 2 , g 4 , g 5 ) 

+ "4 (5o , 5i , 53 ) f" 2 (53 , 54 , 55 ) + "4 (5o , 52 , 53 ) "4 (53 , 54 , 55 ) + "4 (5i , 52 , 53 ) "4 (53 , 54 , 55 ) mod 2 

=0 mod 2 (D10) 



Combine the results from the above three steps, we can 
show: 

^494 = // + /f + // W mod 2 (Df f ) 

Thus, we prove eqn. (C22) for c? = 3. 

We like to mention that the induced mapping 
H d (G b ,Z 2 ) -+ H d+2 {G b ,Z 2 ) by the n d -> / d+2 mapping 
(48), appears to be the Steenrod square 5g 2 . 64,65 This 
realization will allow us to generalize eqn. (48) to higher 
dimensions. 

Appendix E: Calculate T-L d (G b ,Z 2 ) - the graded 
structure 

1. Calculate H d (Z n ,Z 2 ) 

The cohomology group H (Z n ,M) has a very simple 
form. To describe the simple form in a more general 
setting, let us define Tate cohomology groups T-L d {G,M). 

For d to be or — f , we have 

n°(G,M) = M G /lmg(N G ,M), 
%- x {G,M) = Ker(N G ,M)/I G M. (Ef) 

Here M G , Img(iV G , M), Ker(7V G , M), and I G M are sub- 
module of M. M G is the maximal submodule that is 



r 



invariant under the group action. Let us define a map 
N G : M -> M as 

a^^g-a, a € M. (E2) 

gee 

lmg(N G , M) is the image of the map and Kei(N G , M) is 
the kernel of the map. The submodule I G M is given by 

I G M = n g9 ■ a\ £ n g = 0, a e M} (E3) 

In other words, I G M is generated by g • a - s, V g e 
G, a £ M. 

For d other then and — f, Tate cohomology groups 
H d (G,M) is given by 

H d (G, M) = H d (G, M), for d > 

•H d [G,M] = H- d -i[G,M], for d < -1. (E4) 

For cyclic group Z n , its (Tate) group cohomology over 
a generic Z n -module M is given by 51 ' 66 

U\Z n ,M) = {^r M \ if ^° m ° d2 ' (E5) 
v 7 \U- l {Z n ,M) ifd = lmod2. v ; 
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where 



H°(Z n ,M) = M z -/lmg(N Zn ,M), 
U-\Z n ,M) = Ker(N Zn ,M)/I Zn M. 



(E6) 



For example, when the group action is trivial, we have 
M z n _ M and I Zn M = Z x . The map 7V z?i becomes 
Nz n '■ a — > na. For M = Z2, we have lmg(N Zn , Z2) = 
nZ2 and Ker(iVz ji , Z2) = Z(„ 2 ) (where (n,m) is the 
greatest common divisor of n and to). So we have 



•H d (Z„,Z 2 



-2 

"(«,2) 



if d = 0, 
if d > 0. 



(E7) 



What is the non-trivial cocycle in - H d [Z 2 ,Z 2 ]? In 
Ref. 50, it is shown that a d-cocycle can be chosen to 
satisfy 

"dig, g,ff2,ff3, ■■■,gd-2,9d~i,gd) = 1 
^d(ffo,g, g,33, ■■■,9d~2,9d-i,9d) = 1 



vd(go,gi,g2,g3,-,gd-2,g,g) = h 



(E8) 



by adding a proper coboundary. 

Let us denote the two elements in Z 2 as {E, a} with 
a 2 = B. If we choose the cocycle Vd(go, 9i, 92, ■■■) in 
W d (Z2,^-2) to satisfy the above condition, then only 
Vd(E,a, E, ...) = Vd(cr, E, <r, ...) can be non zero. Since 
H d (Z 2 ,Z 2 ) = Z a , So v d (E,a,E,...) = u d (a,E,a, ...) = 
must correspond to the trivial class and v d {E, a, E, ...) = 
!/d(cr, E, er, ...) = 1 must correspond to the non-trivial 
class in U d {Z 2 , Z 2 ). 



2. Calculate H d (G,Z 2 ) from H d (G,Z) using Kiinneth 
formula for group cohomology 

We can also calculate H d {G,Z 2 ) from U d {G,J.) using 
the Kiinneth formula for group cohomology. Let M (resp. 
M') be an arbitrary G-module (resp. G'-module) over a 
principal ideal domain R. We also assume that either M 
or M' is i?-free. Then we have a Kiinneth formula for 
group cohomology 67 ' 68 



Zi -> II M ) ® fl U d - p {G\ M') 

-> "H d (G x G',M® R M') 

d+l 

^ ]jTorf(H p (G,M),H d - p+1 (G',A/')) ^Zi 

(E9) 



If both M and M' are B-free, then the sequence splits 
and we have 

H d {G x G',M® R M') 

d 

= j H?F(G,M)®rf- p (G',M') x 

p=0 

J n TOT f(^ P (G,M),H d - p+1 (G',M'))]. (E10) 

p=0 

If R is a field if, we have 

H d {G x G',M® K M') 

d 

= ]"[ [H p (G,M)® K H d - p (G',M')]- (Ell) 

p=0 

Let us choose R = Z, M = Z and M' = Z 2 . Note 
that M = Z is a free module over R = Z but M' = Z2 
is not a free module over R = Z. We also note that 
M ® z M' = Z« Z Z 2 = Z 2 . Let us choose G = Z x the 
trivial group with only identity. We have 



n d {z u z 2 



v -2 if d = 0, 
?! if d > 0. 



(E12) 



Thus we obtain 

Z a -> H d (G, Z) ® z Z 2 -> H d (G, Z 2 ) 
->Torf(ft d+1 (G,Z),Z 2 )->Zi. 



(E13) 



The above can be calculated using the following simple 
properties of the tensor product <S>z and Torf functor: 

1_®j_M = M <g) Z Z = M, 

Z„ ® z M = M «i z Z„ = M /nM, 

Z m ®Z Z n = Z( m n ), 

(A x B)(g> R M = (A M) x (B <g> fl M), 

M0_r (A x B) = (M ® R A) x (M® R B); (E14) 

Torf(Z, M) = Torf(M, Z) = Zi, 
Torf(Z n , M) = Torf (A/, Z„) = M/nM, 
Torf (Z m , Z„) = Z( m ,„), 

Torf (A x B, M) = Torf (A, Af) x Torf (B, M), 

Torf (M, A x B) = Torf (M, A) x Torf (Af, B). (E15) 

Here (to, n) is the greatest common divisor of to and n. 
Using 



n d (z n ,z) = { z r , 



if d = 0, 

if d = mod 2, 

if d = 1 mod 2, 



(E16) 
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we find that for d = even 

Zi -> Z„ ® z Z 2 -> H d (Z„, Z 2 ) -> Tnrf (Zi, Z 2 ) -> Z 1} 
or Zi -)■ Z (n ,a) ->■ H d (Z„, Z 2 ) Zi ->• Z x ; (E17) 
and for d =odd 

Zi -> Zi ® z Z 2 -> H d (Z„, Z 2 ) -> Torf (Z», Z 2 ) -> Z x , 
or Z a -> Zi -> ^ d (Z„, Z 2 ) -> Z ( „ )2) -> Zx. (E18) 

This allows us to obtain eqn. (E7) using a different ap- 
proach. 
Using 

' :, d = 0, 



H d [Z„ x Z 2 ,Z] < < 



Z 2 X Z n X Z (2 2 n) , d = mod 4, 

d — 1 



^(2 

Z 2 x Z 



(2,n)> 



d = 1 mod 2, 

d = 2 mod 4, 
(E19) 



we can show that for d = mod 4, d > 0, 

Zi -> (Z 2 x Z n x zff n) ) ®z Z 2 -> H d [Z„ x Z 2 , Z 2 
^Torf(z| !n) ,Z 2 )^Z l7 



or Zi — >• Z 2 x 



-(2,n) 



?T[Z„ x Z 2 ,Z 2 ] 



(E20) 

d 

(2,n) Zl ! 



for d = 1 mod 4, d > 0, 



1 -+ Z (2%) ®Z Z 2 ->■ W[^ n X Z 2 , Z 2 ] 



Torf(Z 2 xZ f 2 ,Z 



or £i -> 



"(2,'. 



(2,n)' 

H d [Z„ x Z 2 ,Z 2 ] Z 2 x 



(E21) 



for d = 2 mod 4. d > 0, 



Ii ->■ (Z 2 x Z ( 2 2 n) ) ® z Z 2 -► ft d [Z n x Z a , Z 2 ] 

^Torf(z| n) ,Z 2 )^Z l7 (E22) 



or Zi — > Z 2 x 



-(2,n) 



H a [Z n x Z 2 ,Z 3 ] ^ Z* 2tn) ^ Zn 



for d = 3 mod 4, d > 0, 

Zi -> zg } (g> z Z 2 -> ft d [Z„ x Z 2 , Z 2 ] 
-> Tor? (Z 2 x Z„ x Z ( ^ ); Z 2 ) -> Zj 



(E23) 



or £i -» 



?nz„ x z 2 



'-2 X 



-(2,n) 



Combining the above results, we find that, for d > 0, 

H rf [Z„xZ 2 ,Z 2 ] = Z 2 xZ d 2rl) . (E24) 

Using 



n d [u(i),z] 



Z if d = mod 2, 
Zi if d = 1 mod 2, 



(E25) 



we find that for d = even 

Zj -> Z ® z Z 2 -> H d [C/(l), Z 2 ] -> Traf (Zj, Z 2 ) -> Z x , 
or Zi -)• Z 2 -H d [C/(l), Z 2 ] ->• Zi Zi; (E26) 
and for d =odd 

Zi -> Zx ® z Z 2 -> H d [?7(l),Z 2 ] -> Toif (Z,Z 2 ) -> Zi, 
or Zj ->■ Zi H d (U(l),Z 2 ) Zi Zi. (E27) 
This allows us to obtain 



H d [U(l),Z 2 



l 2 if d = mod 2, 
?i if d = 1 mod 2, 



(E28) 



Note that in H d [U{l), Z] and H d [U(l),Z 2 }, the cocycles 
^ d : Z or : [U(l)} d+1 ->■ Z 2 are Borel 

measurable functions. 

Let us choose R = M = M' = Z 2 . Since R is a field 
and Z 2 (8>z Z 2 , we have 



H d (G x G', Z 2 ) = f[ [n p (G, Z 2 ) ®z 2 H d - p (G', Z 2 )' . 

(E29) 

Note that 



Z 2 0z 2 Z 2 = Z 2 , Z 2 ®z 2 Zi = Zi. 



Wc find 



H d [Z 2 x Z 2 ,Z 2 ] = Z% 



d+1 



(E30) 
(E31) 



Appendix F: Calculate the group super-cohomology 
classes J^ d [G f ,U T (l)] 

1. A general discussion 

From the discussions in Appendix C, we see 
that to calculate J$? d [Gf, J7t(1)] (which is formed 
by the equivalent classes of the fermionic cocycles 
[vd{9di— ,9d),nd-i(9a>— >ffd-i)])> we need to go through 
the following steps. 

First step: Calculate T~L d ~ 1 (Gb 1 Z 2 ) which give us differ- 
ent classes of dD-graded structure nd-i{go, — 5 ffd-i)- For 
each n d _ i (g , — , Sd-i ) , find an allowed x (gr , — , 9d- 1) ■ 
Second step: Next we need to find out B'H d ~ 1 (Gb,Z 2 ) 
from 'H d ~ 1 {Gb, Z 2 ). We note that for each graded struc- 
ture nd-i(go, •••) 9d-i) in W (Gfe,Z 2 ), the correspond- 
ing is a (d + 1) is a cochain in C d+1 [G&, f7y(l)]- 
Then, BH^ 1 {G b ,Z 2 ) denote the subset of H*- 1 ^, Z 2 ), 
such that for each graded structure rid-i(go, <7d-i) in 
B'H d ~ 1 (Gb,Z 2 ), the corresponding is a (d + 1)- 
coboundary in 2? d+1 [Gt,, £/r(l)]- Thus we have 

2-fd+i = AiS'd mod 1, 
for aJlnd-ito jd-dsBM^^Zj), (Fl) 
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where fd+i is obtained from rid-i(go, ffrf-i) through 
eqn. (48). We note that for n d -i(g , ffrf-i) ^ 
B'H d ~ 1 (Gb, Z 2 ), the corresponding will not be a 
coboundary in B d+1 [G b , U T (l)] and eqn. (C17) will not 
have any solution. 

Third step: We need to calculate Vd(go, gd) from 
eqn. (C17). For each fixed n^-iXgoi g d -i) € 
BH d - 1 {G bl Z 2 ), eqn. (C17) is reduced to 

D d (v d +g d ) = 0modl. (F2) 

The solutions v d of eqn. (F2) have a one-to-one corre- 
spondence with the cocycles in Z d [Gb, Ut(1)], provided 
that u^_i(fl , o, — ,9d-i) = 1 for 9 G G h . 
Fourth step: From the equivalence relation v d ~ v d + 
D d -iv d -i + \ fd, we see that the equivalent classes of 
the solutions can be labeled by H d [Gb, Ut(1)]- How- 
ever, the labeling is not one-to-one, since in obtain- 
ing H d [Gb, Ut(1)], we only used the equivalence rela- 
tion v d ~ v d + D d -\v d -\. The different classes in 
H d [Gb, Ut(1)] differ by \f d are regarded as equivalent 
under the group super-cohomology. The labeling is one- 
to-one only when \f d happen to be coboundaries in 

B d [G b ,U T (l)]. 

Through the above four steps, and combining with re- 
sults from the standard group cohomology, we can calcu- 
late the group super-cohomology classes J4 ad {Gf, Ut(1)]- 
We like to point out that a subset of J$? d [Gf, J7t(1)] that 
has trivial e?D-graded structure n d _ i(g , g d _ i ) = is 
given by H d [G f , U T {1)]. Therefore, the fermionic SPT 
phases contain the bosonic SPT phases with symmetry 
Gb as special cases. In the following, we will discuss some 
simple examples. 

2. Calculate ,ff d \Z x x z£, 17(1)] 

Here we choose Gb = Z\ (the trivial group) which cor- 
respond to fermion systems with no symmetry. Since 
H d - 1 {Z 1 ,Z 2 ) = Z 2 for d = 1 and n d ~\Z x ,T % ) = T x = 
for d > 1. We find that there are two ID-graded 
structures and only one trivial dD-graded structure for 
d > 1. For each <iD-graded structure, we can choose a 
u d -i(9o,~ ,9d-i) that satisfies u^^go, —,g d -i) = 1 for 
g G G b - 

For d > 1 and for the only trivial graded structure, the 
equivalent classes of v d are described by H d [Zi,U(l)] = 
Z x . Therefore JT d [Z l x z£, £7(1)] = Z x for d > 1. 

For d = 1, H°(Z 1 ,Z 2 ) = Z 2 . The trivial graded struc- 
ture is given by n^(E) = and the other non-trivial 
graded structure is given by nW(£) = 1. The corre- 
sponding f 2 a \g Q , ...,5d) are given by / 2 (a) = 0, a = 0, 1. 
Thus for the each graded structure, the equivalent classes 
of v d are described by U(l)] — ~L\. Therefore 

Jf 1 ^ x Z£,U(1)] = Z 2 . To summarize 

^ 1 x4,C/(l)] = |^ ^ (F3) 



x Zl,U(l)} = Z 2 implies that there are two 
possible fermionic SPT phases in d sp = spatial dimen- 
sion if there is no symmetry. One has even number of 
fermions and the other has odd number of fermions. 

Jff 2 \Z\ x U(l)] = Zi implies that there is only one 
trivial gapped fermionic SPT phases in d sp = 1 spatial 
dimension if there is no symmetry. It is well known that 
there are two gapped fermionic phases in d sp = 1 spatial 
dimension if there is no symmetry, the trivial phase and 
the phase with boundary Majorana zero mode. 5 One 
way to see this result is to use the Jordan- Wigner trans- 
formation to map the ID fermion system with no sym- 
metry to an ID boson system with Z 2 symmetry. ID bo- 
son system with Z 2 symmetry can have only two gapped 
phase: the symmetry unbroken one (which corresponds 
to the trivial fermionic phase) and the symmetry break- 
ing one (which corresponds to the Majorana chain). 48 
Both gapped phases can be realized by non-interacting 
fermions (see Table II). However, the non-trivial gapped 
phase (the bosonic Z 2 symmetry breaking phase) has 
non-trivial intrinsic topological orders. So the non-trivial 
gapped phase is not a fermionic SPT phase. This is why 
we only has one trivial fermionic SPT phase in ID if there 
is no symmetry. 

Jf 3 [Zi x z£,U(l)] = Zi implies that there is only 
one trivial gapped fermionic SPT phases in d sp = 2 spa- 
tial dimension if there is no symmetry. It is well known 
that even non-interacting fermions have infinite gapped 
phases labeled by Z in d sp — 2 spatial dimensions if there 
is no symmetry (see table II). Those phases correspond 
to integer quantum Hall states and/or p + ip, d + id 
superconductors. Again, all those gapped phases have 
non-trivial intrinsic topological orders. So we only have 
one trivial fermionic SPT phase in d sp = 2 spatial dimen- 
sions if there is no symmetry. 

3. Calculate ,7f d [Z 2 k+i x Z(,U(i)] 

Next we choose G/ = Z 2 u+i X Z|. Again, since 
U d - 1 {Z 2k+1 ,T 2 ) = Z 2 for d = 1 and H d - 1 (Z 2k+1 , Z 2 ) - 
Zi for d > 1, there are two ID-graded structures and 
only one trivial dD-graded structure for d > 1. For each 
c?D-graded structure, we can choose a u^_ x (go, g d -i) 
that satisfies u^ l _ 1 (go, ...,g d -i) = 1 for g € Gb- 

For d > 1 and for the only trivial graded struc- 
ture, the equivalent classes of v d are described by 
H d [Z 2k+1 ,U(l)} = Zx for d=cvcn&ndn d [Z 2k+1 ,U(l)} = 
Z 2k+1 for d=odd. Therefore Jf d [Z 1 x z£, U(l)] = Z x for 
d=even and x Z f 2 , U(l)] = Z 2k+1 for d=odd. 

For d = 1, H°(Z 2k+1 ,Z 2 ) = Z 2 . The trivial graded 
structure is given by n^(go) = and the other non- 
trivial graded structure is given by n^'(go) = I- The 
corresponding f 2 a \g , g d ) are given by f { 2 a) = 0, a = 
0, 1. Thus for the each graded structure, the equivalent 
classes of v d are described by H 1 [Z 2k+ i, U(l)] = Z 2fe+1 . 
Therefore ^ x \Z 2h+x x z{, 17(1)] = Z 2 x Z 2k+1 = Z ik+2 . 
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To summarize 



jf d [z 2k+1 x zi,u(i)] = \ z 



i/, 2- d — 1, 
i . d = even, d > 0, 
d = odd. 

(F4) 



-2k+l, 



For d = 3, H 2 (Z 2 ,Z 2 ) = Z 2 . For the non-trivial 3D- 
graded structure, n 2 (e, o, e) = 722(0, e, o) = 1 and others 
ri 2 = 0, the corresponding / 4 has a form 

fi{e, o, e, o, e) = ^4(0, e, o, e, o) = 1, others = 0. (F6) 



<%> X \Z<ih+\ x 4, = Z 2 A;+i x Z 2 implies that there 

are 4fc + 2 possible fermionic SPT phases in d sp = 
spatial dimension if there is a Z 2k+ \ symmetry. 2k + 1 of 
them have even number of fermions and the other 2k+l of 
them have odd number of fermions. The 2k+l phases are 
separated by the 2fc+l possible Z 2k+ i quantum numbers. 



is a cocycle in 

Z lt (— )/*(ffO,—,S4) 



One can show that (-)Mffo,-,ff4) 
Z 4 [Z 2 ,77(1)]. Since H 4 [Z 2 , 77(1)] = 
is also a coboundary in B 4 [Z 2 ,U(1)]. Therefore, 
BH 2 (Z 2 ,Z 2 ) = Z 2 . 

For d = 4, -H 3 (Z 2 ,Z 2 ) = Z 2 . For the non-trivial 4D- 
graded structure, ri3(e, o, e,o) = 713(0, e, o, e) = 1 and 
others n 3 = 0, the corresponding / 5 has a form 



^ 2 [Z 2k+l x z£,U(l)] = Zi implies that there is no / 5 (e, o, e, o, e, o) = / 5 (o, e, o, e, o, e) = 1, 



non-trivial gapped fermionic SPT phases in d sp = 1 spa- 
tial dimension if there is a Z 2k+1 symmetry. We can use 
the Jordan- Wigner transformation to map ID fermion 
systems with Z 2k+1 symmetry to ID boson systems with 



Z 



f _ 



2k+ i x2 2 '= Z 4k+2 symmetry. Since H 2 [Z Ak+2l [7(1)] = 
Zi, ID boson systems have only two phases that do 
not break the Z 2k+ i symmetry: one does not break the 
Z 2 symmetry and one breaks the Z 2 symmetry. The 
fermionic SPT phase has no intrinsic topological order 
and correspond to the bosonic phase that does not break 
the Z^ symmetry (and the Z 2k+ \ symmetry). Thus ID 
fermion systems with Z 2k+ \ symmetry indeed have only 
one trivial fermionic SPT phase. 

Jf z [Z 2k+ i x zl,U(l)] = Z 2/ r. +1 implies that our con- 
struction gives rise to 2fc+l gapped fermionic SPT phases 
in d sp — 2 spatial dimension if there is a Z 2k+ i symme- 
try. The constructed 2fc+l gapped fermionic SPT phases 
actually correspond to the 2k + 1 gapped bosonic SPT 
phases with Z 2k+ i symmetry [since H 2 (Z 2k+ i, Z 2 ) = Zi 
and there is no non-trivial graded structure]. 



others = 0. 

(F7) 



One can show that (— )/s(flo>"-,95) } s a cocycle in 
Z 5 [Z 2l U{\)\. Since H 5 {Z 2 ,U{1)} = Z 2 and (-)/st?o,-,fls) 
corresponds to a non-trivial cocycle in H 5 [Z 2 , U(l)] — 
Z 2 , we find that RH 3 {Z 2 ,Z 2 ) = Z x . 

Following the third step in section Fl, we can show 
that the elements in j4? d {Z 2 x zl,U(l)} can be labeled 
by H d [Z 2 x Z(, U(l)] x BH d - 1 (Z 2 , Z 2 ). 

Following the fourth step in section F 1, we like to show 
that the above labeling is one-to-one. For d = 1,2,3, 
fd(go, 9d) = 0. Thus the labeling is one-to-one. For 
d = 4, a 3D-graded structure in BH 2 (Z 2 ,Z 2 ) gives rise 
to a fi{go,.—,gA). Since (— )^^a, —,34) j s a cocycle in 
Z 4 [Z 2} U(1)} and since H 4 {Z 2 , U(l)] = Z u we find that 
(_)/4(»o.-.fl4) i s a i so a coboundary in B A [Z 2 ,U{1)] and 
the labeling is one-to-one. 

So, using 



n d [z n ,u(i)} 



d = mod 2, d > 0, 
d = 1 mod 2, 



(F8) 



4. Calculate Jt? d [Z 2 x Z{, [7(1)] 



Now, let us choose G/ = Z 2 x Z^. Following the 
first step in section F 1, we find that / H d ~ 1 (Z 2 ,Z 2 ) = 
Z 2 and there are two graded structures in all dimen- 
sions. For each dD-graded structure, we can choose a 
Ud-iiao,— >9d-i) that satisfies u^igo, — ,fld-i) = 1 for 
9 G G b . 

Following the second step in section Fl, we 
note that the trivial graded structure is given by 
rid_i(<70) •••) gd-i) = 0. The non-trivial graded structure 
is given by 



we find that 



o() 



n<j_i(e,o,e, ...) = n d -i(o, e,o, ...) = 1, 



others = 0. 

(F5) 



where e represents the identity element in Z 2 and o rep- 
resents the other element in Z 2 . 

For d = 1,2, U d - 1 {Z 2 ,Z 2 ) = Z 2 . For each 
graded structure rid-i(go, ...,gd-i), the correspond 
f d+1 {g ,...,g d+1 ) = 0. Thus M^H^,^) = Z 2 . 



J? d [Z 2 x Z|,C7(1)] 



•2 

2: 


d = 


1, 


2) 


d = 


2, 


2 
2' 


d = 


3, 


1- 


d = 


4. 



(F9) 



Jf 1 ^ x Z/f, 77(1)] = Z\ implies that there are 4 pos- 
sible fermionic SPT phases in d sp = spatial dimension 
if there is a Z 2 symmetry. 2 of them have even number 
of fermions and the other 2 of them have odd number of 
fermions. The 2 phases with even number of fermions are 
separated by the 2 possible Z 2 quantum numbers. The 
2 phases with odd number of fermions are also separated 
by the 2 possible Z 2 quantum numbers. 

J$? 2 {Z 2 x Z 2 ,U(1)] = Z 2 implies that there are two 
gapped fermionic SPT phases in d sp — 1 spatial dimen- 
sion if there is a Z 2 symmetry. We can use the Jordan- 
Wigner transformation to map ID fermion systems with 
Z 2 symmetry to ID boson systems with Z 2 x z£ sym- 
metry. Since H 2 [Z 2 x Z 2 , 77(1)] = Z 2 , ID boson systems 



:-!S 



have two phases that do not break the Z 2 x Z 2 symme- 
try. ID boson systems also have another phase that do 
not break the Z 2 symmetry: the phase that break the Z| 
symmetry. The fermionic SPT phases have no intrinsic 
topological order and correspond to the bosonic phases 
that does not break the Z| symmetry (and the Z 2 sym- 
metry). Thus ID fermion systems with Z 2 symmetry 
indeed have two fermionic SPT phases. 

^ 3 [Z 2 x Z|, U(l)] — 2-2 implies that our construction 
gives rise to 4 gapped fermionic SPT phases in d sp = 2 
spatial dimension if there is a Z2 symmetry. 

5. Calculate J^ d [Z^ x Z S 2 ,U T {1)} 

Now, let us choose Gf = Zj x Z%. Since 
H d ~ 1 (Z 2 r , Z 2 ) = 2.2, there are two graded structures in 
all dimensions. The trivial graded structure is given by 
nrf-i(s'O) -••> 9d— 1) = 0- The non-trivial graded structure 
is given by 

ra d _i(e, o,e, ...) = n d -i(o, e,o, ...) = 1, others = 0, 

(F10) 

where e is the identity element in Zj and o is the 
non-trivial element in Zj. For each dD-graded struc- 
ture, we can choose a u^_ 1 (go, gd-i) that satisfies 
Ud-iCffo, —,9d-i) = 1 for g £ G b . 

For d = 1,2, H d - 1 (Z 2 r ,2 2 ) = Z 2 . For each 
graded structure Tid-i(go, 9d-i)> the corresponding 
fd+i(d0t •••) 9d+i) — 0; an d the equivalent classes of Vd 
are described by H d [Zj,U T {l)], Therefore Jf^Zj x 
Zl, U T {1)] = Z 2 and 2 [Zj x Z|, U T (1)} = 2\. 

For d = 3, H 2 (Z 2 r ,Z 2 ) = Z 2 . For the non-trivial 3D- 
graded structure, n 2 (e, o, e) = n 2 (o, e, o) = 1 and others 
n 2 =0, the corresponding / 4 has a form 

/4(e, o, e, o, e) = /^o, e, o, e, o) = 1, others = 0. 

(Fll) 

One can show that (— ~)f*(Bo>—,gt) j s a cocycle in 
Z 4 [Z 2 T ,[/ T (1)]. Since H^Z? , U T (l)} = Z 2 and 
(_)/4(»o,...,fl4) i s a non-trivial cocycle in H^Zj, U T (1)], so 
for the non-trivial 3D-graded structure, eqn. (C17) for v d 
has no solutions. For the trivial 3D-graded structure, the 
equivalent classes of v d are described by 'H 3 [Z 2 r , f7<r(l)]. 
Thus JT^Zj x Zl, U T {1)] = Zi. 

For d = 4, n 3 {Z 2 r ,2 2 ) = Z 2 . For the non-trivial 4D- 
graded structure, 713 (e, o, e,o) = 713(0, e,o,e) = 1 and 
others 77,3 = 0, the corresponding /g has a form 

/ 5 (e, o, e, o, e, o) = / 5 (o, e, o, e, o, e) = 1, others = 0. 

(F12) 

One can show that (— ^/sCsoviSb) j s a cocycle in 
2 5 [Z 2 T ,C/ T (1)]. Since H 5 [Z 2 T ,f/ T (l)] = Z x , (-)ACa.,..-,fl.) 
is also a coboundary in B 5 [Zj, t/V(l)]. So for the 
non-trivial 4D-graded structure, the equivalent classes 



of v d are labeled by U A [Z^ ,U T {1)] = Z 2 . Since 
(_)/4(ffo,— ,34) i s rea i anc j can kg viewed as an el- 
ements in "H 4 [Z 2 ,[/(1)]. Since H 4 [Z 2 ,C/(1)] = Z x , so 
(_)/4(so,---,94) i s a i wa ys a coboundary in B 4 [Zf, E7r(l)], 
the labeling is one-to-one. For the trivial 4D-graded 
structure, the equivalent classes of Vd are labeled by 
•H 4 [Z 2 r ,C/ T (l)] = Zi. Since (-JAGjo.-.m) is a i way s a 
coboundary in fr^Z-f, [/^(l)], the labeling is one-to-one. 
Thus jr 4 [Zj x Z| , f7 T (l)] = 2\. To summarize 

{Z 2 , d=l, 
J T=\ (F13) 
Zi, d = 4. 

'H 1 [Z 2 r x zI,Ut{1)] = Z 2 implies that there are two 
possible fermionic SPT phases in d sp = spatial dimen- 
sion if there is a Zj symmetry. One has even number of 
fermions and the other one has odd number of fermions. 

■H 2 [Zj x Z 2 f ,C/ T (l)] = 2\ implies that there is four 
gapped fermionic SPT phases in d sp = 1 spatial dimen- 
sion if there is a Zj symmetry. We can use the Jordan- 
Wigner transformation to map ID fermion systems with 
Z2 symmetry to ID boson systems with Zj x Z| sym- 
metry. Since 'H 2 [Z 2 r x Z 2 ,U T (1)} = Z\, ID boson sys- 
tems have four phases that do not break the Zj x z! 2 
symmetry ID boson systems also have four phases that 
do break the Zl symmetry: two do not break the time- 
reversal T symmetry and two do not break the TPt sym- 
metry. The fermionic SPT phases have no intrinsic topo- 
logical order and correspond to the bosonic phases that 
do not break the Zl symmetry (and the Zj symmetry). 
Thus ID fermion systems with Zj symmetry indeed have 
four fermionic SPT phases. 

'H 3 [Z 2 r x zI,Ut(1)] = Zi implies that our construc- 
tion only gives rise to one trivial gapped fermionic SPT 
phases in d sp = 2 spatial dimension if there is a Z 2 sym- 
metry. H 4 [Z 2 r x zI,Ut(1)] = Z| implies that our con- 
struction gives rise to three non-trivial gapped fermionic 
SPT phases in d sp = 3 spatial dimension if there is a Zj 
symmetry One of the non-trivial SPT phases is actially 
the purely bosonic SPT phase. The other two are intrin- 
sically fermionic. We note that, using non-interacting 
fermions, we cannot construct any non-trivial gapped 
phases with Zj symmetry (the T 2 = 1 time-reversal sym- 
metry) . 

Appendix G: Properties of the constructed SPT 
states 

1. Equivalence between ground state wave 
functions 

In Section VI, we have demonstrated how to construct 
an ideal ground state wave function $({<7i}, {ny...^}) in 
d spatial dimension from a (d+ l)-cocycle (vd+i, n d , u 9 d ). 



39 



However, a natural question is: how do we know these 
wave functions describe the same fermionic SPT phases 
or not? Here we will address this important question 
based on the new mathematical structure - group super- 
cohomology class invented in Appendix C. 



a. A special equivalent relation 

Let us choose arbitrary {yd, Jid-i, it ji) in one lower di- 
mension. Here (y d , n d -i, u d _ x ) * s n °t a <i-cocycle. They 
just satisfy the symmetry condition eqn. (46). In Ap- 
pendix C, we have shown that {vd+i,hd,u 9 d ) obtained 
from {vd+i,nd,u 9 d ) through eqn. (C23) is also a {d + 1)- 
cocycle. Here we will assume that rid-i — and u 9 d _ 1 = 
1, and the resulting 

Vd+i = v d +i5v d , n d = n d , u 9 d = u% (Gl) 

is a (d+l)-cocycle. The new (d+l)-cocycle (Od+i,nd> u 9 d ) 
gives rise to a new ideal wave function $'({<?i}, {nij...k\)- 
Should we view the two wave functions $({<7i}, {nij...k}) 
and $'({<7j}, {n^...*,}) as wave functions for different SPT 
phases? 

From the definition of 8fd (49), we can show that the 
two wave functions $({3i}> {n^}) and $'({#;}, {n ijk }) 
are related through 

{nijk}) = $({.9*}, {n ijk })x 
Y[ (-)^(»'«'*-"-'») I /*«"-*( ffi , fc -,... sSfc ), (G2) 

(ii-fc) 

where Yluj fc) * s t ne product over all the ei-simplices that 
form the d dimensional space [in our (2+l)D example, 
Yltij k) IS t ne product over all triangles]. Also Sjj.„fc = 
±1 depending on the orientation of the d-simplex {ij...k). 
Since v s d 3 " k {gi, gj, ■■■,gk) is a pure U{1) phase that satis- 
fies v d iS ' k (ggi,ggj,.-,ggk) = v d 3 k {g t ,gj,...,g k ), it rep- 
resents a bosonic symmetric local unitary (LU) transfor- 
mation. Therefore the new wave function <&' induced by 
{v d 7^ 0, n d -i = 0,u d _ 1 = 1) belong to the same SPT 
phase as the original wave function <I>. 

Also since fd+i for d < 3 is always zero. Therefore, 
below 3 spatial dimensions, new wave function belongs 
to the same SPT phase as the original wave function 

In 3 spatial dimensions and above {d > 3) , if we choose 
{vd = 0,n<i_x 7^ 0), the new induced wave function $' 
will differ from the original wave function $ by a phase 
factor Y[( k ^{—)f<i+^9o,gi,gj,---,gk)_ (_)/d+i(so,9«,s.j,-",Sk) 

represents a LU transformation, and the two wave func- 
tions, <&' and $, have the same intrinsic topological order. 
But since {—\f<i+^ao,ggi,99j,—,99h) ^ r_)/d+i(90ifl«.9j ) ---,9*,) ) 

{—y<i+i(go,gi,9j,'--,9k) does not represent a bosonic sym- 
metric LU transformation. So we do not know whether 
$' and $ belong to the same fermionic SPT phase or not. 

We would like to point out that in order for $' and 
$ to belong to the same fermionic SPT phase, the two 



wave functions can only differ by a fermionic symmet- 
ric LU transformation (which is defined in Ref. 25). 
The bosonic symmetric LU transformations are a sub- 
set of fermionic symmetric LU transformations. So even 
though (— yd+i(go,9i,9j,---,9k) j s n0 ^ a bosonic symmet- 
ric LU transformation, we do not know whether it is 
a fermionic symmetric LU transformation or not. The 
structure of group super-cohomology theory developed in 
Appendix C shows that /<j+i(<7o, g%, gj, <?&) is a d+1- 
cocycle with Z2 coefficient, suggesting the sign factor 
(_)/<i+i(ffo,9».9.}V",sj ! ) can be generated by fermionic sym- 
metric LU transformation. 



b. Generic equivalent relations 

In the above, we discussed a special case of 
{vd+i, nd, u 9 d ) where hd and u d are the same as nd and u 9 d . 
However, the equivalent class of group super-cohomology 
has a more complicated structure. In Appendix C, we 
have shown that {vd+i, nd, u 9 d ) obtained from: 

Vd+i = v d +i{-) 9d , n d = n d + Sm'd-i, u B d = u 9 d , (G3) 

(with g d {g , ••• ,3d) = (5d)( So ,... , 9d ) mod 2 ) also belongs 
to the same equivalent class of group super-cohomology. 
We note that here the graded structure changes into 
hd- It would be much harder to understand why the 
corresponding new wave function still describes the same 
fermionic SPT phase. Let us consider the 2 + ID wave 
function on a sphere as a simple example. Since the 
wave function derived from the fixed point amplitude 
Eq. (30) is a fixed point wave function with zero cor- 
relation length, it is enough to only consider a mini- 
mal wave function with 4 sites on a sphere(labeled as 
1,2,3,4). Such a minimal wave functions will contain 4 
triangles 123, 124, 134 and 234. From the definition of Vd, 
we can show that the two wave functions §{{g^\, {nijk}) 
and <&{{gi}, {hijk}) are related through {n' ijk }) = 

(_)94(si,92,s3,S4)$({ 5 .} { niifc })_ To see this explicitly, let 
us construct the ground state wave function by adding 
one more point inside the sphere. In this way, the wave 
function can be formally expressed as: 

*(<7l> 32, 33, 34, #234, #124, #134, #123) 

={-) mi(9u93) J V 8 - (30,31, 32,33)V+(3o,3i, 32,34) 

x V3" (30 , 31 , 33 , 34 ) V3 1 " (30 , 32 , 33 , 34 ) 
=(-r i(9l ' 93) V 3 + (3i,32,33,34) 

="3 (91 ,32, 33, 34) X 

/)»2 (92 ,93 ,94 ) fl «2 (91 ,92 ,94) /)"2 (92 ,93 ,94) /)" 2 (gi ,92 ,93 ) / ri a \ 
W 234 °124 y 234 °123 ■ l" 4 J 

Again, the symbol J is defined as integrating over the 
internal Grassmann variables with respecting to the 
weights toi on the internal links(see Eq.(18)). Similarly, 
the wave function corresponding to a different solution 
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(vd+i, fid, u 9 d ) takes a form: 

* (5l , 52 , 93 , 54 , #234 , #124 , #134 , #123 ) 

=^(5i,52,33,34)(-) 94(9l ' 52 ' ! ' 3 ^ ) x 

fj n 2 (92 ,33 ,94)n n 2 (9i ,92 ,34) 7jn' 2 (92 ,93 ,94) 7j n 2 (fli >fla>fls) /pi r\ 
"234 (7 124 p 234 t7 123 • K^ ) 

We note that 34(331,332,333,334) = 34(51,52,53,34) 
implies the phase factor (—y4(gi, 92,93,94) j s a S y mm et- 



with n 2 (g l ,g J ,gk) = n 2 (g l ,g j ,g k ) + m[(gi,gj) + 
TO 'i(5j,5fc) + m 'i(5fe,3i)- We note that both n 2 and 
m'j are invariant under symmetry transformation. Here 

( c fefe)) is defined as (c\ ijk) J = c (ijk) . ~ means 
equivalent up to some symmetric sign factors which will 
arise when we reorder the ferminon creation/anihilation 
operators in U'. We also note that U' is a generalized 
LU transformation which only has non-zero action on 
the subspace with fixed fermion occupation pattern n 2 . 
The above discussions can be generalized into any di- 
mension. Thus we have shown the wave function con- 
structed from the above (vd+i, ?V, uf) describes the same 
fermionic SPT phase. 

Finally, we can consider more general case where u 9 d 
also changes into u d . It is very easy to see that such 
changes in the wave function can be realized by symmet- 
ric LU transformation since those u d evolve 50 are all 
canceled in the wave function. 

In conclusion, we have shown that (vd+i,rid,u 9 d ) 
and (ud+i,nd,u d ) belonging to the same group super- 
cohomology class will give rise to fixed point wave func- 
tions describing the same fermionic SPT phase. 

2. Symmetry transformations and generic SPT 
states 

In the Section VI B, we have constructed an ideal 
fermion SPT state |\?) = C/|$o) using the the fermionic 



ric LU transformation. However, we still do not know 
whether the two wave functions with different patterns 
of fermion number {fiijk} and {n^} describe the same 
fermionic SPT phase or not. 

Indeed, we find that up to some symmetric phase fac- 
tors, the state |$) and |$) can be related through a sym- 
metric LU transformation |<E>) ~ U'\$), where 



(G6) 



I 

LU transformation U constructed from a cocycle (1/3, n 2 ). 
In this section, we are going to discuss how to construct 
a more generic SPT state \$f') that is in the same phase 
as the ideal state \^). 

Clearly, and \*S?) has the same symmetry. So to 
construct \^'), we need first discuss how symmetry trans- 
formation changes under the fermionic LU transforma- 
tion U. 

The symmetry G}, that act on |^) is represented by 
the following (anti) unitary operators 



W(g) ^^2Y[[u 9 2 {gi,g 3 ,g k )} 1 Y[i4(gi,9j,9k)x 

{<?*} V A 

ni5,)(55d^ 1 ^ M , (G7) 



where Wi is the symmetry transformation acting on a 
single site i, and K is the anti unitary operator 



KiR- 1 = -i, KcR- 1 = c, Kc^R- 1 = c f . (G8) 



We find that 



r m'i(gi,9j)-m' 1 (gi,g j ) , 
"(Oij) 



(Oij) 



n 

A 

n 



Jm'i (gi ,gj ) t"*i (gj ,g h ) -fmi (g< ,gu) 



'(Oij) 



-{Ojk) 



-(Oik) 



( C fefc)) 



"2 (gi ,9j ,9k ) ~n 2 (gi ,9j ,9k ) 



t \ "2 ta ,9 3 ,g k ) - n a ( 9i , gj , n ) Jm[ {g . in j fm , (g . >gh j fm , {g . >g . } 



'(ijfc) 



■'(Oife) 



'(Ojfc) "(Oij) 
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W° (g) = t/t W(g) U = £ II \9i) (99i I II ^ (ft , ft , 3* )] " 1 II u 2 (ft , ft , ft ) 

{a,} * 

IT 



(G9) 



V A 

"2 (Si ,Sj ,flfc ) n 2 (g ,9i,9k) -™2 (go ,flj ,9k ) -«2 (90>ff» ,flj) TT "2 (go ,9i ,9j ) ™2 (go ,9j ,9fc ) -™2 (go ,9i ,9k ) -™2 (Si ,Sj ,flk ) , 



(ijfc) 



'(Oik) 



'(Ojfc) 



'(Oij) 



n 

A 



'(Oij) 



'(Ojfe) 



-'(Oik) 



-(ijk) 



n ^(so, ft, ftnfton ^ x {9o,9i,9j,9k) n 



„t™2 (go, 9;, 9j)-t™2 (go, 9;, 9i) . 



'(Oij) 



"(Oij) 



A 



(Oij) 



n "(ou) 

(Oij) 

n 



"2 (90 ,99i .99j ) -ni (go ,99i ,99j ) 



'(Oij) 



II [^3 1 (ft , 3ft , 33j . 99k )] s{s) II [^3 (ft , 3ft , 3ft :99k)] s(9) K- 



2 x 



'(Oij) 



'(Ojfc) 



t"2 (go ,ggi ,ggj ) t«2 (go ,gg 3 ,ggfc ) -fn 2 (go ,gg; ,ggfc ) -t«2 (gg; ,ggj ,ggjc ) TT t«2 (gg< ,ggj ,ggi<, ) Jn 2 (g ,ggi ,ggk ) -t™2 (go ,ggj ,ggk ) -t«2 (go ,gg» ,ggj ) 

Wn„'„-'\ Wri-'!^ Wn„l,\ W-.»\ I W. 



(Oifc) 



(ijfc) 



'(Oifc) 



'(Ojft) 



'(Oij) 



A 



Using the fact that C/ is independent of go , we can change 
the go in the second half of the right-hand-side of the 
above expression to gg$. Then we can use the symmetry 
condition (46) on (1/3, ria, ttf) to show that the above is 
reduced to 

w°(g) = tfw{g)u = £ JJ hHggilK 1 -^ 1 . (GIO) 

{9,} i 

Note that W°(g) acts on |$o)- If we choose a new \^f' ) 
that is symmetric under the symmetry G&: ^"(flOl^o) = 
|*q), then the resulting |*') = U\^ ) will be symmet- 
ric under W(g): W(g)\W) = |*')- Since W°(g) = 
Tli \9i)(99i\K 2 " nas a simple form, it is easy to con- 
struct the deformed |$o) that has the same symmetry 
under W°(g) = Yli lftX33i|- Then after the fermion LU 
transformation U, we can obtain generic SPT states that 
are in the same phase as | \P) . 

Here we would like to stress that, only on a system 
without boundary, the total symmetry transformation 
W is mapped into a simple on-site symmetry W° the 
LU transformation U . For system with boundary, under 
the LU transformation U, the total symmetry transfor- 
mation W will be mapped into a complicated symmetry 
transformation which do not have an on-site form on the 
boundary. 

3. Entanglement density matrix 

The non-trivialness of the SPT states is in their 
symmetry-protected gapless boundary excitations. To 
study the gapless boundary excitations, in this section, 
we are going to study entanglement density matrix pe 
and its entanglement Hamiltonian He- Pe = g~ Hb , from 
the above constructed ground state wave functions for 
SPT states. The entanglement Hamiltonian He can be 
view as the effective Hamiltonian for the gapless bound- 
ary excitations. 

To calculate the entanglement density matrix, we first 
cut the system in to two halves along a horizontal line. 



To do cutting, we first split each site on the cutting line 
into two sites: \gi) — > \gi) ® \gi>) (see Fig. G4a). The 
ground state lives in the subspace where 3, = g^. We 
then cut between the splitted sites. 

Because the entanglement spectrum do not change if 
we perform LU transformations within each half of the 
system, we can use the LU transformations to deform the 
lattice into a simpler one in Fig. G 4b, where the sites are 
labeled by 1, 2, .., L, 1', 2', .., L' , along the edge which 
forms a ring, and the only two interior sites are labeled 
by a and a'. 

On the deformed lattice, the ground state wave func- 
tion can be written as U\$o), where 



l*0> - 10a) (4>a\ ® \<Pa')(<t>a'\ ®£=1 \<Pi){<t>i\ ®iL V \<t>i>){<l>ii\ 

(Gil) 



a' 




(a) a (b) 



FIG. 17: (Color online) (a) We divide the system into two 
half alone a line by doubling the sites on the line \gi) — > 
Ifli) ® |ft) and then cut between the pairs, (b) The interior of 
each region is deformed into a simpler lattice without affecting 
the spectrum of the entanglement density matrix. 
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and 



U = JJ v 3 1 (50 1 Si' 1 Sj' ! 5a' 



A 



n 

ij 

n 

i'a' 

n 



,«2 (go ,g; ,gj ) ^112 (go ,9i' ,Qji ) j-j _n 2 (go ,Sa J12 (go ,g a ,gi ) 



(Oij) C (Oi']>) 



(Oa 



'(Oa 



jij (go ,g;' ,g a ' ) "2 (go ,g;' ,g a ' ) 

C (Oi'a') L (Oi'a') A 



t«2 (ga ,gi ,gj ) t«2 (go ,ga ,g, ) -t"2 (go ,g; ,g 3 ) -t«,2 (go ,ga ,gi ) 

C (a«) C (Oaj) C (Oij) C (Oai) 



nt"2 (go ,g ; < ,gj' ) t«2 (go ,gj< ,g a < ) _fn 2 (go ,g ; < ,g a ' ) -t"a (ffi' ,9j' .9a' ) 
c Cni'oM c Cn.,'„n ( '/ , rvi/„/'i 



-(Oi'f) u (Oj'a') "((MV) ""(i'j'a') 



A 



Now, let us rewrite U as U = U up Udown, where U up acts 
on the upper half and J7down acts on the lower half of the 
system: 



(G13) 



"(0,t'+l',t') 



n-l/ \ TT n2(go.gi/ + i',ffj') .. 

^3 I9 , 0,fl , i',5j',5 , o')ll c '"---'-^ 

A 

c m . ''in;',,- 1 



J-i 2 (go ,9,' .So' ) ™2 (go ,Si' ,9a' ) . 
'(Oi'a/) c (0i'a') 



nt"2(»oifli',%/) t™2(go,g j ',g a ')-t™2(go,gi',g a ')-t™2(gi',g : ,',g a ') 
C (0i'.j') C ((W) C (0i'a') C (*'j'a') 

A 

n^(po,5 0) ^ ) ft)n 5 Fof + Mr ,ff<)x 

V » 

n 



-n2(go,g a ,gi) n 2 (g ,g a -gi) 
C (0ai) L (0ai) A 



n 



t"2 (ga ,9. ,gj ) t«2 (go ,g a ,gj ) -t«2 (go ,gt ,g, ) -t«2 (go ,ga,g>) 

C (aij) C (0aj) C (0ij) C (0ai) ' 



where the pure sign factor (— )/(si''—>9i,') arises 
from rewriting F| ■ cj 



■p-2(.go,9i,9j) n 2(go,gi> ,gj>) 

c frw<„-''i 



*J (Oii) u (0i'j') 



as 



and we have 



n nalso.ft'+i'iSc) T-r na(ffo,ffi+iifll) 
i' L (0,i' + l',i') 1-U L (0,i+l,i) ' 

used the fact that = g^ in the ground state subspace. 
Thus the express eqn. (G13) is valid only when U acts 
within the ground state subspace. So we find that 

t>u P = (-) /(Sl "-^' ) n^ 1 (9o,ft',5,',3a')x (G14) 

A 



n 



"2(go,gi' + i',gi') tt -n 2 (g ,gii ,g a ') "2(30,9*' ,g a r) 



C (0,i' + 1'A') 



n 



C (0i'a'j C (0i'a') 



n 

A 



t"2 (go, gi',gj>) tn2(go,gj',g a ')_tn 2 (go,g i ',g ')-t™2(g i ',g j ',g a ') 

^(Ol'j'^ C (n-\'n'\ ^(Ci'n'\ ^(i'i'n' 



(Ofa') °(0»'a') u (»'i'a') 



J7dow„ = I] ^O, 5a, ft, 9j) II ^ f + °M) +1 ' 9,) X (G15) 
V * 

n-n 2 (go,9a,9i) n 2 (ga,g*,gi) v 
C (0ai) C (0ai) 

n« 



t«2 (ga ,g> ,flj ) t"2 (go ,g a ,9j ) -t«2 (go, g« ,9j ) -t"2 (go, 9a ,g») 

C (aii) L (0aj) L (0ij) L (0aj) 



Now the entanglement density matrix pe can be writ- 
ten as 



pE = Tr nppet U\$ )(* \fr 

= Tr uppcr <7 up <7down I $ ) ($ 1 f>i own ^ P 



Tr uppcr C> down I $0) (*o I t^L , 



(G16) 



where Tr upper is the trace over the degrees of freedom on 
the upper half of the system. Since gi = g^ along the 
edge in the ground state, pe (as an operator) does not 
change Qi along the edge (ie pe is diagonal in the |<?j) 
basis). So it is sufficient to discuss in the subspace of 
fixed \gi)'s: 



Pe(9i, -,9l) 
ga,g' a v 

L 

c 



n-"2 (go 1 ,9i ) TT -»2 (go ,9a ,9i ) "2 (go ,9a ,9i ) v /pi 1 7 \ 

C (0,i+1,») ll C (0ai) C (0a 4 ) X l^ 1 '^ 

z— 1 az 



n 

V 

i5a)( 5 ;in 



t"2 (ga ,9i ,Sj ) t«2 (go ,9a ,9j ) -t™2 (.90 ,9i ,9j )--\n 2 (g ,g a ,gi) 
C (aij) C (0aj) C (0ij) C (0ai) 



' I TT c n 2(9o,9 a ,gi) -n 2 (g ,gi,gj) c n 2 (g ,g' a ,gj) c n 2 (g' a ,gi,9j) x 



'(Oai) u (0ij) °(0aj) , '(aij) 



n 



c t«2(go,gl,gi)gt«2(go,g a ,gi) TT gt«2(go,g»+i,g») 



'(Oai) l '(0ai) 11 (0,'i+l,i) 

i=L 



We see that pe has a form 



{g.} 

= X] l{fi}edge)<{5i}edge| (G18) 

{gi} 



where 



|{.9 4 }cdgc) = |G fc | 1/2 ^]J^3(5o,5a,5 l ,5j): 



ga v 



n-«2(go,gi+i,g;) TT -n 2 (go,ga,gi) n 2 (g ,ga,gi) 
C (0,i+l,i) 11 c '"~^ c '"--^ 



'(On 



'(Oa 



!=1 



„t"2 (g a ,gi+i ,g*) t"2 (go ,g a ,g0 ^"2 (go ,gi+i ,9.) -t"2 (go ,9a ,9«+i ) , 

"(o,i+l,i) (0,a,i) L (0 : i+1,«) C (0,a,i+1) 

\9a) ® (G19) 



We can simplify the above expression (G19). First we 
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note that 

L 

n-«2(go,3i+i,3i) TT 
C (0,i+M) 11 



-n2(go,9i+i,gi) TT ^n 2 (so ,g a ,9i) n2(go,3a,s») 



'(Ooi) u (0ai) 



L 

n 



^"2 (SO, ffa, Si) "2 (SO, Sa, Si) -"2(S0, ffi + l,9i) 

L (0ai) L (0ai) ""(0,1+1, i) 



_ ^_-)"2(so,Sa,Sl)+n2(so,Sa,Sl) Z)f =1 ™2 (so , Si+ 1 ,Si ) x 
L 



r «2(so,3a,Si) -«2(so,3i+l,Si)-«2(so,3a,3i + l) 

Therefore. 

L 

n-"2(so,Si+i,Si) TT 
C ((M+M) 11 



n "2 1,30 ,Sa ,Si J /'2 (.SO ,Si + l ,Si j -«2 (.SO ,Sa ,Si + l J 
L (0ai) C (0,J+l,i) L (0,a,i+1) l^ZUJ 



-"2(so,Si+i,Si) TT -"2(so,Sa,Si) c n 2 (so,Sa,Si) x 



"(Oai) °(0oi) 



i=l 



n 



t«2 (Sa, Si + l, Si) t"2 (SO, Sa, Si) -t"2 (SO, Si + l,Si)^t«2 (SO, Sa,Si+l) 

(a,i+l,i) (0,a,i) L (0,i+M) C (0,a,i+1) 



(-) 



n2(S0,Sa,Sl)+"2(S0,So,Sl) ]C*=1 ™2(so,3i+l,3i) . 



n 



«2(S0,Sa,Si)-n2(S0,Si+l,Si)-"2(S0,Sa,Si + l) /pni N 

C (0ai) C (CM+M) C (0,a,i+l) X l^ ZJ J 



t«2(Sa,Si + l,Si) t"2 (SO, Sa, Si) -t"2 (SO, Si + l, Si) -t«2 (SO, So, Si + l) 

L (a,i+l,i) (O.a.i) C (0,i+l,i) C (0,a,i+1) 

_ ^_yi2(S0,Sa,Sl)+"2(S0,Sa,Sl) ]Cf=l ™2 (SO ,Si + 1 ,Si ) x 

f_\Y*i=i «2(so,Sa,Si)n2(so,Sa,Si+l)+n2(so,Si + l,Si)n2(Sa,Si+l,Si) x 



2 = 1 



J"2(Sa,Si+l,Si) 

"(a,i+l,i) 



We find 

|{ft}cd g o) = |G fc | _1/2 ^]J^3(go,5a,g l +i,g l )x 

Sa I 

(__)™2(S0,Sa,Sl)+"2(S0,Sa,Sl) Z]f =1 "2 (SO ,3i + 1 ,3i ) x 
f_)Z)f=l n2(S0,Sa,Si)n2(S0,Sa,Si + l)+™2(S0,Si+l,Si)«2(So,Si+l,Si) x 

L 



n 



J"2(Sa,Si+l,Si) 

'(o,i+l,i) 



IfiO <8 (Oilfli)) 



(G22) 



From the above calculation, and using the relation be- 
tween the entanglement density matrix and gapless edge 
excitations, we learn two things. 1) The low energy edge 
degrees of freedom are label by {51, on the bound- 

ary, since Pe(9i, ■■■■,9l) has one and only one non-zero 
eigenvalue for each {g±, ...,g L }. 2) The low energy edge 
degrees of freedom are entangled with the bulk degrees of 
freedom, since the states on the site-a (in the bulk) are 
different for different {<7i, <7l}. Using the expression 
(G22), we calculation the low energy effective Hamilto- 
nian H e f{ on the edge from a physical Hamiltonian .Hedge 
on the edge: 

i H es){g' l },{g i } = ({pledge (Hedge I {#i}edge) • (G23) 

4. Low energy effective edge Hamiltonian 

Now let us assume that .Hedge only changes g k and only 
contains fermion operators C( a .fc+i,fe) and C( a .k.k-i)- We 
find 



( H cf{){g'J,{g t } = \G b \ (ilA^) 1 ^- 90 ' Sa > 9k + 1 > 9k) v 3(90, 9a, fffc , 3fc- 1 ) ^3 (^0 , 5a , fffe+1 , Sfe)^ (ffO , ffa , 5fc , fffc- 1 ) X 

Sa i=£k 

^_y»2(so,Sa,Sl)[n2(so,Sfc + l,Sfc)+"2(so,Sfc,Sfc-l)+»i2(so,Sfc + l,Sfc)+"2(so,Sfc,Sfc-l)] x 

(_y>2(3o,3a,3*,)n2(3o,3a,3fc + i)+™2(3o,3fc + i,3k)«2(3a,Sfc + l,3fc)^_yi2(^^^^ x 
^_y»2(so,Sa,Sfc-i)n2(so,Sa,Sfe)+n2(so,Sfc,Sfc-l)"2(Sa,Sfc,Sfc-i) ^_^™2(so,Sa,Sfc-i)n2(so,Sa,Sfc)+«2(so,Sfc,Sfc-l)™2(Sa,Sfc,Sfc-i) x 
"™2(Sa,Sfc,S/ t -i)+n2(Sa,Sfc+i,Sfc)+«2(Sa,Sfc,Sfc-i)+n2(Sa,SA ; + i,Sfc)=oven x 

/ J I "2(Sa,Sfc + l,Sfc) «2(Sa,Sfc,Sfc-l) TT t "2 (Sa ,Sfe ,3fc - 1 ) t™2 (Sa ,Sfc + 1 ,3fc ) I \ ( C A\ 

\yfcl C (a,fe+l,fc) L (a,fc,fc-1) ^ cd g cL (a,fc,fc-l) L (a,fe+l,fe) If*/' V^^J 



Since 



n 2 (ga,gk,gk-i) + n 2 (g a ,g k+ll g k ) + n 2 (g a ,g k ,gk-i) + n2(g a ,g k+1 ,g k ) 

= n 2 (go,gk,9k~i) + n 2 (g Q ,g k+1 ,g k ) + n 2 (g ,g' k ,g k ^i) + n 2 (g , g k+1 , g' k ) mod 2, 



(G25) 
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the factor (— ) n ^(9o,ga,gi)ln2(go,gk+i,gk)+n2(.go,gk,gk-i)+n2(go,gk+i,g'k)+n 2 (go,g' k ,gk-i)] can dropped: 
( H cs){ g 'j,{ 9l } = IGbT 1 ^ ( Il^i'sO u * (.9o, 9a, gk+i,g' k )^{g , g a , g' k ,gk-i)v3 (go, 9a,. gk+i, 9k)v 3 (go, 9a, gk,9k-i)x 

g a i^k 

/_\n2(go,ga,gk-i)n2(go,ga,gk)+n2(go,gk,gk-i)n2(.g a ,gk,gk-i) f_\n2(go,ga,gk-i)n2(go,ga,g'k)+n2(go,g'k,gk-i)n2(ga,g'k,gk-i) x 
^■n2{ga-gk-gk-i)+n2(ga,gk+i-gk)+n. 2 {ga,g' k -gk-i)+n2(ga,gk+i-g' k )=cvcn x 

IJ I n 2 (g a ,gk+i,g'k) n2(g a ,g' h ,g k -\) tt Jn 2 (ga,g k ,gk-i) Jn 2 (g a ,gk+i,gk)i \ (nofi\ 
\9k\ c (a,k+l,k) C (a,k,k-1) jned S cC (a,fc,fe-l) C {a,k+l,k) \9k)- 

Now let us consider that following three form of H cx \ gc 

TtOO _ TtOQ TT + + _ TT + +J J TT + - _ fH J „ (r07\ 

-"edge--" . U c dgc- U C {a,k+,k) C {a,k,k-iy -"edge ~ C (a,fc+l,fc) C OAfc-l) ' l^'J 

where H oa , H ++ , and are bosonic operators that only act on \gu)- The corresponding H c q are given by 

( H cff) Wj,{g.} = IGbT 1 ^ {]^ s g' v gi)^Ugo, g a , gk+i, g' k )^(go, g a , g' k , gk-i)^3(go, g a , gk+i, gk)Mgo, g a , gk, gk-i)^< 

g a i^k 

(_)™2(go,3a,3fc)n2(gO:5a,5fc + i)+«2(5o,3fc + l:5fc)™2(ga^ x 

^_^™2(S0,9a,9fc-l)™2(S0,9a,9fc) + "2(S0,Sfc,Sfc-l)n2(Sa,Sfc,3fc-l) ^_^™2(ff0,9a,Sfc-l)»l2(50,9a,s! c )+™2(90,Sfc,Sfc-l)™2(9a,Sfc,Sfc-l) x 

^n2(ga,gk,gk-i),n2(ga,g' k ,g k -i)^n2(g a ,g k +i,g k ),n2(ga,g k +i,g' k )-^g' k ,g k i (G28) 

^_y»2(ff0,9a,gfc)n2(g0,9a,fffc + l) ^ „ "J "2 (SO , Sa ,9fc ) » l 2 (SO :9a ,Sfc + 1 ) +« 2 (ff , Sfc + 1 ,9fc ) x 

^_y n 2(ffo,9a,gfc-i)n2(so,ga,gfc)^_^n2(9o,ga,9fc-i)n2(9o,9a,9fc)+n2(so,Sfc,Sfe-l) x 

^ri2(ga,g k ,g k -i),0^n2(ga,g k +i,gk),0^n2(ga,g' k ,gk-i),^n2(ga,gk+i,g' k ),l-^g' k ,g k > (G29) 
(- ff eff _ ){9a.{Si} = l Gb l _1 X! (n^9i) y 3XSO>Sa,0fc+l,SjD^3(SO,^^ 

(•_^™2(ffO,Sa,9fc)n2(gO,Sa,fffc + l) ( „ ) «2 (90 . 9a , 9fc ) "2 (90 :9a ,9fc + 1 ) +«2 (SO , 9fc + 1 ,9fc ) x 

^_^™2(90,9a,9fc-l)n2(90:9a:9)c)+™2(90,9fc,9fc-l) ^_^"2(go,9a,gfc-l)ri2(go,9a,9fc) x 

^n2(ga,g k ,g k -i),l^n2(ga,g k +i,g k ),0^ii2(ga,g' k ,g k -i),0^n 2 (g a ,g k+ i,g' k ),lH g ^ gk . (G30) 

Those expressions allow us to constructed simple edge effective Hamiltonians for the corresponding SPT phases. 

To obtain more generic edge effective theory, let us assume that -Hedge changes both and gk+2 and contains 
fermion operators C( aife+1)fe ), C( 0)fc)fc _i), C( aife+3)fe+2 ), and C( a , fe+2 ,fc + i): 

9a 

v t (go , g a , gk+3 , g' k +2 ) ^ (go , g a , g' k +i , 9k+i ) ^ (go,g a ,gk+3,gk+2) »3 (go ,g a ,gk+2,9k+i)x 

f_\n2{go,ga,gk)n2(g ,g a ,gk+i)+n2(go,gk+i,gk)n2(ga,gk+i,gk)/_\n2(go,ga,gk) n ^ x 
f_\n2(go,ga,gk-i)n2(go,ga,gk)+n2(,go,gk,gk-i)n2(,g a ,gk,9k-i) {_y*2{go,ga,gk-i)n2{go,ga,gk)+™v(go,gk>9k-i)™2(ga,g'k,gk-i) x 

^_^«2(S0,Sa,Sfc + 2)n2(S0,ff<i,Sfc + 3)+«2(S0,Sfe + 3>Sfc + 2)«2(Sa,fffc + 3,Sfe + 2) y« 2 (90 ,9a ,9fc + 2 )«2 (90 ,9a ,9fc + 3 ) + «2 (90 ,9fc + 3 ,9fc + 2 )™2 (9a ,9fc + 3 ,9^ + 2 ) X 

f\»2(9o,9a>9fc-i)n2(9o,9a>9fc)+™2(90i9fc,9fe-i)™2(9a,9*^ x 

^n2(9a,9fc,9k-l)+n2(go,9k+l>9fe)+n2(9a,9£,9k-l)+«2^ 

/„' I n 2(9a,9k + 3,9fc + 2 ) "2(9a,9fc + 2'9fc + l) "2 (ffa ,9fc + l ,9fc ) "2 (ffa ,9fc ,9fc - 1 ) fT w 

\yfcl C (a,fc+3,fe+2) L (a : fe+2,fe+l) L (a,fe+l,fe) L (a,fc,fc-1) ^edge* 

t«2(9a,9fc,9fc-l) t«2(9a,gfc + l,9fc) t "2 (ga ,9fc + 2 , 9fc + 1 ) t ™2 (fia , 9fc + 3 ,9fe + 2 ) I \ fP^I ~\ 

C (a,fe,fc-1) t "(a,fc+l,fc) C (a,fc+2,fc+l) C (a,fe+3,fe+2) If*!/- V J ' ix ) 



r 



Appendix H: Ideal Hamiltonian tonian 

After constructed the fermionic SPT state |\&) in the 
Section VI A, here we would like construct a local Hamil- 



H = -Y i H i (HI) 
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FIG. 18: (a) The matrix elements of the Hamiltonian 
term Hi can be obtained from the fermionic path inte- 
gral of the fixed-point action amplitude V3 on the complex 
(i?i92ff3ff435i765iffi}- The branching structure on the complex 
is chosen to march that of the trianglar lattice in Fig. 10. 
(b) The action of Hi can be obtained by attatch the complex 
in (a) to the trianglar lattice in Fig. 10. The Grassmann 
numbers on the six overlapping triangles are integerated out. 



that has the Gf symmetry, such that is the unique 
ground state of the Hamiltonian. 

We note that the path integral of the fixed-point action 
amplitude not only gives rise to an ideal ground state 
wave function (as discussed above), it also gives rise to 
an ideal Hamiltonian. From the structure of the fixed- 
point path integral, we find that Hi has structure that 
it acts on a seven-spin cluster labeled by i, 1 - 6 in blue 
in Fig. 10, and on the six-fermion cluster on the six 
triangles inside the hexagon (123456) in Fig. 10. 

The matrix elements of Hi can be obtained by evaluat- 
ing the fixed-point action amplitude on the complex 
in Fig. 18a, since the action of Hi can be realized by 
attaching the complex Fig. 18a to the triangular lattice 
and then performing the fermionic path integral on the 
new complex (see Fig. 18b). The evaluation of V 3 on 
the complex give us 



\m, 1 (g i ,g' i )+m 1 (g' i ,g 1 )+m 1 (g i ,g 2 )+m 1 (g 3 ,g i )+m 



nv 3 



u(ff5,9i)+TOl(3;,96) / J~J d9 U2 } 

" A IOC 



"2 (gi,g'i,9j) j H«2 (gt ,g' t ,gj ) 
3') 



3=1,2,6 

mi(g4,gi)+m 1 (g 5 ,g' i )+m 1 (gi,g 6 )+m 1 (g' i ,g 1 )+m 1 (gi,g 2 )+m 1 (g3,g' i ) f\rn 1 (g i ,g' i ) x 

V 3 + (.95, 9i, g'i, 9eWt (54, 55, 9i, 5i)V 3 "(54, 5a, 9h 9i) V 3 (93, 9i, g[, ff2)V 3 " {gi,g'i,g2,gi)V${gi,g'i, g 6 ,gi 

^m 1 (g4,g' i )+m 1 (g' i ,g 1 )+m 1 (gi,g 2 )+m 1 (g 3 ,gi)+m 1 (g 5 ,gi)+m 1 (gi,g 6 ) 

3=1,2,6 

^rn 1 (g 4 ,gi)+m 1 (g 5 ,g' i )+m 1 (gi,g 6 )+m 1 (g' i ,g 1 )+m 1 (gi,g 2 )+m 1 (g 3 ,g' i ) ^_^m 1 (g i ,g' i ) x 



11 (j« ) (3" ) 

3=3,4,5 



/ n < 

" A IOC 



"2 (si ,g\ ,gj ) j Hn 2 (ffi ,g'i -.gj ) 
3) (H'i) 



nj o«2 (gj,gi ,g'i )jn n 2 (gj ,g> 

3=3,4,5 



^3(95, 9i, 9i, 96)0 (u> 



"2 (fli ,g'i,9e) n"2 (s5 ,gi,gn) /)«2 (95 ,gi ,ge) fl «2 (gs ,g; ,g'i ) , 



./ <l fl™2 (35 ,9. ,g'i ) fl«2 (S4 ,95 ,9^ ) /)"2 (94 ,9i ,9i ) /1«2 (94 ,95 ,9i ) 



("/i ^ n n '\ff l 2yy^-gi-gi!Q ll 2Kgi-g&^i)ff i 2\gi,gi,yi)ff' 

(ii'6) "(5<6) u (5i'6) v (5ii') fS) i/ii 9iP(fiU') {7 (45i') I7 (4ii') P (45i) 

,/— ^, ^ „'\/j™2(94,93,9i)/i II 2(94,gi,g-)fl n 2(94,g3,9i)fl n 2(g3,9i,9i),,-i/„ „ _/ „ \ a n 2 (g 3 ,gi,g'i) a n 2 (g 3 ,g[,g 2 )n n 2{g 3 ,gi,g 2 )n n 2(gi,g'i,g2 

V Z (34,53,^,^)^(43^ (4ii ,) f^j,) f 3 (93,^,^,52)^(3^,) (3i , 2) «( 3l2 ) ^-2) 

..-l/'^ ^ „ Afl™2(9i,9i,g2) a n 2 (gi,g 2 ,gi)^n 2 (g i ,g' i ,g 1 )^n 2 (g' i ,g 2 ,g 1 ) , / ^ „n 2 (g( ,ge,9l ) /i™2 (g, ,g- ,91 ) fl»2 (g« ,96 ,gi ) J»2 (g; ,sj ,ge) 

^ 3 [9i,9i,92,9l)0( W 2) V( i2 l) V(U>1) V21) ^3(5i,5i,56,5l)^ 6 i) ^1) %61) y («'6) 



'(ii'6) 

(H2) 



Just like the wave function, we have include an extra factor (-)™i(9i,90+»H(gl,gi)(_)m 1 (g l ,g 2 )+m 1 (93,g l )+™i(g5,g,)+m 1 (g l ,g 6 ) 
in the path integral. With this extra factor, multiplying eqn. (H2) to the wave function can then be 
viewed as attaching the complex Fig. 18a to the triangular lattice (see Fig. 18b). Note that the factor 

^rn 1 (g i ,g 2 )+m 1 (g 3 ,gi)+m 1 (g 5 ,g i )+m 1 (g i ,g e ) plus the factor that we mc l u d e d in the wave function (-)™i(fl4,fl*)+"»i(ff*,fli) 

give us a factor (-) m i(s4,g ! )+™i(9 ! ,gi)+™i(g s ,g2)+mi(g3,g I )+™i(95,g 1 )+™i(g.,9 6 ) on a n t h e six interior edges {i, 1}, {i,2},..., 
{i,6}. This factor is needed in the path integral on the complex in Fig. 18b obtained by attaching the complex 
in Fig. 18a to the triangular lattice. Also note that the factor (-) m i(94'9*)+ m i(s<>9i) j s w h a t we neec i j n the 
new wave function after the action of attaching the complex Fig. 18a to the triangular lattice (see Fig. 18b). 
By combining all the mi factors on the right-hand-side of eqn. (H2), we find that all those mi factors become 
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f—\n 2 (g3,gi,g' i )+n2(gi,gi,gi)+n 2 (g5,gi,g'i)_ g we fi nc J 

^_yn 1 (g 4 ,g' i )+m 1 (g' i ,g 1 ) ^_^m 1 (gi,g2)+m 1 (g 3 ,gi)+m 1 (g s ,gi)+m 1 (gi,g 6 ) / J~J (H3) 

_ C V l2 (93,3.,ffD+«2(94,S,:SD+™2(95,Si,sD / 1 r f j/)»2(fl<jSi.fli)jfln2(ai,flJ,flj) 1 f r] /in2(ffi,fli,Si) r lS«2(ffj 

- I J / 11 at ^'j) at/ (ii'j) 11 a %w) aw (i«') 

•* J. IOC „■ Q/IC 



s 

ui g"2 (gj,gi,g'i) A ^n 2 (g, ,gu9i) 

3=1,2,6 3=3,4,5 
^_y i 2(94,g5,Si)n2(S4,ffi,ffO+ n 2(94,9i,sO™2(ff4,S3,gO+"2(ff^92,gi)n2(gi,S^Sl)+n2(gi,9^ffl)n2(s-,S6,gi) 

f„ „ „ \ a n 2 (gi,g' i ,ge) a n 2 (g-i,gi,ge)n n 2(g5,gi,g6)nn2(g5,gi,gi),, i„ „ „ „/\a n 2{95,9i,gi) An 2 (94.,g5,g'i) a™2(s4, ss,9i) n n 2(gi,gi,g'i) 
V3{95,9i,9i,96)V( U , 6) V (5i6) f/ (5i , 6) y (5ii ,) V3{94,95,9i,9i)0( 5W ) f (4S i') 6 '(45i) y (4«') 

„ „ „r^ a n2(g4.,gi,g' i ) a n2(g4.,g3,gi)n n 2(g4,g3,g'i)n n 2(g3^i,g' i )_ : -if n / s/)«2 (93,94,34) a n 2 (g3,g'i,g2) ^"2(33, 3^,92) ;p2(g4,9i,92) 

^ 3 \,9i,93,9i,9iW {iiil) v (m) (m >) v s {93,gi,gi,g2)v (3ii , ) f (3i , 2) tf (3l2) f (ii , 2) 

, .— 1 z'- _/ _ _ \ fl"2 (g<, g 92 ) /)n 2 (gi, 92, 91 )n™2 (9 i,g2, 91 )n n 2(gi,g'i ,91 ) ( „ J „ „ ^/) n 2(gi,g-,gi) /^(g-.ge.gi) nn.2{gi,g&,gi)n n 2(gi,g'i,gf>) 
v z [9i,9i,92,gi)0 (ii , 2) C (i21) (i , 2 i) V3K9i,9i,96,9i)0 {ull) (m) f (i61) ( ^, 6) 

In the above, we have also rearranged the order of the Grassmann numbers within some simplices to bring, say, 
gnvfa^g'i) next t Wff4,ffi,fli) N integral out d(9d#'s and obtain 

(_y™i(94,g0+ m i(ffi'9i)^_) m i(ffi.92)+"il(93,9i)+mi(95,9,)+m 1 (9 i ,g 6 ) / TT Vg = (_)™2 (9i ,gj ,9e) x (H4) 

(_y i 2(g3,gi,gO+"2(g4,gi,gD+™2(g5,g;,gO+™2(g4,g5,g;)n2(g4,9i,gO+"2(g4,g^ 

r« ^, „' o \n n 2(g5,gi,ge)n n 2(g5,g' i ,g6) ll / „/\/)™2(g4,g5,gj)a«2(g4,g5,g4) -1/ „/\a n 2{94,,93,9i) A«2(g4,g3,g0 

V3{g5,gi,g i ,g6)v {5i6) f (5i , 6) m9^g5,gi,g i )V( 45i , ) f/ (45i) f 3 (54,53,5i,5d0( 43i ) 

i,- 1 (n n n' n ^"2 (93 ,9- ,g 2 ) ^"2 (g 3 ,9> ,92) -1 / / \ fl ii2 (g< ,92 ,9l ) n"2 (94,92 ,91 ) / / \ /i™2 (gj ,96 ,91 ) 5"2 (g, ,96 ,91 ) 

^3 {93,gi,gi,g2)V( 3i , 2 ) y ( 3 i2) ^3 w J 5i,52,5ij w (i 2 i) y (i'2i) v3{gi,gi,ge,gi)V( V61 ) ^ (i6 i) 

where the factor (— ) Tl 2(gi,g i ,ge) come s from bringing 0?. 2 ,^' Si ' 96 ^ in eqn. (H3) all the way from the back to the front. 



(ii'6) 

Let us rearrange the order of the Grassmann numbers to bring, say, ^g 2 j/g)' S< m front of 97^\' 9% " "'' 



(_)™ l i(g4,gD+ m i(9i,9l) ^_^mi(9i,92)+mi(93,9i)+mi(g 5 ,gi)+mi(9i,96) / J~J y g 



(_^ n 2(gi,g^g6)+n2(g3,gi,g-)+™2(g4,gi,g-)+ rl 2(g5,gi,g-) x (H5) 
(„) n 2(g4,g5,gi)n2(g4,gi,g0+™2(g4,9i,9i) rl 2(g4,g3,gi)+n2(g-,g2,gi)n2(gi,g' i ,gi)+n2(gi,9i,gi)n2(g-,g6,gi) x 

^3(g4,g5,g»,gQ^3(ff5,g»,gj,g6)^3(g»,g^g6,gi) , N«2(g5,g, , ,g 6 )"2(g5,gi,g6)+"2(g4,g3,g, , )»2(g4,g3,gi)+»2(g, , ,92,gi)"2(gi,g2,gi) x 

^3 (g», 5^ 52, 51)^3 (53, 5*, 5i, 52)^3(54, 53, 5*, 5*) 

ff^2 (9B,9j ,96 ) /)"2 (95 ,94 ,96) /i™2 (94 ,95 ,9j ) /p2 (94 ,95 ,9; ) n«2 (94 ,93 ,9; ) /jri2 (94 ,93 ,94 ) 
(7 (5i'6) P (5i6) I7 (45i') (7 (45i) l7 (43i') p (43i) X 

/3"2 (93 ,9j ,92) /)n 2 (93 ,94 ,92 ) n"2 (94,92,91) /)«2 (g« ,92 ,9l) /i«2 (9i ,96 ,91 ) /)"2 (g»,96 ,9l) 
P (3i'2) t7 (3i2) P (i'21) P (i21) (7 (i'61) l7 (i61) 



Indeed, the above expression can be regarded as the fermion coherent state representation of the ideal Hamiltonian. 
We can replace, for example, dffi/^tfgft'"'^ as: 

c ( 5l6) J c "m)' 9 - (1 -^2 (55,54,56)) (1 - "2 (55,54,56)) c[ 5i6) C( 5i6) , (H6) 

We note that the fermion coherent state is defined as |0(5»6)) = |0) — ^(5i6) c ( 5 j 6 )|0)- ^ i s eas Y to check 

(^s^)!^)"^ 5 ' 9 ^ 6 ^^) 5 '"' 365 - (! - ^(55,5^,56)) (1 - n 2 (55,5.,56))cj 5i6) c (5i6) |% 6) ) = ^^fl^*^ 

(H7) 

However, a sign factor (— ) n 2(95,gi,ge) j s required when the Hamiltonian acts on the coherent state since 
de M9,,9i,ge) d ^2(g,, gi ,ge) ghould be reordered as diS^' Si ' 9 ' ) d9^' 9i ' B ' ) . Similarly, for the face (45i), the fermion 

coherent state is defined as |^( 4 5j)) = |0) — ^(45i)c| 45i J0). In this case d^^ 2 ^ 4,95 ' 9 *^^ 2 ^ 4 ' 95 ' 3 ^ is the correct ordering 
and we don't need to introduce the extra sign factor. 
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By applying the same discussions to other triangles, we find that 43i and H2 also contribute sign factors 
(_)"2(94,s3,9i) an d (_)n2(si,s2,si). We also note that 



and finally obtain 



^_y»2(9i,g-,96)+n2(S3,Si,30+™2(S4,S»,9-)+™2(S5,Si,S-)(_)™2(S5,Si,96)+n2(S4,S3,Si)+n2(Si, 92, 91 ) 
^_y i 2(S4,93,90+™2(ss,S-,96)+n2(s»,S2,9l) 

^3(54, 55, 9i, 91)^3(95, 9i, g'i, 9a)^3 (5i, g' t , 56, 5i) „ 



(H8) 



H *= 1^ 5*, 5i525354555e) 5*, 515253545556 —7 ; r - ? 

^(5^,5 ,52, 51)^3(53, 5i, 5 ,52)^3 54,53, 



9* ,9; 



5i,5 



(H9) 



(i61) 



«2(95,9i,96) 

C (5i6) 


-(1 


- "2 (55 


5i,5e)) (1 


"2(94, 95, 9i) 
C (45i) 


-(1 


- «2(54 


55,50) (1 


"2(94, 93, 9i) 
C (43*) 


-(1 


- n 2 (54 


53,50) (1 


«2(93,9i,92) 

C (3i2) 


-(1 


- "2(53 


5O52)) (1 


"2(9i,92,9l) 
C (i21) 


-(1 


-n 2 (5i,52,5i))(l 


"2(9i,96,9l) 
C (i61) 


-(1 


- n 2 (g' i7 


56, 5i)) (1 



^_yi2(94,95,9i)™2(94,9i,9i)+™2(94,9i,9i)"2(94,93,9i)+n2(9i,92,9l)™2(9i,9i,9l)+™2(9i,9i,9l)™2(9i,96,9l) x 

(_y i 2(95,9^96)™2(95,9i,96)+«2(94,93,90"2(94,93,9i)+ri2(9i,92,9l)™2(9i,92,9l) j-_yJ2(94,93,9D+"2(95,9i,96)+n2(9i,92,9l) x 



( C (43 4 )) 
( C ( 3l 2)) 
( C ( 4 21)) 

(< 



(H10) 



The above expression is valid in the subspace where the fermion occupation number n^k on each triangle (ijk) satisfies 
riijk = n 2 (gi,gj,gk)- 

Next, we would like to show that Hi is hermitian. This property is very important and it makes the whole theory 
to be unitary. In the above, we express the matrix element H^g. ^ (515253545556) as: 



^_^m 1 (g i ,g' i )+rn l (g[,g l )+m l {g i ,g2)+m l (g 3 ,g i )+m-i(g 5 ,g i )+m l {g 

hence, its hermitian conjugate HT g , (5i5253545556) reads: 

^_^m 1 (g4,gi)+m 1 (gi,g 1 )+m 1 (g' i ,g 2 )+m 1 (g 3 ,gl)+m 1 (g 5 ,g' i )+m 1 (g' i ,g e ) / J~J 

" a in 



nv 3 , 



(Hll) 



(-) 



n2(g'i,9i,9j) iQ n z{a'i>gi>Sj) TT j/)™2(9j ,9i,9i) A n n 2(9j ,9i,9i) 
(ii'j) aU (ti'j) 11 a °(jii') aV (jii>) 

.7=1,2,6 .7=3,4,5 

m l(94,9i)+ m i(95,9i)+«il(9i,96)+"il(9i,9l)+mi(9^,92)+mi(93,9i)^_y™l(9i,9i) v 



n' r, r, ^fl™2 (95 ,9^9i) fl«2 (95 ,9i ,96 ) fl ™2 (95 ,9^96 ) fl «2 (9; ,9i ,96 ) - 1 / / 

^3 I55,5i,5i,56jy (5ii /) f (5i , 6) ^(5*6) y («'6) ^3 \9i, 95, 9i, 9i r 



"2 (94 ,95 ,9; ) /)™2 (94 ,9; ,9i ) /in2 (94 ,95 ,9i ) /)™2 (95 ,9^ ,9i) 
t7 (4ii') P (45i') ,7 (5*i') 



f„ „ „ \n n 2(93,9i,gi) a n 2 (g4,g3,9i) n n 2(94.,g'i,gi) n n 2(94,93,9i) i„ „l „ „ \ /)™2 (9 ■, 9i , 92) /)"2 (93 ,94,92) ^"2(93, Si ,92)^12(93,94, 9.) 

^3(54,53,54,54)^(3^) *(43i') ^(4™') 6 '(43 l ) I/ 3(53,5i,5i,52)6 l (li /2) y (3i2) y (3i'2) y (3ii') 

(r, 1 r, r, n \ fl™2 (9> ,92 ,91 ) ^"2 (9^9; ,9i ) fl «2 (9^92 ,91 ) fl M9^9> ,92 ) - 1 / / ■. a n2(g' i ,gi,g6)( j n2(g' i ,ge,gi)pn2(g' i ,gi,gi)pn 2 (gi,g6,9i) 
V3{9i,9i,92,9l)V(i>21) U (Wl) y (i21) y (»i'2) ^3 (ft, ft, 56, 5l]<V 6 ) f( i6 i) Vei) 

(H12) 

Note ^* = v~ Y (y is a C/(l) phase factor) and (#i6>2#3#4)* = O48382O1. However, we can not directly com- 
pare the above expression with H^g. ^i (515253545556), since its a function of new pairs of Grassmann variable 

n 2(g5,g'i,9e) ^n2(95,9»,96) /)™2 (94, 95, 9i)/5™2 (94, 95, 9») /)™2 (94 , 93, 9i)fl n 2 (94, 93, 9i) /)"2 (93, 9^,92) /p2 (93, 9i, 92) /)"2 (9i, 92, 9l)/)™2(g>, 92, 9l) 
^i'e) ^e) ' t7 (45i') P (45i) ' P (43i') P (43i) t '(3i'2) (7 (3i2) ' < 7 (»'21) t7 (i21) 

and ^™ i /gi)' 36 ' 91 '' 0^q^ ,96,9i \ Such a difference is simply because the action of maps, for example, #(5,6) to O^q) 
while the action of i/| maps #(5i'6) back to #(5^6) • Thus, to sec whether they are the same mapping in the Hilbert 
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space or not, we need to redefine 0( 5 i6) (0(5*6)) as 0( 5 i'6)(0(5i'6)) and0( 5i / 6 )(0(5i'6)) as 0(5*6) (0(5i6))- A simple way to do 
so is just replacing by in the above expression: 



(_\ni 1 (g i ,g l )+m 1 (g i: g 1 )+7n 1 (g' i:g 2)+rn 1 {g 3 ,g' i )+Tn 1 ( g5 ,g' z )+rn 1 (g' i , g6 ) f TT J a"2 (s < ,9. ,9j ) j«n2 (g-,9i ,9j ) TT j /j"2 (ffj ,9i,9i ) .S»2 (fi 

I J /ll aw (i'u) a Vii) 11 at V*) ay (i*'i 



q™2 (g j ,9i,Si ) A n n 2 fe ,9i ,9i ) 
^(ji'i) 

J=l,2,6 .7=3,4,5 
^_^mi(g4,g;)+mi(g5,g i )+mi(g^g 6 )+"il(Si,9l)+nil(9l,S2)+nil(g3:9i)^_yni(Si,Si) x 

n J n n \n n 2(95,gl,gi) a n 2 (g 5 ,gi,g e ) nn 2 (gs,g[,gs) Qn-zig'^gi.ge) -1/ / \ ,,"2 (94, 95 ,9- ) fl "2 (94 ,9; ,9i) nn 2 (gi ,g 5 ,9i) #2 (95 ,9 ■ ,fli) 

^ 3 {g5,gi,gi,g6)V i5i , i} (5i6) (5i , 6) (i , i6) i/ 3 (34, 35, ft, ft)^^ 0(«'i) ^45;) y (5;';) 

,/ f„ „ „' „ \ fl ™2(93,9i>9i)/in2(94,g3,9i)fl n 2(g4,g-,9i)fl™2(g4,93,9i),, /_ _/ - _ \/)™2(g-,gi,92) a n 2 (g 3 ,g' i ,g 2 ) n n 2(93,9i,92) n n 2(93,9i,9i) 

zMs^, 53, s^sgt^^) 0( 43i ) 0^/^ (43i /) M93,gi,gi,g2)v {i , i2) f (3i , 2) (3i2) (3i 'i) 

o\iJi: an yz> 31/ («21) {I'll) {1 '21) (i'i2) 3 (z'26) (i'6l) (I'll) (*61) 

_ ^_yni(g4,gi)+mi(g,,gi)+mi(gJ,g 2 )+m 1 (g3,gO+"il(g5.gD+ m i(gi ! 96) / TJ ^"2 (fl* ,9. ,9) ) ^Q n 2 (g't ,9 t ,9j ) J~J dO™ 2 ^, 

.7=1,2,6 .7=3,4,5 
^_^™il(g4,g'i)+mi(g5,gi)+m 1 (g' i .g 6 )+m 1 (gi,g 1 )+m 1 (g' i ,g 2 )+m 1 (g 3 ,gi) /_\mi(g<,9i) x 

V 3 ~ (55 , 5i , ft , 5e ) VjT ( 54 , g 5 , ^ , 3, ) V 3 + (34 , 33 , g[ , ft ) V 3 + (.93 , ft', 9i , 52 ) (3*, ft, ft, 51)^(3,', ft, 3 6 , 31) 

— ^_yrn(g4,gi)+rn 1 (gi,g 1 )+m 1 (g^g 2 )+m 1 (g 3 .g' i )+m 1 (g 5 ,g' i )+m 1 (g' i ,g 6 ) I J~J (H13) 



"2 (9j ,9i,9« ) Jfl^Cgj ,9j ,9*) 



We note that in the last line we evaluate the complex Fig. 18a in a different way(by choosing opposite orientation for 
the internal link) . The topological invariance of V3 path integral implies the two different ways must give out the same 
results. It is also not hard to see that the new sign factor (-) m i(a4.^)+"n(g,,gi)+n ll (gig 2 )+ mi (g 3 .g' z )+ mi (g 5 ,g' z )+rn 1 (g' l .g e ) 

is equivalent tO (— ^ m i(94.,gi)+rni(g' i ,gi)+m 1 (g i .g 2 )+m 1 (g 3 ,g i )+rn 1 (g 5 .g i )+m 1 (g i .g 6 ) ^ snlce: 

mi(ftt,ft) +mi(ft,ft) +m 1 (g / l ,g 2 ) +m 1 (g 3 ,g / l ) +m 1 (g 5 ,g / l ) + m 1 (g / l ,g e )+ 

m i(g4,g / i ) +mx{g' i ,g x ) + m 1 (g l ,g 2 ) +m 1 (g 3 ,g l ) +m 1 (g 5 ,g i ) + m 1 (g l ,g 6 ) 
='^2(35, 5*, ft) + n 2 (g 5 ,g l ,g 6 ) +"2(34, #5, ft) +n 2 (34,55,ft) + ^(fttjfts, 94) + n 2 (g 4 ,g 3) ft)+ 

"2(33, 3*, ft) +"2(53, ft, ft) +"2(3*, 32, 3i) + n 2 (g l ,g 2 ,g 1 ) + n 2 (g[, g 6 , 31) + n 2 (g h g 6 , 31) 
=0 mod 2. (H14) 



Thus, we have proved that the matrix el- 
ement ^^(.g^ftftftft) is the same as 
#* g < igi (ftft2ftsftiftft>), implying Hi is a hermitian 
operator. 

Finally, we note that the Hamiltonian only depends 
on v 3 and n 2 . It does not depend on mi. So the Hamil- 
tonian is symmetric under the Gf symmetry. By con- 
struction, the ideal ground state wave function is an 



eigenstate of Hi with eigenvalue 1. The topological in- 
variance of V3 path integral implies that Hf = Hi and 
HiHj = HjHi. Thus Hi is a hermitian projection oper- 
ator, and the set {Hi} is a set of commuting projectors. 
Therefore H = — -H* ^ s an exactly soluble Hamilto- 
nian which realizes the fermionic SPT state described by 
(v 3 ,n 2 ,u g 2 ) E & 3 [G f ,U T (l)}. 
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